MathGroup Archive 2004

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

Search the Archive

How to express ODE?

  • To: mathgroup at
  • Subject: [mg52817] How to express ODE?
  • From: Adam Getchell <agetchell at>
  • Date: Mon, 13 Dec 2004 04:24:15 -0500 (EST)
  • Sender: owner-wri-mathgroup at

Still (unfortunately) not getting another ODE to work, even using Dr. 
Bob's example. I want to express this equation:

m = 10^16;
mp = 1.2211*10^19;
G = mp^(-2);
V[\[Phi]_] := (1/2)*m^2*\[Phi][t]
\[Rho][\[Phi]_] :=
  Derivative[1][\[Phi]][t]^2/2 +
H := Sqrt[((8*Pi*G)/3)*

Inflaton :=
  Derivative[2][\[Phi]][t] +
     t] + D[V[\[Phi]], \[Phi]]

That is, phi double-dot of t + 3 H phi dot of t + V' of phi, where phi 
is a function of t; V is a function of phi; H is a function of rho; rho 
is a function of phi dot and V.

However, solving Inflaton for phi as a function of t:

f = \[Phi] /. First[NDSolve[
     {Inflaton == 0,
        t] == 0}, \[Phi],
     {t, 0, 100*m}]]

Generates NDSolve::overdet and ReplaceAll::reps.

I *think* part of my problem is that dV/d(phi) is not getting properly 
evaluated (it returns 0 when I look at it explicitly).

I thought that Derivative[1][\[Phi]][x] would give the above, but not 
the way I'm doing it.


  • Prev by Date: Re: How to evaluate Exp[I Pi(1+x)]?
  • Next by Date: Re: Solve bug?
  • Previous by thread: Re: rotatelabel in linux
  • Next by thread: Re: How to express ODE?