My equation is
dxdt=N0*sin(omega*t)*x*(1-x/K);N0=1:1:10;omega=1:1:10;
I want to plot 3-D for N0,omega,Xmax and Xmean.
I am trying with the following code but it does not work.
function RunOsciliationsky3DN0all= 1:1:10;N=length(N0all);omegaall= 1:1:10;M=length(omegaall);Pmax=zeros(N,M);Pmean=zeros(N,M);Pall=[Pmax,Pmean];x=size(Pall);for i=1:N for j =1:N [t,x]=ode45(@osciliation,[0 100],0.1,[],N0all(i),10,omegaall(j)); Pall(i,j)=x; endend[N0x,omegay]=meshgrid(N0all,omegaall);h=mesh(N0x,omegay,Pall);
1;
Best Answer