大家看看这个编程，用的是ODE45，不知道什么地方出了问题，有错误，我是完全按照ODE45的格式做的。用的是夏的那篇文章（附件）
M文件：
function xp=led3(t,x)
xp=zeros(2,1);
xp(1)=x(2);
xp(2)=-((3.5+x(1))/(cos(t*pi/180))^2/((1.53419*cos(t*pi/180)-sqrt(cos(t*pi/180)))/(1.53419*sin(t*pi/180)-sqrt(1-cos(t*pi/180)))+tan(t*pi/180)));

end
ODE45:
=ode45('led3',,0)
plot(t,x(:,1),'r-o',t,x(:,2),'b-*');
legend('y1','y2');
xlabel('t');
xlabel('y');

• 近朗伯光型LED透镜的光学设计.pdf