以下代码提供了解决圆锥超音速流动Taylor Maccoll方程的解决方案。用户需指定自由流马赫数(M)、半锥角(theta_cone)和比热比值(g)。该求解器调用四个子程序来求解整个流场,并输出所有关键的流变量。参考资料:空气动力学基础,JD Anderson。