Hi, I am trying to use NDSolve in Mathematica to solve a set of differential-algebraic equations:(adsbygoogle = window.adsbygoogle || []).push({});

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!

**Physics Forums | Science Articles, Homework Help, Discussion**

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Mathematica Mathematica NDSolve initial condition

Loading...

Similar Threads for Mathematica NDSolve initial |
---|

Mathematica Storing Mathematica output |

Mathematica Mathematica to MATLAB |

Mathematica Cannot do the integral of the Hyper-geometric function? |

Mathematica Cannot Plot This Function in Mathematica |

Mathematica While Loop in Mathematica |

**Physics Forums | Science Articles, Homework Help, Discussion**