问题描述
我如何在以下两个体中加入对电荷的时间依赖性 q1 = q2 = q0(2-exp(-t)) 问题？提前致谢
ini=[r1(1) r1(2) r1(3) r2(1) r2(2) r2(3) vx1 vy1 vz1 vx2 vy2 vz2]; % initial conditions
options= odeset('AbsTol',1e-12,'RelTol',1e-12);
F = @(t,r) [r(7);r(8);r(9);r(10);r(11);r(12);...
((q1*q2/m(1))*((r(1)-r(4))/((r(1)-r(4))^2+(r(2)-r(5))^2+(r(3)-r(6))^2)^(3/2)));...
((q1*q2/m(1))*((r(2)-r(5))/((r(1)-r(4))^2+(r(2)-r(5))^2+(r(3)-r(6))^2)^(3/2)));...
((q1*q2/m(1))*((r(3)-r(6))/((r(1)-r(4))^2+(r(2)-r(5))^2+(r(3)-r(6))^2)^(3/2)));...
((q1*q2/m(2))*((r(4)-r(1))/((r(1)-r(4))^2+(r(2)-r(5))^2+(r(3)-r(6))^2)^(3/2)));...
((q1*q2/m(2))*((r(5)-r(2))/((r(1)-r(4))^2+(r(2)-r(5))^2+(r(3)-r(6))^2)^(3/2)));...
((q1*q2/m(2))*((r(6)-r(3))/((r(1)-r(4))^2+(r(2)-r(5))^2+(r(3)-r(6))^2)^(3/2)))];
[t,r]=ode45(F,tspan,ini,options);
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）