% Моделирование стрельбы ядром из пушки
function yadro0
    close all
    %hold on
    alpha= pi/4; % угол стрельбы
    v0= 100;    % начальная скорость
    h= 0.001;   % шаг по времени точек в массиве xy
                % должен быть постоянным, иначе
                % анимация будет неправильной
    axes
    axis equal
    [t,xy]= ode45(@yrp,0:h:14,...
        [0 v0*cos(alpha) 0 v0*sin(alpha)]);
    comet(xy(:,1),xy(:,3));
end
function dxy=yrp(~,z)
    r= 0.3; % Радиус ядра
    mu= 0.3*1.2*pi*r^2/2;   % коэффициент трения для шара
    m= 30;  % масса ядра
    x= z(1); vx= z(2); y= z(3); vy= z(4);
    dxy= [vx;-mu*sqrt(vx^2+vy^2)*vx/m;...
        vy;-mu*sqrt(vx^2+vy^2)*vy/m-9.81];
end
