MathGroup Archive 2005

[Date Index] [Thread Index] [Author Index]

Search the Archive

A bug in the documentation?

  • To: mathgroup at smc.vnet.net
  • Subject: [mg56679] A bug in the documentation?
  • From: Andrzej Kozlowski <akoz at mimuw.edu.pl>
  • Date: Tue, 3 May 2005 05:26:30 -0400 (EDT)
  • Sender: owner-wri-mathgroup at wolfram.com

The following passage is taken from NDSolve Overview in Mathemaitca 5.1:

> Here is a simple DAE.
>
> In[119]:=
> NDSolve[{Derivative[2][x][t] + y[t] == x[t],
>    x[t]^2 + y[t]^2 == 1, x[0] == 0,
>    Derivative[1][x][0] == 1}, {x, y}, {t, 0, 2}]
>
> From In[119]:=
> NDSolve::ndsz: At t == 1.6656481721762024`, step size is effectively 
> zero; \
> singularity or stiff system suspected.
>
> Out[119]=
> {{x -> InterpolatingFunction[],
>    y -> InterpolatingFunction[]}}
>
> Note that while both of the equations have derivative terms, the 
> variable y appears without any derivatives, so NDSolve issues a 
> warning message.  When the usual substitution to convert to first 
> order equations is made, one of the equations does indeed become 
> effectively algebraic.

It is not true that both equations have derivative terms nor does the 
warning message seem to be of the kind that is referred to. If we 
change the equations and initial conditions to something like this:


NDSolve[{Derivative[2][x][t] + y[t] == x[t],
    Derivative[1][x][t]^2 + y[t]^2 == 1, x[0] == 0,
    Derivative[1][x][0] == 0}, {x, y}, {t, 0, 2}]


we get a message more like the one that is referred to in the 
documentation.


NDSolve::pdord: Some of the functions have zero differential order so 
the equations will be solved as a system of differential-algebraic 
equations.


NDSolve::ndsz: At  `1.0525780991716984 step size is effectively \
zero; singularity or stiff system suspected.


{{x -> InterpolatingFunction[],
    y -> InterpolatingFunction[]}}



Andrzej Kozlowski
Chiba, Japan
http://www.akikoz.net/andrzej/index.html
http://www.mimuw.edu.pl/~akoz/


  • Prev by Date: Re: arrange lists side by side
  • Next by Date: Re: debugging
  • Previous by thread: Re: arrange lists side by side
  • Next by thread: named pattern variable scoped as global, should be local