Решить аналитическим и численным методом уравнения движения точки М с использованием современных вычислительных средств.

Код функции для вычисления скоростей и ускорений точки, радиуса кривизны траектории

function[Vx,Vy,V,ax,ay,a,at,an,ro]=kinema(x,y);

%кинематика точки

%входные параметры - х(t), y(t)- уравнения движения точки

Vx=diff(x,'t');

Vy=diff(y,'t');

V=sqrt(Vx^2+Vy^2);

ax=diff(Vx,'t');

ay=diff(Vy,'t');

a=sqrt(ax^2+ay^2);

at=(Vx*ax+Vy*ay)/V;

an=sqrt(a^2-at^2);

...