Курсовая работа: Автоколебательная система. Волны пластической деформации

axis([0 5.5 0 0.6]);

pause;

clear all;

alpha=1;

beta=10;

gamma=1;

f=@(t,y)[-y(1)+alpha*y(2);beta*(1-y(2)-gamma*y(1)*y(2))];

for i=0:1/5:1,

for j=0:1/5:1,

[T,Y]=ode45(f,[0 100],[i j]);

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

hold on;

end

end

hold off;