Загрузка данных


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);