I am taking a class where we have to solve second order differential equations (symbolically) for spherical and cylindrical coordinates. The equations are sometimes written
D(r^2*D(f(r)) (all are differented with respect to r, but it has r^2 between first and second differential.
Can Maple solve something in that form?

Hello everybody!
I again have a problem with "Dsolve numeric".
I solve a system of linear ODEs with boundary and continuity conditions.I got an error, which is associated with any internal procedure of Maple. For the sake of simplification, I reduced the order of equation from two fourth-order ODEs to the eight equations of the first order, but the same error is outputed.

Download 2847_problem-1-1.mwsView file details
Please help me clarify it.
Thank you

Hi everybody,
As a part of some general problem, I need to solve a system of non-linear ordinary differential equations.I try to get a numerical solution, so I use the "dsolve numeric" command. In same cases, the Maple print out two errors that I didn't find in his help instruction.
The first one is:
error, (in dsolve/numeric/bvp) cannot determine if this expression is true or false: 20*nerr

I have a large system of equations that takes a good hour to solve.
Each time I save and reopen the worksheet there is another lengthy recalculation on execution after which other operations are carried out.
I know how to save individual variables for reuse but what I would like to do is save the session so that the I can pick up from where I left off. Is this possible?

**UPDATE!!! ** I found that my blogpost appers in Maple reporter. You should know, that this blogpost was just a question about how to solve my problem better. Alec Mihailovs gave me and aswer and his solution of testing nested objects is much better. If you still want to use my piece of code, you should remove "set object" form "convertAMVStolist" procedure and replace all occurrences of "hastype" with "type" (as mentioned in comments below).

This piece of code should be able to test two objects (not of every type) for equivalence (like *testeq()* does). The benefit is, that it should be able to test also nested objects. Is there any other and more simple way how to do that? How to test nested objects in sets?

Hello!
I try to solve the eq1, eq2, eq3, and eq4.
I simplified the 4 equations to rs1 and rs2
to get the solution t3 and t5.
Then, I check the solution by subsituting Eq(14)
into Eq(10). The result is right. The rs1 equation
was found by numerical tests.
However, it is not correct for the eq1, eq2, eq3,
or eq4. I do not understand why it is not correct.
Please tell me the reason and how I can solve and obtain
correct solution. The worksheet can be found in the

File Manager
Thank you for your time and help.

October 25 2006
PyWiz 8
Hi
Please see Spiderflygrid2.mws
A spider is on the floor of a cubical room of side 1. Its position is given. Where must the fly position itself to be as far as possible from the spider?
My attempt minimizes the square of the twelve possible distances. This is m.
m:=min(seq(cd[i], i=1..12));
Partial derivatives of m are found with respect to xf (the x coord of the fly)
eq1:=diff(m,xf);
... & with respect to yf (the y coord of the fly)
eq2:=diff(m,yf);
These two equations are then solved simultaneously.
solve({eq1,eq2});
Maple gives an error message - is this due to the discontinuities in the expression m? Or (more likely) a fault in my Maple code?

Hello,
I am quite new to Maple and I am trying to solve a system of 3 non-linear equations in 3 variables (A,zp,zh). I have tried the solve command, but after waiting few hours I have decided to stop the computation. Could someone help me with this?
The system is as follows
eq1 := 1/(-138.8888889/(24000000-50/3*zh)^1.5-.50/zh^1.5-122.7222223/(-8640000+47/3*zh)^1.5)*(5.465600000*A/(1+zp)^.5e-1
+(-16.66666667/(24000000-50/3*zh)^.5+1.0/zh^.5+15.66666667/(-8640000+47/3*zh)^.5)*(.3e-1+.3*A))=0;
eq2 := .4879000781e-6/A*(1+zp)^2.05*((-1952000.00*A/(1+zp)^1.05+50/3)*(.3e-1+.3*A)-10.93120000*A/(1+zp)^.5e-1

Hi,
I'm trying to get maple to output the numerical answer, from a call to dsolve, to an external txt file; but I'm having a bit of trouble. I want to do this so that I may then plot the results in an external program like easy-plot, or just analyse them.
I'm using a "read" command in Maple to read in a procedure from a txt file (shown below), which I then execute in Maple:
#Utilise prime notation for derivatives wrt t:
with(PDEtools, declare):
with(plots):
declare(y(t), prime=t, quiet):
#Also shorten output notation:
alias(y=y(t),y0=y(0),yp0=D(y)(0)):
#define the ODE:
ODE:=diff(y,[t$2])-(1-y^2)*diff(y,t)+y=0;

October 17 2006
bjain 8
Download 569_ODE-25 Reactions.doc
Hi,
I want to solve 25 set of differential equation for unknown variable and contants.
Problem Statement: Solve the following set of differential equations for concentrations terms in terms of CA0 and t
CA, CB, CC, CD, CE, CF, CG, CH, CI, CJ, CK, CL, CM, CN, CP are the variables and k1 to k25 and CAO are constants
Initial Values: t=0 CA=CA0 and CB=CC=CD=CE=CF=CG=CH=CI=CJ=CK=CL=CM=CN=CP =0
How do I proceed for the same?
Regards,
Bikash

Hi,
I am trying to use fsolve to find the roots of p non-linear equations. These vary from 2 dimensional to, so far, 5 dimensional.
I'm really trying to minimize the log likelihood of an ARCH time series:
L := (t) -> 1/2*T*log(2*Pi)+1/2*Sum(log(h(t)),t = p .. T)+1/2*Sum(y[t]^2/h(t),t = p .. T):
where
h := (t) -> alpha[1] + alpha[2]*y[t-1]^2 + alpha[3]*y[t-2]^2 ... + alpha[p]*u[y-p]^2:
and
y[t] are observed data.
I've used the p partial derivatives to solve for dL/dalpha=0 using fsolve and intervals alpha=0..1. Occasionally, it works fine and estimates the parameter vector alpha reasonably. However, sometimes it takes about 30-60 seconds and just returns the equation and no answer:

October 10 2006
Rick 28
Hi,
I have a question about Gröbner bases and using the Rosenfeld_Grobner algorithm in diffalg. What is best to use? If the equations only contain time differentiated variables, is it better to differentiate the equations "by hand" first and let the Gröbner package deal with the problem or is it better to just plug the whole thing into diffalg from the start?

Given N-1 independent (not necessarily linear) equations in N variables, how do I tell Maple to find solutions in terms of a specific variable?
For example:
with(RealDomain);
eq1:=x-1=a*cos(theta);
eq2:=x+1=a*sin(theta);
I want to find a and theta in terms of x:
a=sqrt(2*x^2+2), theta=arctan((x+1)/(x-1))
solve({eq1,eq2},[a,theta]) appears to give me solutions based on elimination of x.

Dear All,
I'm working on a problem in differential equations where I need to compute the n-th power of a differential expression. Can Maple be taught to differentiate polynomial functions without carrying out the multiplications / simplifications that it usually does? For example:
diff(-2*c/x^3,x); gives the answer 6*c/x^4.
I want it to return: (-3)(-2)*c/x^4.
Ideally, I also do not want it to simplify sums of such results.
Thank you very much, I appreciate your help!
Andrea Bruder