Листинг 2 Создаем программу для вывода колебаний

global sigma r b

sigma=10;

r=28;

b=8./3;

x0=[10 20 30];

[T,Y]=ode45('lorenzo',[0 500],x0);

subplot(131)

plot(T,Y(:,1));

title('X ot T');

subplot(132)

plot(T,Y(:,2))

title('Y ot T');

subplot(133)

plot(T,Y(:,3))

title('Z ot T')

figure

plot3(Y(:,1),Y(:,2),Y(:,3),'k')

Приведите программу вывода на экран решений в виде колебаний инерционного нелинейного генератора.

dx/dt=mx+y-xz

dy/dt=-x

dz/dt=-gz+I(x)x²