Hi, I'm trying to create a hypersonic waverider through the tracing of streamlines through conical flow. This requires integrating the famous Taylor Maccoll equation to find the radial and angular velocity, at each ray between the shockwave and the generating cone. However, in my coding of the RK45, the angular velocity is supposed to equal zero at around 10 degrees for the shock wave angle i've chosen. however this is not the case and the integration well exceeds the required distance. Can anyone direct me as to where it is that i'm going wrong in my code?