我虚脊世把第一个式差肢子改写了一下野槐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)=Rode15s(dfun,[0 10],[1 1])
