2 курс / Вычислительная физика / Задание 14 / Новый текстовый документ

init xz=1, yz=0, Vxz=0, Vyz=6, xm=1.5, ym=0, Vxm=0, Vym=3;
{System of difurs}
Vxz'=-mc*g*xz/(xz*xz+yz*yz)^3/2+(mm*g*sqrt(xm*xm-xz*xz))/(xm*xm-xz*xz)^3/2;
Vyz'=-mc*g*yz/(xz*xz+yz*yz)^3/2+(mm*g*sqrt(ym*ym-yz*yz))/(ym*ym-yz*yz)^3/2;
Vxm'=mc*g*xm/(xz*xz+yz*yz)^3/2+(mz*g*sqrt(xm*xm-xz*xz))/(xm*xm-xz*xz)^3/2;
Vym'=mc*g*ym/(xz*xz+yz*yz)^3/2+(mz*g*sqrt(ym*ym-yz*yz))/(ym*ym-yz*yz)^3/2;
xz'=Vxz;
yz'=Vyz;
xm'=Vxm;
ym'=Vym;
output xz,yz,xm,ym,xz,yz,xm,ym;