Maple 18 Questions and Posts

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



I was wondering how to go about plotting a Fourier Tranform in Maple.

My assignment is to plot a simple harmonic equation as a Fourier transform, depicting amplitude against fequency.

I've been given: x'' + w^2 x = 0

And want to obtain both the f(x) = a0 sin(wt) + b0cos(wt) form, and a graph of the the amplitude (c^2 =a0^2 + b0^2) against frequency.

I know how to do this on paper but not in Maple, so any help with line commands and layout would be very much appreciated.





I have an issue calculating an electronics circuit with Maple, using units. I have a current source that I know, and I want to determin the voltage in a capacitor by solving an ODE (except that the current source is defined piecewise). And to make sure I have all the units and scales right, I use the standard unit package. All my variables have their units defined.

Except that Maple doesn't want to solve the equation. It seems to me that it assumes that the function I am trying to solve is unitless, and therefore refuses to solve. 

V__out := 3*Unit('kV');

C__out := 2*Unit('nF');
R__blead := 520*Unit('`kΩ`');

I__fly := proc (t) options operator, arrow; Unit('A')*piecewise(t < 3.25*Unit('us'), (1+(-1)*t/(3.25*Unit('us')))*.2, 0) end proc;


dsolve({I__fly(t*Unit('s'))-V__C(t*Unit('s'))/R__blead = C__out*(diff(V__C(t*Unit('s')), t)), V__C(0*Unit('s')) = V__out}, V__C(t*Unit('s')));
Error, (in Units:-Standard:-+) the units `A` and `S` have incompatible dimensions


Is there a way to make Maple assume the unit of what it's trying to solve ? I need it to understands that V__C is in Unit('V') ...




It took me a while to figure out what was causing this unusual sort result;

I loaded the MTM package because I want to convert numbers to integers for indexing. This sort result seem like it could be a bug.

Sort a List of Floats




x := [1/3, 0, 4] = [1/3, 0, 4]``


sort(x) = [0, 1/3, 4]``


f_x_MTM := [seq(convert([x[n], n], float), n = 1 .. 3)] = [[.3333333333, 1.], [0., 2.], [4., 3.]]NULL

sort(f_x_MTM) = [[.3333333333, 1.], [0., 2.], [4., 3.]]NULL



f_x := [seq(convert([x[n], n], float), n = 1 .. 3)] = [[.3333333333, 1.], [0., 2.], [4., 3.]]


sort(f_x) = [[0., 2.], [.3333333333, 1.], [4., 3.]]




Variable exists but is not shown when using save/read with extension .m

My example maple code is



save a, "test1.m"


read "test1.m"



after the read command, I can access the variable, but it is not shown under "Variables".

This lead to some confusion when debugging the worksheet. Can I change this somehow?

Using the input type file format is not a solution, since then reading takes forever for complicated expressions.

Further, in the read command documentation it says "This functionality is not intended for end users" for saving the file as .m. What does that mean?

Hi everyone,


Consider this maple 18 doc:


The code is regular code for Julia sets of the exponential.


To see how the Julia set behaves at infinity, I apply the transform mu(z)=1/z.


The plot3d command correctly plots the Julia set at an appropriate neighborhood of infinity, but:

1) Axes are not transformed

2) Saving as .eps produces an empty plot, sans the axes (plot is saved correctly, when not applying mu(z))


Is there any trick to force the axes to also show correctly with the transformed ranges?

Seems that this misalignment is bothering the .eps renderer, which probably plots the graph in twilight zone, given the false ranges of the untransformed axes.


Any ideas on how to force the saveas .eps to work in this case?


Many thanks,


I would like to plot the following singular double integral, but I cannot due to singularities...


where x>0, t=0.2 and m=0.2.

I defined f(y) function as f:=y->exp(-(y-4.68)^2/0.4):

I attached my file:

Thank you !

i have this problem -> f'^2 -ff''=f'''-k1(2f'f'''-ff''''-f''^2)+Ha^2(E1-f') with boundary conditions f(0)=0, f'(0)=1, f'(∞)=0.

since it is a fourth order equation, but only three bcs, it does not produce unique solution. so the solution of the equation may be seek in form of f=f0(eta)+k1f1(eta).

thus the equation will become 




boundary conditions are 



i had been clueless in solving this problem. please somebody help me with this problem.

Given the sequence defined by the recursive relation a[n+1] = r*a[n](1-a[n])
You need to use the procedure iterate.
Throughout this problem you should choose initial values in the interval 0<a0<1.
(a) Let r=3/2. Calculate a moderate number of terms in the sequence (between 10 and 20). Does the sequence appear to be converging? If so to what value? Does the limit depend upon your choice of initial value? Plot the terms you have calculated
(b) Let r=2.8. Calculate a moderate number of terms in the sequence (between 10 and 20). Does the sequence appear to be converging? If so to what value? Does the limit
depend upon your choice of initial value? Plot the terms you have calculated How does this sequence differ from that in part (a).
(c) Let r=3.2. Calculate a moderate number of terms in the sequence (between 10 and 20). Show that the sequence does not appear to converging. Plot the terms you have calculated and describe how the sequence behaves in this case.
(d) Consider intermediate values between 2.8 and 3.2 to determine more precisely where the transition in behaviour takes place. Provide a few plots (no more than 4) showing the values you have investigated.
(e) Consider the values of r in the range 3.43<r<3.46. Determine as accurately as you can the value of r for which the period of oscillation doubles.
(f) As r increase further period doubling occurs. Try to find the when the sequence appears to oscillate between 8 values.
(g) Let r =3.65 and calculate a considerable number of terms (at least a few hundred) and plot your values.
(h) For r=3.65 choose a0=0.3 and then a0=0.301. Find and plot some terms in the sequence for each initial value. Determine how long the terms in the two sequences remain close together and when they begin to depart significantly from each other.

I have used data and the curve fitting procedure to find a curve and I now have to put error bars on the graph to complete my task? I know that there is an errorplot function but I cant seem to find a way to do it with the curvefitting thing?


I have attached the whole worksheet, I hope someone can help! Thanks


Even if there was some method of exporting the graph to another medium that would also be very helpful




I am pretty new to Maple and Im trying to find the parameters of an equation using some kind of fit routine but I can only find such a routine to fit an expression to data and not vice versa.

my equation is as follows


epsilon(E):= a0 + a1*ln(E) + a2*ln(E)^2 + a3*ln(E)^3


I do have data for E which I imagine I need. 


E:= {121, 244, 344, 411, 444, 778, 867, 964, 1085, 1112, 1212, 1299, 1408} all in keV :)


Any suggestions/help would be much appreciated, although I am new to the program so go easy on me :) 




Hi Mr Preben Alsholm .

i before asked you about a i have another question...

for determine unknown parameter oemega(or eigenvalue) according to line ''newsys2 := subs(omega^2 = omega2*10^5,newsys)'' and also line ''approxsoln = [omega2 = 1, f(x) = x^2*(1-x)^2]''

power (5) in term'' 10^5''  in  (omega2*10^5,newsys) and initial guess (1)  (omega2 = 1) are very effective on final result eigenvalues .sometimes gain different eigenvalues ,which it is not impossible recognize that

which one is correct and convege.however i need Positive

 eigenvalue that is minimum between them,If and only if , it converge after
some iteration in two section of maple file which is attached as

  According to my code converge occuerd until  eigenvalue gained in first section (before line `11` := ((1+6*alpha2)*(1/12))*(int(fy11^2, x = 0 .. 1))...........) are equal with those obtain in second section .(please see below file for example

it is necessary mention that between section 1 and 2 is relations.for example amont of first obtained eigenvalue which obtain from first section must be repleaced in line ''approxsoln = [omega2 = 0.661514014001420, h(theta) = theta^2*(1-theta)^2]...........''in second section,

and this procedure should be continued between section 1 and 2 until convergence with desired accuracy occurred.Another relate is that,the first ODE system  can be solved using the first set of boundary data to obtain first estimate for  second section (fy11). At the next step by repeating in the same manner, 

this time by obtained function (g) at the end of maple code (Solution of the second ODE system  with second set of boundary conditions leads to the first estimate of function g),that at this stage, the first iteration is completed.

Now by replacing omega2 which is determined in this section in to the first section ,fy11 is updated and gained Further.

Next, the updated function , by continuing the iterative procedure.

in matlab bvp4c rule is used for this this impossible in maple software ?if not please help me for solve and gain correct omega..

thanks alot

I am doing a Calculus assignment and I can't find the commands for certain things.

1.Given the function f(x) = ((x+1)^2) / (1+x^2)

i) The domain of continuity of f(x)

ii) The intervals of increase and decrease of f(x) by using test points.


2. Use the IVT to prove existence of a root to the equation x^3 +10x^2 -100x +50=0 in the interval [-20,10]. Use again the IVT to show that there is a 1st root in [-17,-15], a 2nd toot in [0,1] and a 3rd root in [ 5,6]. Find or approximate those roots with Maple. (the bolded is what I need help).



I am working on a Windows machine. In older Maple versions, decreasing the font size for the screen display (via zoom factor) also decreased the font size for the printout. But now the printout comes with a pretty large font size and does not change when using the zoom factor. It practically makes no sense anymore to print. How can I decrease the fonts for the printout to my liking?



Hello altogether,

I want to plot the numerical result of an ODE, which seems to be pretty simple at first sight, but the difficulty is that the boundaries are depending on the solution.

The following pseudo-code describes what I want to have, but it doesn't work. This code fills the RAM pretty fast and you will have to kill the process.

Is it possible to calculate a solution to this problem numerically (or even analytically) and if yes, how?

Since I am new here, I am sorry for any bad-to-read maple code or any noob errors I have made. I would be very thankful, for any response and help.



hi.i have problem for dsolve equation.please help me again


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