Hi, I need to solve 3 coupled first order ODE's using NDSolve (numerical solution).(adsbygoogle = window.adsbygoogle || []).push({});

This is the code I have used ;

NDSolve[{u'[t] == ((1 - (u[t]/t^2))/(-3 - (u[t]/t^2) - (v[t]/t^2))),

v'[t] == (((u[t]/t^2) - (v[t]/t^2))/(-3 - (u[t]/t^2) - (v[t]/t^2))),

w'[t] == (v[t])/(-3 t^2 - (u[t]) - (v[t])),

u[0] == 0., v[0] == 0., w[0] == 0.}, {u, v, w}, {t, 0.1, 1}]

where u,v and w are functions of t and at t=0, u=v=w=0. However I am getting the following error "Infinite expression 1/0.^2 encountered", but I can't see where the division by zero occurs. Originally t was from 0 to 1, but I tried making it from 0.1 to 1 to ensure no division by zero's were taking place.

Any ideas on how to fix the code? or an alternative way of solving these ode?

Regards

**Physics Forums - The Fusion of Science and Community**

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

# NDSolve coupled ode

Loading...

Similar Threads - NDSolve coupled | Date |
---|---|

Mathematica Rescaling equations in Mathematica | May 5, 2017 |

Matlab Coupled System of ODE's with algebraic coefficients | Sep 18, 2016 |

Mathematica NDSolve different solutions same eq's | Jul 23, 2015 |

Solving ODE numerically in Mathematica - I get 'ndnum' error | Mar 17, 2015 |

Mathematica NDSolve initial condition | Mar 5, 2015 |

**Physics Forums - The Fusion of Science and Community**