MaplePrimes Questions

I am trying to solve a transcendetal equation that has no analytical solution (AFAIK).
My question is how to find different branches of solutions with Maple. So far, Maple
only returns one particular solution.

To fix ideas, my equation is of the type,
x-a*(1-x)* exp(bx)+c=0
For suitable values of a and b this equation seems to have, at least, two solutions in
the real axis and I would expect it to have a numerable infinite of solutions in the



I have an expression, say 


v(t) = .5*(sum(abs(q)*(J*exp(q*j*Wc*t)+K1*exp(q*j*phiK1)*exp(q*j*(Wc+Wx)*t)+K2*exp(q*j*phiK2)*exp(q*j*(Wc-Wx)*t)), q = -1 .. 1))


When I raise this quantity to a power, say (v(t))^3, I want the answer as well in the form of a summation. This system is frequently encountered for a voltage source, and I am looking to segregate the harmonics generation through the use of the...

System: Mac OS X 10.8.1 + Maple 16

I am trying to copy + paste a selected region from an image shown in Preview to Maple. The image is blank when it is pasted. Why is this not working as expected?

Chemistry / Chemical & Materials Engineering

He := sum(((-1)^k)*n!/(k!*(n-2*k)!*2^k)*x^(n-2*k), k=0..n/2);
Gen := simplify(sum(He*z^n/n!, n=0..infinity));

convert(Gen, ratpoly);


I have solutions of a 4th-order differential equation as a solutions.dat file having 5 columns. In the same maple worksheet (not necessarily the same worksheet, though), I want to solve a system of differential equations using [columns 2] and [column 3]/[column 2] of soultions.dat as inputs.

Can anyone help me how I can implement this?

Thanks in advance.



I'm trying to implement tan(x) over [-Pi/4,Pi/4] to have < 4ulp error.

How can I generate rational polynomial approximation so that that can be used to implement tan(x) in C.

Please suggest any references.


Thank you.

Hi every one,

I am trying to find a way to plot a dependent variable against a parameter. 

The data I want to use for this ploting should be from the numerical results

obtained by using shooting method.


restart:ShootLib := "C:/MP/Shoot9/":
> libname := ShootLib, libname:
> with( Shoot ):
> with(plots):with(LinearAlgebra):
> A:=0:
> FNS:={ f(z), G(z), H(z)}:
> ODE:={diff(f(z),z) = G(z),

Hey maple gurus

when i solve and equation with two real roots (1 positive and 1 negative roots), how do i automatically pick out the postive solution?

As you can see from the attached file, i used sort command to order it, and it works for examples1 and 2, but not for 3.

i guess i could use Abs cmd.....

I want to find the first term and the common difference of the arithmetic progression (un) with the following condintions u7 + u15 = 60 and (u4)^2 + (u12)^2 = 1170, but i can not the code. Please help me. Thank you.

given differntial equations with boundary conditions.



f(0)=k,  f'(0)=m, g(0)=1, 
f'(infinity)=1, g(infinity)=0

Dual solution exist in the given problem with shoot command i think. how we ll calculate them?


This thread stimulated me into playing with a few ideas,

I pieced together this little newton-raphson procedure...

> newton := proc( f, # the function x0 # the initial guess n, # step count limit tol # error tolerance)

local x, g, k;

g := D(f); # the derivative f’(x)

x[0] := evalf(x0); # initialize the iteration

for k from 1 to n do # loop for newton’s iteration

# Newton’s iteration formula

I can not get any subscription in the forum for any topic. Activating the subscription link directs me to the page:

Without the subscription validity.

The goal is to find generating function of Hermit

Before i try, i try the gFun's paper's example, got error

eq := {z^2*diff(y(z),z)+y(z)=z, (D(y))(0)=1};
Error, (in gfun:-borel) invalid argument, Y

When try Hermit,

eq := {t*diff(y(z),z$2)-z*diff(y(z),z)+y(z)=0,y(0)=1,y(-z)=y(z)};

but have unknown constant coefficient, do not know whether missing initial condition

I am trying to solve the system of equation by maple but its not giving the correct answer.

> e1 := (1/2)*alpha*alpha[3]^2*q[2]^3-V*alpha[4]*q[2]^2-V*alpha[0]+(1/2)*alpha*alpha[0]^2+(1/2)*alpha*alpha[1]^2*q[2]+(1/2)*alpha*alpha[2]^2*q[2]^2+(1/2)*alpha*alpha[4]^2*q[2]^4-V*alpha[2]*q[2]+C+alpha*alpha[3]*q[2]^2*alpha[1]+alpha*alpha[4]*q[2]^3*alpha[2]+alpha*alpha[4]*q[2]^2*alpha[0]+alpha*alpha[2]*q[2]*alpha[0];
> e2 := -gamma*alpha[2]*q[3]*q[2]^2+2*beta*alpha[4...

I was using "eliminate" to combine a set of equations removing some variables and getting a clearly wrong result. So I looked at the docs and the first example, and I find I don't understand its result either. Here is the example:

sol:=eliminate( {x^2+y^2-1=0, x^3-y^2*x+x*y-3=0}, x); # eliminate x from this set & leave one eqn of y

This returns:

First 787 788 789 790 791 792 793 Last Page 789 of 1528