Hi guys,
I was just coding to plot an ODE, but ran into some problems dispalying the result
The code are as follows
function dxdt = ODE(t, x, alpha, beta)dxdt = diag(alpha - beta*[x(1);x(2);x(3)])*[x(1);x(2);x(3)];for n=1:200p=.01*n;alpha=[.1; .5; .3];beta=[.2 .3 .4; p 1.5 .6; .2 .5 .8];time = [0 100];initial=[0.5; 0.2; 0.3];[t,y] = ode45(@(t,y) ODE(t,y,alpha,beta), time, initial);F1=y(100,1);F2=y(100,2);F3=y(100,3);endfigure(1)plot(p,F1);hold onplot(p,F2);plot(p,F3);hold off
Could you tell me what's off?
Many thanks in advance
Best Answer