clc;
clear;
close all;
% Решение дифференциального уравнения
[t, y] = ode45(@primer1, [0 10], [0.01 1]);
% График перемещения
figure(1);
plot(t, y(:,1));
grid on;
xlabel('t, c');
ylabel('x, м');
title('Зависимость перемещения от времени');
% График скорости
figure(2);
plot(t, y(:,2));
grid on;
xlabel('t, c');
ylabel('V, м/с');
title('Зависимость скорости от времени');
% Максимальное перемещение
A = max(abs(y(:,1)));
fprintf('Максимальное перемещение A = %.6f м\n', A);