p >Hello. I'm still kind of new to the practice of using a computer program to manipulate equations. I knew that one way to define an ellipse is as the locus of points P such that for 2 fixed points F1 and F2 the sum PF1 + PF2 is constant. I also know that in Cartesian coordinates, we get a formula of the form Ax2 + By2 = 1. If we set up our coordinate system such that F1 and F2 are on the x axis each a distance c from the origin and the maximum width of the ellipse is 2a, let's try to figure out A and B using MAPLE.

**> **

Hi, I have a system of 3 coupled equations (3 unknowns) with complex numbers. I solved them with solvefor. I'm interested in the modulus (|a+ib|^2 = a^2 + b^2 ) and phase of each solution.
Because of that I used polar and abs (I have more interest in the modulus). In each case I ran simplify, despite of that the solutions are huge.
I would like to work with them in some spreadsheet (in order to set the parameters, evaluate the solutions, compare them with some experimental data and do some plotting).
Given that the solutions are huge I can just write them by hand where I want. In order to automatize that I've tried to use the VisualBasic command to cut&paste the formulas in VBA. The problem is that they are so long that the source code just gets broken.

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.

I'm trying to solve a system of equations in maple and I have no idea why it's doing what it's doing but it certainly isn't working. Anybody have an idea what's going on?

View 2503_numbers2.mw on MapleNet or

Download 2503_numbers2.mwView file details
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;

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:

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?