Re: Is this so hard?
- To: mathgroup at smc.vnet.net
- Subject: [mg95270] Re: Is this so hard?
- From: "Steve Luttrell" <steve at _removemefirst_luttrell.org.uk>
- Date: Wed, 14 Jan 2009 05:50:53 -0500 (EST)
- References: <200901111138.GAA12346@smc.vnet.net> <gkgq3p$4a2$1@smc.vnet.net> <gkh0vt$6ia$1@smc.vnet.net>
Try solving it numerically.
t0 = 1.12;
sol = NDSolve[{x''[t] == 1/x[t]^2, x[0] == -1, x'[0] == 0},
x, {t, 0, t0}]
Plot[(x /. sol[[1]])[t], {t, 0, t0}]
This gives the following error message:
"NDSolve::ndsz: At t == 1.1107203590726262`, step size is effectively zero;
singularity or stiff system suspected."
There is a singularity where x[t] goes to zero, which makes the 1/x[t]^2
term in the ODE blow up.
--
Stephen Luttrell
West Malvern, UK
"cool-RR" <ram.rachum at gmail.com> wrote in message
news:gkh0vt$6ia$1 at smc.vnet.net...
>
>> You didn't want another initial (or boundary) condition for this 2nd
>> order ODE?
>>
>
> You're right.
>
> DSolve[{x''[t] == 1/x[t]^2, x[0] == -1, x'[0] == 0}, x, t]
>
> This still gives no answer.
>
> Any ideas?
>