分别算w=5 和 w=20的时候 求代码求图

分别算w=5 和 w=20的时候 求代码求图

clc
w=5;
f=@(t,y)([w*(10+20*cos(y(3))-y(1))/sqrt((10+20*cos(y(3))-y(1))^2+(20+15*sin(y(3))-y(2))^2);
w*(20+15*sin(y(3))-y(2))/sqrt((10+20*cos(y(3))-y(1))^2+(20+15*sin(y(3))-y(2))^2);
1/sqrt(400*(sin(y(3)))^2+225*(cos(y(3)))^2)]);
ode45(f,[0 1],[0 0 0])
第二个程序把w值改一下即可