       Re: Solution of a non-linear ODE

• To: mathgroup at smc.vnet.net
• Subject: [mg81936] Re: Solution of a non-linear ODE
• From: Jens-Peer Kuska <kuska at informatik.uni-leipzig.de>
• Date: Sun, 7 Oct 2007 05:29:43 -0400 (EDT)
• References: <fe7hgc\$ok8\$1@smc.vnet.net>

```Hi,

FullSimplify @@ odesys

gives you a single second order equation for the
two variables u[t] and v[t] and you need a second
second order equation for v[t]. You can try

deq3 = FullSimplify @@ odesys;

odesys1 = {deq3, D[#, t] & /@ eq2}

and

DSolve[
odesys1, {u[t], v[t]}, t]

Regards
Jens

sigmundv at gmail.com wrote:
> Distinguished group members,
>
> I would like you to try the following:
>
> First define
>
> f[u_, v_] := 1
> eq1 = D[ArcTan[v'[t]/u'[t]], t] + 1/2 (D[Log[f[u[t], v[t]]], u[t]]
> v'[t]^2 - D[Log[f[u[t], v[t]]], v[t]] u'[t]) == 0;
> eq2 = u'[t]^2 + v'[t]^2 == 1/f[u[t], v[t]];
> odesys = {eq1, eq2};
> ic = {u == 0, v == 0, u' == 0, v' == 1};
>
> Then try to solve 'odesys', with the initial conditions 'ic', either
> using DSolve or NDSolve, or both. What do you get? For NDSolve I get
> the warning
>
> Solve::svars: Equations may not give solutions for all "solve" \
> variables.
>
> and the the input is returned unevaluated.
>
> When I use NDSolve on the initial value problem, the system runs
> several minutes, without returning any result.
>
> Can any of you explain this behaviour?
>
> When doing the same calculations with another similar software
> package, I get solutions easily. Therefore it seems strange that
> Mathematica does not return any solution.
>
> Best regards,
> Sigmund Vestergaard
>
>

```

• Prev by Date: Why is FindRoot[] HoldAll?
• Next by Date: Re: How to remove an artifact from a plot
• Previous by thread: Solution of a non-linear ODE
• Next by thread: RE: Re: BarChart - intensity of colour increases up the