- #1

- 131

- 2

NDSolve[{-6250 f[t] + 0.025 p[t]^2 f[t] + 0.1 f[t]^3 + 3 q[t] f'[t] +

f''[t] == 0,

1.5230870989335428*^-35 p[t]^4 + q'[t] ==

0, -4.32*^36 q[t]^2 + 10/3 \[Pi]^2 p[t]^4 + 0.0125 p[t]^2 f[t]^2 +

0.025 (-62500 + f[t]^2)^2 == 0, f[0] == 1, f'[0] == 0,

q[0] == 6.899183459150538*^-13}, {f, p, q}, {t, 0, tf}]

q[0]/.%

After solving the equations, I tried check my solution. I found that q[0]=4.75446*10^-15 rather than the initial value 6.899183459150538*^-13 given in the code.I check q[0] because when I plot my solutions, q[t] is somehow constant at 4.75446*10^-15. So could anyone tell me what is going wrong in my code?

Thank you so much!