Листинг 2

% Листинг программы для нахождения численного решения

% уравнения движения маятника и визуализации

% объявление глобальной переменной

global omega;

% Задание начальных условий х(0)=2, р(0)=3

x0=[2 2];

% нахождение значений х и р вывод на экран семейства кривых

% фазового портрета нелинейного маятника при w0=[1;4]

for omega=1:0.3:4

[T,Y]=ode45('mayat',[0:0.05:4*pi],x0);

% mayat – имя файла-функции, cодержащий определение

% функции; [0:0.05:4*pi]– вектор, определяющий интервал

% интегрирования; x0 – вектор начальных условий

hold on

% построение графиков р(х)

plot(Y(:,1),Y(:,2),'k');

grid on

% задание диапазона осей

axis([-5 20 -15 15]);

title('Фазовый портрет нелинейного маятника для omega=[0;4]');

xlabel('x');

ylabel('y');

end