Re: NDSolve problem
- To: mathgroup at smc.vnet.net
- Subject: [mg31817] Re: NDSolve problem
- From: Roland Franzius <Roland.franzius at Uos.de>
- Date: Wed, 5 Dec 2001 06:51:39 -0500 (EST)
- Organization: Uni Osnabrueck
- References: <email@example.com> <firstname.lastname@example.org>
- Sender: owner-wri-mathgroup at wolfram.com
If the system of ODE's has a critical point somewhere, say at x=a, it should
always be possible, to extract the singularilty by the ansatz of critical
y_i = (x-a)^(s_i) z_i(x), i=1,..n
in such a way that the resulting system for z is smooth or L^1 or something like
that. At least one solution should be well behaved at this point then.
Jens-Peer Kuska wrote:
> your system of equations has a singularity at this point
> and one of the components goes to infinity. The stepsize
> control detect the singularity and reduce the stepsize
> until it is to small.
> You should care full analyse the differential equation
> if the singularity is not expected/obvious.
> With an ordinary differntial equations solver it is
> impossible to jump over the singularity, because
> the most ode solver use a polynomial approximation
> to the solution and you can't approximate 1/x with
> a polynomial of positive powers.
> deniz.seker at arcelik.com.tr wrote:
> > Hi,
> > I have a problem with NDSolve. I'm trying to solve 2 coupled differential
> > equations. But dur,ng the solution, I saw a message that,"At t=.... step
> > size is effect,vely zero, singularity suspected." Then, I'm trying to
> > decrease StartingStepSize option up to 10^(-14) or increase MaxSteps to
> > 100000. But, for these two cases, I couldn't get solutions for 6 hours run.
> > What should I do?
> > Best regards,
> > Deniz
> > PS : If anybody wants to see the notebook, I can send via e-mail.
Prev by Date:
WTB: Control Systems Professional for Mathematica 2.2 for the Macintosh
Next by Date:
scope all wrong? in Mathematica 4.1
Previous by thread:
Re: NDSolve problem
Next by thread:
Mathematica for Mac OS X now shipping