Maple 13 Questions and Posts

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

with(Physics) :
t:=Intc(Dirac(k1+k2+k3)*phi(k1)*phi(k2)*phi(k3),k1,k2,k3) ;

# how to force

# to find
3*Intc(Dirac(k2+k3-k)*phi(k2)*phi(k3),k2,k3) ;


I am thinking about plotting a star chart of the whole sky. While it might not be difficult to get the (x,y)-coordinates of fairly good number of stars from star catalogues and to plot them accordingly in a 2d-plot, it might be more difficult to plot the stars according to their brightness. Also, I would like to plot the chart in cylindrical coordinate with astronomical coordinate axes for x and y. Moreover, only the brighter stars should be visible when the whole chart is shown, but when zooming in (while still having the axes visible) to have a closer look at some specific locations where I wand to make some square to show some fields, fainter stars should appear.

Is this possible in Maple ? I'd prefer to do a plot as I can really draw several squares where I want them which is not possible in my planetarium software.

   Hi, there

 I work with maple 13 and made a file that its contet is plotting curves and its size is 

151 MB.Now when I want to open it the file cannot be opened completely only about a quarter of its contents rises.


Thanks for your help



Hello everyone!

I'm pretty new with Maple. I think I've understood the way Maple handles differentiation fairly well, but upon a specific request from my PhD tutor I have to perform a task which is giving me a hard time. 

My question is: is in any way possible to express a derivative of a function or expression in terms of the function itself?
I'll try to explain myself with an example: let f(x) and f'(x) be the function and its first derivative:

Instead of expressing f'(x) in the way shown, I'd like to express it as a function of f(x), such as in the following:

I would apply the same process to the higher order derivatives, if possible.

A huge thank you to whoever will help me!

Dear all

I want to know, how one can install third party package into Maple13, the package is "wkptest" i downloaded it from link If anyone knows how to do this please help.

  Hi there,

  I want to use maple 13 for calculating mean value theorem for differentiable function  f:=piecewise(-2≤x≤0, -x2 ,0≤x≤2,x2

on the interval [-2,2]. But an error occured, that is,"function must be continuous".Any help will be appreciated.




Has there been any progress with new commands for Maple 13 that allow one to auto-resize the plotting grid?

I am generating fractals with Maple 13 and resizing the plot output grid manually is not an option because it distorts it. Neither is the option to resize it once and recalculate, because the final grid contains many points.

Any pre-processing or massaging code prior to executing the main code is very welcome (if it can be done).

Many thanks,


Matlab seems to be pretty strong at doing color plots with separate color bars, e.g.

Is this also possible in Maple and somehow in combination with `plots[surfdata](...,color=zhue,...)`?

In a 3d coordinate system I have a circular spacecurve with z-minimum -4 and z-maximum +4. In the same 3d coordinate system I have a 3d surface plot with z-minimum -0.5 and z-maximum +1.3 . When I choose the color option "Z(Hue)" in order to color-code the z-values on the 3d surface and make the topography more clear, I mostly get a totally green 3d surface. It seems that the color scaling is coupled with the spacecurve with z-values of +-4 . How can I uncouple the color scaling from the spacecurve and couple it with the z-range of the surface, while the color-limits shall be at +-1.3 ?

Hello i want to solve the differentiel equation but i have these problem i don't understund  why !?


ode1 := (K[Q]*T*R[a]/K[kol]+R[a]*B[m]*sqrt(T/K[kol]))/K[i]+K[b]*sqrt(T/K[kol]) = 0

(K[Q]*T*R[a]/K[kol]+R[a]*B[m]*(T/K[kol])^(1/2))/K[i]+K[b]*(T/K[kol])^(1/2) = 0


ode2 := (1/2)*(-(4*(diff(theta(t), t)+theta(t)))*M+2*B1*(diff(theta(t), t)+theta(t))-2*w[2]*sin(theta(t))+m1*g*sin(theta(t))*l[kol]-2*B1*theta(t)+2*w[1]*sin(theta(t))-2*m1*g*sin(theta(t))*l[2])/l[1]

(1/2)*(-4*(diff(theta(t), t)+theta(t))*M+2*B1*(diff(theta(t), t)+theta(t))-2*w[2]*sin(theta(t))+m1*g*sin(theta(t))*l[kol]-2*B1*theta(t)+2*w[1]*sin(theta(t))-2*m1*g*sin(theta(t))*l[2])/l[1]



ode3 := subs(T = (1/2)*(-(4*(diff(theta(t), t)+theta(t)))*M+2*B1*(diff(theta(t), t)+theta(t))-2*w[2]*sin(theta(t))+m1*g*sin(theta(t))*l[kol]-2*B1*theta(t)+2*w[1]*sin(theta(t))-2*m1*g*sin(theta(t))*l[2])/l[1], ode1)

((1/2)*K[Q]*(-4*(diff(theta(t), t)+theta(t))*M+2*B1*(diff(theta(t), t)+theta(t))-2*w[2]*sin(theta(t))+m1*g*sin(theta(t))*l[kol]-2*B1*theta(t)+2*w[1]*sin(theta(t))-2*m1*g*sin(theta(t))*l[2])*R[a]/(l[1]*K[kol])+(1/2)*R[a]*B[m]*2^(1/2)*((-4*(diff(theta(t), t)+theta(t))*M+2*B1*(diff(theta(t), t)+theta(t))-2*w[2]*sin(theta(t))+m1*g*sin(theta(t))*l[kol]-2*B1*theta(t)+2*w[1]*sin(theta(t))-2*m1*g*sin(theta(t))*l[2])/(l[1]*K[kol]))^(1/2))/K[i]+(1/2)*K[b]*2^(1/2)*((-4*(diff(theta(t), t)+theta(t))*M+2*B1*(diff(theta(t), t)+theta(t))-2*w[2]*sin(theta(t))+m1*g*sin(theta(t))*l[kol]-2*B1*theta(t)+2*w[1]*sin(theta(t))-2*m1*g*sin(theta(t))*l[2])/(l[1]*K[kol]))^(1/2) = 0



(1/2)*(-4*R[a]*K[Q]*M*(diff(theta(t), t))-4*R[a]*K[Q]*M*theta(t)+2*R[a]*K[Q]*B1*(diff(theta(t), t))-2*R[a]*K[Q]*w[2]*sin(theta(t))+R[a]*K[Q]*m1*g*sin(theta(t))*l[kol]+2*R[a]*K[Q]*w[1]*sin(theta(t))-2*R[a]*K[Q]*m1*g*sin(theta(t))*l[2]+R[a]*B[m]*2^(1/2)*((-4*M*(diff(theta(t), t))-4*M*theta(t)+2*B1*(diff(theta(t), t))-2*w[2]*sin(theta(t))+m1*g*sin(theta(t))*l[kol]+2*w[1]*sin(theta(t))-2*m1*g*sin(theta(t))*l[2])/(l[1]*K[kol]))^(1/2)*l[1]*K[kol]+K[b]*2^(1/2)*((-4*M*(diff(theta(t), t))-4*M*theta(t)+2*B1*(diff(theta(t), t))-2*w[2]*sin(theta(t))+m1*g*sin(theta(t))*l[kol]+2*w[1]*sin(theta(t))-2*m1*g*sin(theta(t))*l[2])/(l[1]*K[kol]))^(1/2)*l[1]*K[kol]*K[i])/(l[1]*K[kol]*K[i]) = 0


eol := (1/2)*(-4*R[a]*K[Q]*M*(diff(theta(t), t))-4*R[a]*K[Q]*M*theta(t)+2*R[a]*K[Q]*B1*(diff(theta(t), t))-2*R[a]*K[Q]*w[2]*sin(theta(t))+R[a]*K[Q]*m1*g*sin(theta(t))*l[kol]+2*R[a]*K[Q]*w[1]*sin(theta(t))-2*R[a]*K[Q]*m1*g*sin(theta(t))*l[2]+R[a]*B[m]*sqrt(2)*sqrt((-4*M*(diff(theta(t), t))-4*M*theta(t)+2*B1*(diff(theta(t), t))-2*w[2]*sin(theta(t))+m1*g*sin(theta(t))*l[kol]+2*w[1]*sin(theta(t))-2*m1*g*sin(theta(t))*l[2])/(l[1]*K[kol]))*l[1]*K[kol]+K[b]*sqrt(2)*sqrt((-4*M*(diff(theta(t), t))-4*M*theta(t)+2*B1*(diff(theta(t), t))-2*w[2]*sin(theta(t))+m1*g*sin(theta(t))*l[kol]+2*w[1]*sin(theta(t))-2*m1*g*sin(theta(t))*l[2])/(l[1]*K[kol]))*l[1]*K[kol]*K[i])/(l[1]*K[kol]*K[i])



R[a] := 1.42;



K[b] := 0.649e-2;



K[i] := 0.649e-2



K[Q] := 0.1051618298e-6



l[kol] := 1



l[1] := .7



l[2] := .3



K[kol] := 0.1168464776e-5



B1 := 0.955e-3



B[m] := 0.955e-3



J := 0.475e-6



M := 0.91e-2



m1 := 0.726e-1



w[1] := 0.72e-1



w[2] := .45



g := 9.81



a1 := 0




-.4851223862*(diff(theta(t), t))-.5119876735*theta(t)-6.626549550*sin(theta(t))+.1077211171*2^(1/2)*(-42167.66273*(diff(theta(t), t))-44502.83918*theta(t)-575990.9557*sin(theta(t)))^(1/2)



-.4851223862*(diff(theta(t), t))-.5119876735*theta(t)-6.626549550*sin(theta(t))+0.1523406647e-3*(-0.4216766273e11*(diff(theta(t), t))-0.4450283918e11*theta(t)-0.5759909557e12*sin(theta(t)))^(1/2)



ic1 := theta(0) = a1

theta(0) = 0



dsol1 := dsolve({eol, ic1}, numeric, output = listprocedure, range = 0 .. 10)

Error, (in DEtools/convertsys) unable to convert to an explicit first-order system


plots[odeplot](dsol1, [t, theta(t)], 0 .. 10)

Error, (in plots/odeplot) input is not a valid dsolve/numeric solution



thanks for your help


It seems a frequent issue that exported 3d plots are not shown as wished. I experience the same problem. Although I exported in the .eps format into a .tex latex-file the resulting .pdf-file shows a somewhat pixelated image of my 3d plot as if it was created in "Paint". Is there a solution for this in Maple13?

Take a 3d plot of some uneven surface in the xyz space and you want to have the surface colored according to local z-coordinates (e.g. a valley is blue, a peak is red with all rainbow colors inbetween). For such a color-coding one can in principal select the plot option "Color->Z (Hue)". What do you do if the valley and the peak are still more or less green colored? How can you force the valley and peak to have different colors?

I have this surface which I display with the following commands:



It all works, and the surface has nice default rainbow colors, but when I choose Color->Z (Hue) everything turns to green. Also when I manually write in the first of the above lines "color=["Blue","Red"]" the surface turns all red. The ranges of the x and y axis are larger than the z axis by a factor of about 50. Might this be the problem? How can I adjust the color ranges for my small z axis?

How would you insert a label-coordinate in this command:




What is the command for filling the space between two circular spacecurves? Something like a Saturnian ring?

   Hi, there

Hoe can we plot the floor function by Maple13? I used plot command and implicitplot command. But the output was wrong on integer values.

Thanks for your help



First 16 17 18 19 20 21 22 Last Page 18 of 52