coupled Ordinary Differential Equations

Can someone help me solve the following set of coupled equations

Inline form:

Display form:

With initial conditions ; ; r ranges from 0 to 30. Take .

needs to be interpolated with calculated pressure at every iteration (which I can handle with).

[Torsten]

What problems do you encounter when you try to use ode45 to solve ? "interp1" can be used for interpolation therein.

https://uk.mathworks.com/help/matlab/ref/ode45.html

https://uk.mathworks.com/help/matlab/ref/double.interp1.html

[Sam Chak]

Hi @Arunkarthiheyan

If r can range from 0 to 30, then singularity will occur at . Since both the numerator and denominator have division-by-zero terms, then will become an indeterminate form. The response of the variable m will depend on the behavior of the function .