Maple 16 Questions and Posts

These are Posts and Questions associated with the product, Maple 16
Hello,When I plot graphs I can simply add colors, gridlines, labels, etc., but when I want to do this with the command 'Explore' it doesn't work.The commad 'Explore' sees the words 'color', 'labeldirections' and other commands as extra variables.For exmaple is does recognize 'labels' and 'title', which I can fill in, but it doesn't some of the other commands I would like to insert.I have a formula for deriving the eigenfrequency of a concrete beam, where I want to add some things to (like the things mentioned above).Somebody knows what the difference is between 'normal' plotting and plotting with 'Explore'?Explore_-_Adding_Col.mwGreetingsFrank

Hi,

i wrote the following code

restart:
zet:=0.5:
rhop:=3880;
rhobf:=998.2;
cp:=773;
cbf:=4182;

eq1 := diff(u(eta), eta, eta)+0.133762025280e-2/(0.993000000000e-3+0.388362300000e-1*phi(eta)+.530162700000*phi(eta)^2)+(1/(eta-1)+(0.388362300000e-1+1.06032540000*phi(eta))*(diff(phi(eta), eta))/(0.993000000000e-3+0.388362300000e-1*phi(eta)+.530162700000*phi(eta)^2))*(diff(u(eta), eta));
eq2 := diff(T(eta), eta, eta)+(.63267672*((2.66666666666*(-1.1752324*10^6*phi(eta)+4.1744724*10^6))*u(eta)/p2+7.04876575828*(diff(phi(eta), eta))-(1.58058605349*(.597+4.45959*phi(eta)))*(diff(T(eta), eta))/(1-eta)))/(.597+4.45959*phi(eta));
eq3 := diff(phi(eta), eta)-5.00000000000*phi(eta)*(diff(T(eta), eta));
p:=proc(pp2) if not type([pp2],list(numeric)) then return 'procname(_passed)' end if:

res := dsolve({eq1=0,subs(p2=pp2,eq2)=0,eq3=0,u(0)=0,u(1-zet)=0,phi(0)=0.008,T(0)=0,D(T)(0)=1}, numeric);
F0,F1,F2:=op(subs(res,[u(eta),phi(eta),T(eta)]));
res(parameters=[pp2]);
evalf(2/(1-zet^2)*Int((1-eta)*(F1(eta)*rhop+(1-F1(eta))*rhobf)*( F1(eta)*rhop*cp+(1-F1(eta))*rhobf*cbf )/(F1(eta)*rhop+(1-F1(eta))*rhobf)*F0(eta),eta=0..1-zet))-pp2
end proc;

fsolve(p(pp2)=0,pp2=(20000)..(8000000));

 

it has an error and i cannot find that. would you please help me?

 

thanks

Good day, Is it possible to plot multiple (more than two) functions with 'dualaxisplot'?I've found some information about it on the Internet, which helped me out a lot, but I still don't know how to put in multiple functions in one dualaxisplot.For example I would like to plot a figure which consists of four functions, all of them having their own domain(!).Three of them are related to the 'left' y-axis, and the fourth function is related to the 'right' y-axis.Comment: I wanted to insert the formula over here, which I use to plot the dualaxisplot, but the 'dualaxisplot-command' isn't working ("You have entered an invalid Maple plot command").So for example, I do know how it works with two formulas, like:dualaxisplot(plot(2x,view=[0..10,0..20],labels=[x,left=2x],legend=2x),plot(x^2,view=[0..10,0..100],labels=[x,right=x^2],legend=x^2,color=blue)); Now I can independently change the left and right y-axis like I want to. Now I would like to have the shape of the next graph:The trhee formulas which belong together, and that I want to relate to the left axis (making one line), consist of: horziontal line, linear line, horizontal lineThe other line I want to relate to the right y-axis, which is a linear line. So I can write the four functions for example like:restart;with(plots):eq1:=50,x=0..3:eq2:=(50/(5-3))*(5-x),x=3..5:eq3:=0,x=5..10:eq4:=10x,x=0..10: I thought I could 'join' the (first) three functions by using the '{' and '}' symbols, but that isn't working either. Somebody has a suggestion?Or that it just isn't possible at all.... Thanks anyway for thinking along. Greetings, Frank   

hi

please help me to solve the following system of ODE

 

diff(u(eta), eta, eta)+.5699919300/(0.9930000000e-3+0.3883623000e-1*phi(eta)+.5301627000*phi(eta)^2)+1/(eta-1)+(0.3883623000e-1+1.060325400*phi(eta))*(diff(phi(eta), eta))*(diff(u(eta), eta))/(0.9930000000e-3+0.3883623000e-1*phi(eta)+.5301627000*phi(eta)^2) = 0;

 

diff(T(eta), eta, eta)+(5.05659*((-1.1752324*10^6*phi(eta)+4.1744724*10^6)*u(eta)/(int((1-eta)*(-1.1752324*10^6*phi(eta)+4.1744724*10^6)*u(eta), eta = 0 .. .5))+7.47*(diff(phi(eta), eta))/(1+7.47*phi[0])-(.1977617327*(.597+4.45959*phi(eta)))*(diff(T(eta), eta))/(1-eta)))/(.597+4.45959*phi(eta))=0;

 

diff(phi(eta), eta)-5.000000000*phi(eta)*(diff(T(eta), eta)) = 0;

 

with following boundary condition

u(0)=0,u(0.5)=0,phi(0)=1,T(0)=0,D(T)(0)=0

 

thanks for your attention in advance

I consider a PDE system (?pdsolve,numeric)

PDE := {diff(u(x,t),t)=-1/20*diff(v(x,t),x,x),
diff(v(x,t),t)=diff(u(x,t),x,x)}:

IBC := {u(x,0)=sin(Pi*x), u(0,t)=0, u(1,t)=0,
v(x,0)=1-x, v(0,t)=1, v(1,t)=0}:

pds := pdsolve(PDE, IBC, numeric, spacestep=1/50);

Suppose, I need to use the data diff(u(x,t),x,x)  or diff(v(x,t),x,x) at x=0.

So here is what I need to do within Maple.  I have to define a piecewise function and then determine left and right hand derivatives at a certain x value.  Then determine the whole derivative and graph both parts.  I am just pretty new to Maple so it is still fairly confusing and I cannot determine the exact commands that I need.

 

Hello everybody !

Do you know any mean to find the number of iterations of the simplexe algorithm used in the Minimize pocedure?

Thank you very much for your help

Martin 

In the attached worksheet there are 2 numbers. The first is the answer a got, the second is the offical correct answer.Ther look entirely different. Numerically they are the same. I wondering is there a way to simplify my answer?Rationalize_ans.mw 

Hello,

I have currently trouble running existing C-code in a new Linux (64 bit) system with Maple16. I have code that used work just fine, now it fails to evaluate. As an example, statements like

EvalMapleStatement(kv,"1+1+2^3;");

still work as expected and return a number. However, code like

EvalMapleStatement(kv,"int(x, x=0..1);");

returns int(x, x=0..1), i.e. a string, which causes my program to crash, because a number is expected.

Hello,

Please when i solve a list of equation with fsolve i get a list of result:

I am trying to solve a three-point BVP, but Maple returns the error:

Error, (in dsolve/numeric/process_input) boundary conditions specified at too many points: {0, 100, 150}, can only solve two-point boundary value problems

Is Maple just not able to solve 3-point BVP's or is there another way?

I created a Maplet that calculates ODE (Ordinary Differential Equations) however I would like to create a button that can allow printing on paper

I am using an old package, Janet, I got from http://wwwb.math.rwth-aachen.de/Janet/janet.html with Maple 16.

I changed the libname to point to the directory where the package files are.

Few days ago I was able to call help(Janet); without any problems but now everytime I try to call help(Janet)

Maple front-end loses connection to the ihe Maple kernel.

The package loads fine...

Hello,

Please, i want to resolve a non linear eaquation with Lagrange multiplier with 4 variables.

But it dosen't work.

Can you help me in this.

 

Thank you

Multiplicateur_de_La.mw

First 12 13 14 15 16 17 18 Last Page 14 of 37