funtion sim_system tspan = 0:0.001:50; z0 = [1,0,sqrt(2),0]; [t,z]=ode45(@problem,tspan,z(0)); plot(t,z(:,1)); figure(); end
function zdot = problem(t,z) zdot = zeros(4,1); zdot(1)=z(2); zdot(2)=((2*sin(z(1))+(sin(z(1)-z(3))*z(4)*z(4))+(cos(z(1)-z(3))*sin(z(1)-z(3))*z(2)*z(2))-(cos(z(1)-z(3))*sin(z(3))))/((cos(z(1)-z(3))*cos(z(1)-z(3)))-2)); zdot(3)=z(4); zdot(4)=((-2*sin(z(1)-z(3))*z(2)*z(2))+(2*sin(z(3)))-(cos(z(1)-z(3))*sin(z(1)-z(3))*z(4)*z(4))+(2*cos(z(1)-z(3))*sin(z(1)))/((cos(z(1)-z(3))*cos(z(1)-z(3)))-2));
Best Answer