[Date Index]
[Thread Index]
[Author Index]
Mathematica question : Stepsize control
*To*: mathgroup at smc.vnet.net
*Subject*: [mg83357] Mathematica question : Stepsize control
*From*: "Yannick Bartocci" <yannick.bartocci at gmail.com>
*Date*: Sun, 18 Nov 2007 04:48:07 -0500 (EST)
Hy,
I would like to ask a question on your forum bout some Mathematica code I
cannot solve.
Here is a Riccati type equation I want to solve numerically :
tend = 0.5*10^4;
ktest = 10^(-8);
h[t_, k_] = k + 10^(-6)*Cos[[Pi]/(2*tend)*t];
equ = {f'[t]+(f[t])^2+h[t,ktest] == 0, f[1] == 0.1};
{ndsol, steps} = Reap[NDSolve[equ, f, {t, 1, tend}, MaxSteps -> \[Infinity],
StepMonitor :> Sow[{t, f[t]}]]]
this code gives diverging solution for f[t], but note that if
tend = 10^2
the solution is not diverging, which is what I want to reproduce for tend =
0.5*10^4 or larger.
Actually the system I want to solve uses a more complicated h[t,ktest] as
Interpolating function input which is slowly varying in the beginning but
rapidly oscillating in the end. Initially it behaves like the system
described above and gives NDSolve message ::ndsz.
This clearly is a problem of stepsize control but which has to be adapted on
the fast or slow behavior of the input function h[t,ktest].
Can anybody help me to find out how to control stepsize for this problem?
The mathod in NDSolve is by default Automatic.
Thank you very much,
Yan
Prev by Date:
**Re: Trouble with Rotate**
Next by Date:
**Plot to Plot3D problem**
Previous by thread:
**Re: trying to Import[] USGS-resident DEM data**
Next by thread:
**Re: Mathematica question : Stepsize control**
| |