Hello,I just can't find out whether it's possible to plot the x-axis 'below' the graph? Normally speaking the x-axis is plotted at y=0, but I want my x-axis to be plotted for example at y=-5 or y=-10. Can that be done?Greetings,Frank

Hello every one,

While plotting contours with grid size [200,200] my computer consume too much time to

execute the commond. To have an idea, the current process which I am running right now already

took 2380.22s and 68.60M and still its on.

Any suggestion to make the contour plot process faster?

Thanks

Good day,

I have plotted a family of lines in one graph, and now I want to shade an area between two horizontal (parallel) lines, for example 2 and 5.
So what I did was to use the command 'filled'. First I shaded the area in vertical direction from 0 to 2 white. After that I filled the area in vertical direction from 0 to 5 for example blue.
This results in a shaded area between 2 and 5 in the color blue, what I like to have.

Question: Is there a way (a command) to fill an area between two horizontal lines more easily? So that I shade it between 2 and 5 in just one commandline. Now I first make a part white, and then the next 'layer' behind is in the color I'd like to have.

Another thing is that when I shade an area in the color white, my gridlines disappear! Is there a way to avoid this?

I have found information about shading area between two intersecting lines (can be done relatively easy), but between two horizontal (parallel) lines I still can't find a simpler way to do it.

Somebody got a suggestion?

Greetings,
Frank

 

I recently purchased a new MacBook Pro (retina) with Mountain Lion and found that my Maple 12 programs would not run.  So I updated to Maple 16.  Initially all seemed to work well.  But now when i try to open a Maple .mw file, I get the message

"File acm23test.mw could not be found.  Please select a different file."

What has happened and how can I get my programs working again?

DavidR

Hi,

I have the following code which solves a system od ode (3 ode) with constraint parameter. actually, this algorithm is time consuming and is not effective as I see that some people solve this odes with matlab in an efficient (qucik) way. So, I ask you if it is possible help me to improve the efficiency of the following code

thanks

restart:
eq1 := diff(u(eta), eta, eta)+diff(u(eta), eta)/ (eta+zet)+Gr*Pr*T(eta)-Nr*phi(eta)+sigma1;
eq2 := diff(T(eta), eta, eta)+diff(T(eta), eta)/(eta+zet)+Nb*diff(T(eta), eta)*diff(phi(eta), eta)+Nt*diff(T(eta), eta)^2-u(eta);
eq3 := diff(phi(eta), eta,eta)+diff(phi(eta), eta,eta)/(eta+zet)+ Nt/Nb*(diff(T(eta), eta, eta)+diff(T(eta), eta)/(eta+zet))-Le*u(eta);
zet:=0.1:
Nt:=0.05;
Nb:=0.05;
Le:=10;
Pr:=1;
Gr:=5;
Nr:=0.1;
Ree:=5;
#sigma:=2;



p:=proc(pp2) option remember; global res,F0,F1,F2;
if not type([pp2],list(numeric)) then return 'procname(_passed)' end if:
res := dsolve({subs(sigma1=pp2,eq1)=0,eq2=0,eq3=0,u(0)=0,u(1-zet)=0,phi(0)=0,phi(1-zet)=0,T(0)=0,T(1-zet)=0}, numeric,output=listprocedure);
F0,F1,F2:=op(subs(res,[u(eta),phi(eta),T(eta)]));
evalf((2/(1-zet^2))*Int(F0(eta)*(eta+zet),eta=0..1-zet))-Ree*Pr;
end proc;

fsolve(p(pp2)=0,pp2=(-80)..(80));   # I dont know what exactly this parameter has
#res := dsolve({eq1=0,eq2=0,eq3=0,u(0)=0,u(1-zet)=0,phi(0)=0,phi(1-zet)=0,T(0)=0,T(1-zet)=0}, numeric);
 

sigma1:=%;
                             
res := dsolve({eq1=0,eq2=0,eq3=0,u(0)=0,u(1-zet)=0,phi(0)=0,phi(1-zet)=0,T(0)=0,T(1-zet)=0}, numeric);
with(plots):
odeplot(res,[[eta,u(eta)],[eta,phi(eta)],[eta,T(eta)]],0..1-zet);

 

Hi all,

This might be more a mathematical problem than a Maple problem, but still maybe somebody could help me out.

I am trying to solve the differential equation of an Euler-Bernoulli beam with an distributed load (q1) and a varying axial force (N(x)). The following equation applies:

 

 

Hello guys,
I'm newb with Maple and trying to solve the following nonlinear recurrence system:

rsolve ({L(0)=1,L(n)=1.025*L(n-1),h(n)=L(n)/(1.1*L(n)-0.1*v(n-1)),v(0)=0.5,v(n)=(1+h(n))*v(n-1)},{L,h,v});

But Maple returns me nothing, no error. What's might be the problem?

Thanks in advance.

Hello,

I've just asked maple to plot a(s) from s=0 ..10 where a(s) is the absaolute value of a complex piecewise function. When I enter the command to plot:

> plot(a(s), s = 0 .. 10, discont = true);

maple just shows 'Evaluating' and fails to produce anything. I've left it for up to 20 minutes without any result. Maple has successfully plotted very similar things before in a matter of seconds so I don't think time is an issue here. The full script is as follows:

Hi everyone,

I am new to Maple, and require the use of grtensorii in my research. The package is installed correctly and will open fine, but it has troubles in finding the metric files to use. When I try to qload the metric I want I get the error message:

Error, (in qload) Could not find metric perturbed.

I have spoken to some colleagues who ran into similar issues making grtensor compatible with Maple 15, and they ran the commands

unprotect(sqrt);

Hello, Is there a way to assign mutiple variables at once using the 'if-then' statement?For example: if A:=3, then B:=434 and C:=52 and H:=2039So that I fill in these variables B, C and H one time and then when A is chosen to be a value, all other values are assigned automatically. Like filling in a table actually.In a simple 'if-then' statement I assign one variable, which works well.a:=10if a:=10 then b:=4elif a:=20 then b:=3end if:b ends up being 4I don't know how to add variable C and H, without getting an error.It should be very simple in my opinion, but I just don't see it.Greetings,Frank

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.