Неработещ код на решение на задача на MatLab

[sbozhinov]

Здравейте,

Имам следната задача по диференциални уравнения, която трябва да реша на MatLab и съм я решил но не ми работи кода. Може ли, който разбира да ми каже къде е проблемът за да ми тръгнат решенията спрямо задачата по условие.

Благодаря предварително.

Задачата съм прикачил най-отдолу

Равновесните точки са (3,-3) и (3,2)

А това е моят код :

function phasemovie
tmax = 15;
a0 = 3;
b0 = -3;
a1 = 3;
b1 = 2;
x = 1 : 0.5 : 4;
y = -5 : 0.2 : 5;
[X,Y] = meshgrid(x,y);
hold on
plot(a0,b0, 'r*', a1,b1,'r*')
axis([0,5,-6,6])
deltax = -X+3;
deltay = Y^2+Y-6;
R = sqrt(deltax.^2 + deltay.^2);
quiver(X,Y,deltax./R,deltay./R,0.5,'k')
[x0,y0] = ginput;
plot(x0,y0,'g+')

function Z = ff(x,y)
Z = [-y(1)+3;y(2)^2+y(2)-6];
end
[T,Z] = ode45(@ff,[0,tmax],[x0,y0]);
for k=1:length(T)
xx(k) = Z(k,1);
yy(k) = Z(k,2);
plot(xx,yy,'r')
M(k) = getframe;
end
movie(M,2)
end