MATLAB实现下列

Question

使用MATLAB解决图片内容，要求源代码，可以完整运行即可。使用Python语言若能实现也可以

补充帖子：下面是参考答案的代码，但是我无法运行，从答案的基础上修改也行

Answer 1

使用matlab完成

ff = odeset;
ff.Mass= @exc7eqm;
[t, x] = ode45(@exc7eqb, [0,10],[45;30;0;0],ff);
plot(t, x);


function M = exc7eqm(t, u)
M = [1 0 0 0; 0 1 0 0; 0 0 2 cos(u(1) - u(2)); 0 0 cos(u(1) - u(2)) 1];
end


function du = exc7eqb(t, u)
du = [u(3); u(4); 9.81 * sin(u(1)) - sin(u(1)-u(2))*u(4)^2; 9.81 * sin(u(2)) + sin(u(1)-u(2))*u(3)^2 ]
end

运行结果如下