Maple 18 Questions and Posts

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

Just like the title described, I have encountered an error when I use the command "expand". Actually, I just follow the example, but it doesn't work. Please help me or tell me how can I solve it in other commands.

alias(epsilon = e, omega = w, omega[0] = w0, t[1] = t1, t[2] = t2); e := proc (t1, t2) options operator, arrow; e end proc; w0 := proc (t1, t2) options operator, arrow; w0 end proc; a := proc (t1, t2) options operator, arrow; a end proc; f := proc (t1, t2) options operator, arrow; f end proc; mu := proc (t1, t2) options operator, arrow; mu end proc;
ode := (D@@2)(u)+2*mu*e*D(u)+w0^2*u+e*w0^2*u^3-e*f*cos(omega*t) = 0;
     @@(D, 2)(u) + 2 mu epsilon D(u) + omega[0]  u

                          2  3                             
        + epsilon omega[0]  u  - epsilon f cos(omega t) = 0
e_oredr := 1;
ode := simplify(subs(D = sum('e^(i-1)*D[i]', 'i' = 1 .. e_oredr+1), ode), {e^(e_oredr+1) = 0});
 / 3         2                                                 
 \u  omega[0]  + 2 (epsilon D[2] + D[1])(u) mu - cos(omega t) f

                  \                   2                   
    + 2 D[1, 2](u)/ epsilon + omega[0]  u + D[1, 1](u) = 0
simplify(collect(%, e), {e^(e_oredr+1) = 0});

u := sum('v[i]*e^i', 'i' = 0 .. e_oredr);
                      epsilon v[1] + v[0]
ode := simplify(collect(ode, e), {e^2 = 0});
for i from 0 to e_oredr do eq[i] := coeff(lhs(ode), e, i) = 0 end do;
               omega[0]  v[0] + D[1, 1](v[0]) = 0
       3         2           2                       
   v[0]  omega[0]  + omega[0]  v[1] + 2 D[1](v[0]) mu

      - cos(omega t) f + 2 D[1, 2](v[0]) + D[1, 1](v[1]) = 0
remove(has, lhs(eq[1]), cos); convert(%(t1, t2), diff);
eq[1] := %-convert(f*cos(sigma*t2+t1*w0), 'exp');

v[0] := A(t2)*cos(w0*t1+B(t2)); convert(%, 'exp'); v[0] := unapply(%, t1, t2);
       (t1, t2) -> A(t2) |- exp(I (omega[0] t1 + B(t2)))

            1                              \
          + - exp(-I (omega[0] t1 + B(t2)))|
            2                              /

Error, (in property/ConvertProperty) invalid input: PropRange uses a 2nd argument, b, which is missing
collect(%, exp(I*w0*t1));
Error, (in collect) invalid 1st argument proc (t1, t2) options operator, arrow; A(t2)*((1/2)*exp(I*(w0*t1+B(t2)))+(1/2)*exp(-I*(w0*t1+B(t2)))) end proc
coeff(%, exp(I*w0*t1));
map(proc (x) options operator, arrow; x*exp(-I*B(t2)) end proc, %);
combine(%, 'exp');
subs(I*B(t2) = I*sigma*t2-I*C(t2), B(t2) = sigma*t2-C(t2), %);
conds := combine(%, 'exp');


Hi all,

I have the following expression,

Typesetting[delayDotProduct](b, K[P], true)+k*K[D]-(Typesetting[delayDotProduct](b, m[1], true)+Typesetting[delayDotProduct](b, m[2], true)+K[D]*m[2])*k*K[P]/(Typesetting[delayDotProduct](b, K[D], true)+k*m[1]+k*m[2]+K[P]*m[2]-m[1]*m[2]*(Typesetting[delayDotProduct](b, K[P], true)+k*K[D])/(Typesetting[delayDotProduct](b, m[1], true)+Typesetting[delayDotProduct](b, m[2], true)+K[D]*m[2]))

where K and K_P are the controller gains, k is the stiffness, b is the damper, m1 and m2 are masses.

How can I find the condition on variables such that the numerator of this expression is greater that zero?

The conditions should appear as inequalities.









I have four spindle tori governed by the following equations:





I could plot the surfaces using implicitplot3d and I can imagine the volume common to these surfaces but I could not visualize it. So, am looking for a way to plot the volume covered by the surfaces such that f1<0, f2<0, f3<0 and f4<0. I know that it's easy in case of a 2D filled plot but is there any way this could be done for the 3D case? Any mathematical advice as to how to characterize or calculate this volume would also be great.



Using Maple's native syntax, we can calculate the components of acceleration. That is, the tangent and normal scalar component with its respective units of measure. Now the difficult calculations were in the past because with Maple we solved it and we concentrated on the interpretation of the results for engineering. In spanish.



Lenin Araujo Castillo

Ambassador of Maple



I want to find the i-th root of this equation s := BesselJ(1, alpha)+a1*BesselY(1, alpha)+b1*BeselJ(1, alpha)+c1*BesselK(1, alpha)


When I use the Nextzero it doesnt work



I couldn't find a command or tool for searching a value in a multi-dimensional list/array/matrix. Does it exist?. 

I wrote this:




pSearch: proc for searching

A: list, array, vector or matrix


Examples with a list:

> ss := [[1, 2, 3], [2, 3, 5], [2, 4, 5]]
> j:=1:
># search first index i such that ss[i, j]=2

> nPos := 

# search next index i such that ss[i, j]=2

> nPos := 


But I am not sure if this procedure is efficient for large lists or arrays and I would suggest (if this tool does not exist yet) a more efficient one.

César Lozada




Pretty much the title.

At first I thought something's wrong with my vectors, but even something like <1,1,1>.<2,2,2> produces this error.

I have restart in the very beginning of my worksheet and with(VectorCalculus): right after.

I'm using the 64 bit version of Maple 18 on Windows 10, if this matters.

how can save of resultes as matrix in file.dat?


How to convert maple file into the pdf file?For example the attacahed file is maple file.I want to convert this to a pdf file.

hello my friends

I have a critical problem to solve below the system of differential equations

-(9/14)*R(t)^(19/14)-(285/196)*(diff(R(t), t, t))/R(t)^(9/14)+(2565/2744)*(diff(R(t), t))^2/R(t)^(23/14) = -k*(4*lambda+1/a(t)^3)

and R(t) = 6*((diff(a(t), t))^2/a(t)^2+(diff(a(t), `$`(t, 2)))/a(t))

I need to find form of a(t) numerically from set of differential equatios

plz help me as soon as possible

This is an example, the program gave me functions _F as outputs and i want to see what F is.

This is an example.



Hi, I am using fsolve to find the roots x of a nonpolynomial expression f(x,a) = 0, using a loop for do to vary the values of the parameter a in small identical steps. At some stage, fsolve cannot find a solution. According to the Maple help, in this case fsolve returns a NULL sequence. I would like to know how to use the fact that fsolve cannot find a solution to break the loop with an if-then statement. 

Dear all,

In Maple 18, I type the following three lines and I wonder how I can then resize the resulting displayed window:


Any idea?

Hi all, since I cannot use Latex over here I try to formulate as good as possible. Suppose we have the standard integral from a to b of a function f(x), which is equal to:

h * sum (k=1 to N) [ck f(a+kh)] +O(hp), where h=(b-a)/N+1, an p(N)>= N+1;

Would it be possible to write a procedure to evaluate the coefficients ck , k=1,2,3..N and P(N)

Perhaps a stupid question, because I don't know if it is possible.

Best wishes, Math


I am a beginner with Maple. I had some basic tutorials. In fact I want to solve a first order differential equation with non-constant coefficient:

Can you really tell me what are the steps to follow in Maple 18 in order to solve this equation. I didn't find any tutorial or work that deals with variable coefficients.

Thank you in advnace.

1 2 3 4 5 6 7 Last Page 1 of 56