我把第一个式子改写了一下
v=0.5;
dfun=@(z,x)[((1-v)*((cos(x(1))+1).^2-1)-v)/x(2);-cot(x(1));];%x(1)=theta,x(2)=R
ode15s(dfun,[0 10],[1 1])
我把第一个式子改写了一下
v=0.5;
dfun=@(z,x)[((1-v)*((cos(x(1))+1).^2-1)-v)/x(2);-cot(x(1));];%x(1)=theta,x(2)=R
ode15s(dfun,[0 10],[1 1])