Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

I am running a for loop which for some loop index may give error message but it is fine for me to go with next loop indexes. but maple stop after the error message. How can I make maple to run even after the error message. Thank you

Can I do something like it?


I am trying to solve 1 equation with goniometric functions so I know the range (0..2*Pi), but it just doesnt solve it (though it solved similar ones before).

fsolve(3.*10^6 = sqrt((-4.2*10^14*(-1.*10^7*u*cos(u) + 3.84*10^8)/((-1.*10^7*u*cos(u) + 3.84*10^8)^2 + 1.*10^14*u^2*sin(u)^2)^(3/2) + 4.2*10^21*u*cos(u)/((-1.*10^7*u*cos(u) + 3.84*10^8)^2 + 1.*10^14*u^2*sin(u)^2)^(3/2))^2 + 7.056*10^43*u^2*sin(u)^2/((-1.*10^7*u*cos(u) + 3.84*10^8)^2 + 1.*10^14*u^2*sin(u)^2)^3), u, u = 0 .. 2*Pi)

Is there some way to get the solution?

Thank you in advance



I have a problem in solving a differential equation in maple. I got an error which is "unexpected occurence of the variables..." .  I have uploaded my full problem as a pdf and maple worksheet (mw) files.



What is the derivative of the composed Bessel function BesselJ(alpha, sqrt(u^2+v^2-2*uv*cos(phi))) with respect to u?

BesselJ(alpha, sqrt(u^2+v^2-2*uv*cos(phi))) is a composed function of the form f(g(u)) with f=BesselJ(u) and g=sqrt(u^2+v^2-2*uv*cos(phi)).

Best regards

I can't figure out the source of the problem.  Any ideas?




"/usr/local/opt/maple2019/lib/update.mla", `2019, May 21, 20:18 hours, version in the MapleCloud: 445, version installed in this computer: the "Physics Updates" is not installed.`

3.0*_i + 4.0*_j;


3.0*_i + 0.0*_j;

Error, (in Physics:-Vectors:-+) wrong sum of a vector with the scalar 0.


How do we handle this when it happens within a computation burried inside a proc?



How could one simplify the following expression, preferably over z in large parts of the complex plane, and with few calls under different assumptions, and preferably gracefully.

expr := 2*Pi*(-z)^(1/2)-z^(1/2)*(2*ln(-z^(1/2))-ln(z));

 Dear all

Is there a nice idea to determine the Lyapunov function of the following system

Some parameters used in the code:

b fixed parameter

a any parameter in R

Many thanks for your help



I was trying to learn more about the commands in this package and found it to be someone non satisfying:




In the attached maple worksheet, I am trying to compute an integral involving hyperbolic cosines. While not straightforward to compute, it should be completely well behaved at zero and both infinities as long as rho is real:

Much to my surprise, maple gives "undefined" as the answer, even if I attach an "assuming rho > 0" to make completely sure that the denominator never vanishes.

If I fix rho to 1 (or any other value) and compute the integral numerically, maple happily gives a numeric value. Moreover, I can scale the R, use the fact that the integrand is an even function of R, factor the denominator, and use partial fraction decomposition to transform it into the equivalent form

which maple integrates happily and gives a completely well-defined function in terms of polylogarithms, which agrees with the numeric answer.

Now, I am not expecting maple to just solve every integral I throw at it, but why does it give "undefined" as the answer to a completely well-defined integral?


I have maple code saved as a text file f.txt.

It resides here


what is the most expedient way to convert to f.mpl ?

Is there any easy way to use logical "operatives"(not sure what they are called)?


e.g., x^2[x=3] = if x = 3 then 3^2 or 0

f(x)[g(x)=x^2] = f(x) when g(x)=x^2

or whatever. I don't care about the syntax but would like it to be short and distinct. There are obviously ways to get the same behavior but they are too verbose. I like the bracket syntax but I doubt maple works with it. It could be done with a binary operator that does the comparison and returns if true...


But curious if maple has such a thing built in?

Maple gave Lie algebras of a system of PDE in which some of them do not leave the system invariant. Dont know whether the mistake is maple's or mine. File attached.


I have an arc length parametrization problem. I got the right answer for the speed. The lines of code before the long dividing line I successfully got to work. The main problem I am having is with the code underneath that. It is producing weird answers and just returning the same words without computing any mathematical calculation.




T:= 4:

r := t -> <t^2 + t, sin(t^2)*(t + 1), cos(t^2)*(t + 1)>;

speed := Norm(diff(r(t), t));

evalf(Int(speed, t = 0 .. T)); (I got 62.98633182 for this part)

----------------------------This is where I started running into problems with the arc length parametrization.

L := b -> int(speed, t = 0 .. b);

speed := t -> subs(c = t, Norm(diff(r(c), c)));

speed2 := t -> sqrt(factor(simplify(speed(t)^2)));

solve(s = L(t), s);

assume(b > 0 'real');

g := s -> solve(s = L(b), b, useassumptions = true);

newr := s -> r(g(s));



First 9 10 11 12 13 14 15 Last Page 11 of 1604