r0 = 10；z0 = 1；%自己改
w0 = 0.0942？
a = 0.875
[r，z]=ode45(@(r，z)-tan(w0*tan(w0*r/a))，[-r0，r0]，z0)；
繪圖(r，z)
註:夾角φ(或圖上的θ)=-(正水平方向到正切線方向的夾角)，這個夾角的切線就是這個函數的導數，所以這個方程就轉換成下面的形式DZ/Du = Tan(-φ)=-Tan(W0 * Tan(W0 * R/A))，所以可以直接用ode45求解。