использование ODE45 в MATLAB для решения трех PDE

Question

Я пытаюсь решить уравнения 3, 10 и 11 из следующего источника: http://web.math.ku.dk/~moller/e04/bio/ludwig78.pdf

Ниже приведен код, который я создал в MATLAB, однако он не производит график, отображающий три решения на одном графике.

%% x=(S,N)^T
%%initial value:
x0= [1;1/2];
vopt= odeset('RelTol' ,1e-3, 'AbsTol', 1e-3, 'NormControl' , 'on', 'InitialStep',0.1,'MaxStep',0.5);
beta=20;E=1;r=0.1;K =4*S;P=10;N0=S/2;alpha=1;

%%define the right hand side:
f = @(t,x) r.*x(1).*(1-x(1)./K)-beta.*(N0.*N0./(alpha.^2-N0.^2));
        r.*x(2).*(1-(x(2)./K).*(K./x(3)));
        r.*x(3).*(1-x(3)./K)-P.*(x(1)./x(2))

[t,y] = ode45(f,[0 30],x0,vopt);
figure(1);
plot(t,y(:,1),'N',t,y(:,2),'S',t,y(:,3),'E');
legend('Population of Budworms','Available foliage','Energy Reserve');