Taylor-Maccoll equations for axis-symmetric conical flow using the ode45 solver

I wrote the differential equations a system of equations. I have reproduced the function as follows - function dvrdtheta= TMequation(theta,vr) gamma=1.4; dvrdtheta=zeros(2,1); dvrdtheta(1) = vr(2); dvrdtheta(2) =(((gamma-1)/2)*(1-vr(1).^2-vr(2).^2)*(2*vr(1)+cotd(theta).*vr(2))-vr(1).*vr(2).^2)./(vr(2).^2-((gamma-1)/2)*(1-vr(1).^2-vr(2).^2)); I just don't know where I got this wrong. Please help me fellas!

[Rajan Bhandari]

I also have been working to solve this equation but I got the problem with integration tolerance. What is your problem actually?

Answer 1

[Rohan Kokate]

Attaching the link to the code which solves the Taylor Maccoll equation numerically. Hope this helps!

https://www.mathworks.com/matlabcentral/fileexchange/84490-taylor-maccoll-solution-for-supersonic-flow-over-a-cone