Maple 15 Questions and Posts Maple 15 Questions and Posts Feed

These are Posts and Questions associated with the product, Maple 15


I have a problem with the last dsolve in the following code. I am asking it again since nobody answerd me.

The last dsolve in the code " dsolve({ss[ii], icss[ii]}); " connot integrate theta[1] twice and some "Z1" apears on the screen !BUT when I integrate theta [1] seprately, maple gives the best answer but it does not work in this code !
please help me

B := 1+(2*(1-zzeta))*Br*(int(tau*(diff(u(tau), tau))^2, tau = (1/2)*zzeta*(1-zzeta) .. 1/2*(1-zzeta)));
eq4 := 4*B*u(tau)-(1+zzeta)*(diff(tau*(diff(theta(tau), tau)), tau))/tau-(1+zzeta)*Br*(diff(u(tau), tau))^2;

HH:= p*((4*(1+(2*(1-zzeta))*Br*(int(tau*(diff(u(tau), tau))^2, tau = (1/2)*zzeta*(1-zzeta) .. 1/2-(1/2)*zzeta))))*u(tau)-(1+zzeta)*(diff(theta(tau), tau)+tau*(diff(theta(tau), tau, tau)))/tau-(1+zzeta)*Br*(diff(u(tau), tau))^2)+(1-p)*(diff(theta(tau),tau$2)):

convert(series(collect(expand(eq5), p), p, nn+1), 'polynom');

for ii to nn do
ss[ii] := (coeff(eq7, p^ii)) ;
print (ii);
end do;

ss[0]:=diff(theta[0](tau), tau, tau);

icss[0]:=theta[0](zzeta/(2*(1-zzeta)))=0, D(theta[0])(1/(2*(1-zzeta)))=1;

dsolve({ss[0], icss[0]});
theta[0](tau):= rhs(%);

for ii to nn do
icss[ii]:=theta[ii](zzeta/(2*(1-zzeta)))=0, D(theta[ii])(1/(2*(1-zzeta)))=0;
dsolve({ss[ii], icss[ii]});
end do;

I would be most grateful if you help me to find this problem.

Thanks for your attention in advance

I'm trying to display some flow lines for a vector field. The vector field I have is:

What I put into Maple is:

I've tried this with a number of points, but I always get an error message along the lines of:

"Error, (in Student:-VectorCalculus:-FlowLine) cannot determine if this expression is true or false: (2.*(1)(0)-4.)^2+(-2.+4.*(1)(0))^2 < 0.9e-9"

Could someone please shed some light as to what's going wrong for me? I would be very grateful for any help.


Hi there!

I wrote a piece of code which spits out the numerical datapoints (x,y(x)) corresponding to a function y(x). So that the result is accurate, I need quite a lot of data points - currently I am working with 5k.

In order to work with this function later, I interpolated it with a Spline. For instance, I would like to sample the function values on a fifferent grid, etc.. However the evaluation of this function really takes up hell of a lot of time, and the reason seems to be, that it, being a spline on 5k nodes, is simply a huge expression.

Is there a better way to do this? Are other fitting functions than a spline maybe better suited?

Thanks for help!


I need to get rid of the type "constant" for gamma. In Maple, gamma is defined as Euler's constant by default. While it is easy enough to unprotect(gamma) and then get rid of its value, Maple will refuse to solve an equation for gamma, as it remains of type constant even after deassigning it. So I need to regain gamma as a variable.

Some may feel this is an unwise thing to do. But it actually is not: I am writing a document involving physics, and gamma is the accepted symbol for the relativistic energy. I cannot avoid using that, lest mass confusion ensues (this involves students). I really don't want to write gammar instead. Euler's constant, otoh, does not figure at all in my document.

Note that I need a solution that works in Maple 15 and later as I am working in a heterogeneous environment as far as Maple versions are concerned.

Thanks in advance,

Mac Dude


i have got alot of mixed and high degree derivatives. For example:

u[x]*u[x,t]*eta[x,t]+u[]^2*u[x]*eta[x]+kis(x,y)u[x,t]^2*u[]+eta(x, y)*u[]*u[x]^2+ksi[x,t]*u[x]^2*u[x,t]+......

like this alot of terms

my question is how can i solve divided by the derivative of the u(x,t) partial differential equations system and so  how can i find eta(x,t,u) and ksi(x,t,u) 

I have the following expression (obtained from an earlier calculation):

I want to collect all the terms under one summation. So I define a rule:

end proc:

and then


I get

Error, (in +) unable to identify A::algebraic

I used similar constructs before so I think the rule is constructed correctly. I should, however, mention that I use the Physics:-Vectors package and in fact the expression I start up with here reads, in 1-d Maple inputform:

Physics[Vectors][`+`](Physics[Vectors][`+`](Physics[Vectors][`+`](-y*(Sum((diff(a[n](r), r))/(exp(I*Pi*n/L))^2, n))/r, (2*I)*(Sum(a[n](r)/(exp(I*Pi*n/L))^2, n))*k0), y*(Sum(a[n](r)/(exp(I*Pi*n/L))^2, n))*k0^2), -y*(Sum((diff(a[n](r), r, r))/(exp(I*Pi*n/L))^2, n)))

Is my problem related to the use of Physics:-Vectors? If so, how can I get around that?


Mac Dude

Hi everyone

I'd be pleased if you could give a hand with the exploration assistant.


1. I want the exploration assistant to appear on the same document I am working on, but everytime I use it (either by right-clicking or by the explore command) it automatically appears on a new document.


2. can I manipulate a piecewise function when using embedded components? i.e.: plot the function and varying the parameters using sliders.




a error and code in this attachment mw

i can pdsolve it, but numeric pdsolve it get error

I have this kind of problem. When I try to solve my command:


is there any way to avoid this? in order to get the answer directly...

my question may be simple but I don't find the answer in any help guide.
when I define a function I cannot use a linearalgebra expression such as Trace.
Here is an example of what I would like to do:

If anyone can help me...
Thank you


Say, I got an expression that depends on two variables, x and y. How can I tell Maple, that y is actually just a (real) constant, so y does not depend on x?

Because when I apply a differentiation with the "D" - command, it would always also write out expressions, where y is differentiated w.r.t. x.



this question has been asked several times here already accoruding to the search function, but I didn't find an actual answer to it.

My Maple 15 gives outlarge numbers without scientific notation, so even 10^20 would be displayed as 1000....000  (20 zeroes). It seems like this is not the standard, so at one point I might have changed the way Maple displays the output, and cannot remember anymore how.

How can I tell Maple to display the output in scientiffic notation, and how can I set up the threshold of powers of ten, when this is done. (so that e.g. 100 still gets disaplyed as 100, and not as 1*10^2, but anything above gets displayed in scientific notation).

I know that I can right click on the output to change it's display form. But I want to change this generally, so I dont have to do this for any single numerical output I get.

I use Maple 15, running in Mac OS.

Thanks for help!

I have to plot 4 vectors in one graph. I have the following structure:

plot(Vector([12, 12.5, 13, 13.5, 14, 14.5, 15]), Vector([1.622712644, 1.265443137, 1.028604736, .8605013333, .7352916667, .6386248233, .5618945274]), style = line, symbol = asterisk, color = blue)


plot(Vector([12, 12.5, 13, 13.5, 14, 14.5, 15]), Vector([5.483608580, 4.289400489, 3.496793877, 2.933480578, 2.513320599, 2.188469637, 1.930230220]), style = line, symbol = asterisk, color = blue)

Could you help me to plot these two curves in one graph.


please help me to find a solution for this issue...

I would like to thank you in advance

Best regards,


I tried to get the maximum and minimum values of the following function. From the plot I get them but its not accurate. Please advise me to get them accurate.





implicitplot(((N^2+B^2-F*N)*K=N*(N^2+B^2-F*N+F*N),K=0..10,N=0..10,view=[0..5,0..4],numpoints=90000,axes=boxed,thickness=2,color=black,font=[1,1,20],tickmarks=[3, 3],linestyle=1));


5 6 7 8 9 10 11 Last Page 7 of 43