I have this code which I am getting a warning message when I am trying to run it ,,could anyone help me please??
function RunlogOscilnumerfishfixedn0omega
omega=1;
N0=1;
k = 10;
A = 1;
p0 = 0.1;
tspan=(0:0.01:4);
[t,p] = ode45(@logOscilnumeric,tspan,p0,[],omega,k,N0);
figure(1)
plot(t,p)
P = @(T) interp1(t,p,T);
f = @(t) ( ( A.*( ( N0.* (sin(omega.*t)).^2 .*(1-(2.*P(t)./k))+(omega.*cos(omega.*t) ) ).^2 ) ./( (N0.^2).*((sin(omega.*t)).^4).*((P(t)-P(t).^2./k)).^2 )) ) ;
I1 = integral( f1, 0.01,1,'ArrayValued',true)./1
I2 = integral( f1, 0.01,2,'ArrayValued',true)./2
I3 = integral( f1, 0.01,3,'ArrayValued',true)./3
I4 = integral( f1, 0.01,4,'ArrayValued',true)./4
I=[I1,I2,I3,I4];
T=[1,2,3,4];
figure(2)
plot(T,I)
1;
% function dpdt = logisticOscilnumerical(t,p,omega,k,N0)
% dpdt = N0*sin(omega*t)*p*(1-p/k);
% end
Kindly save as the following codes in a separate file named with name logisticOscilnumerical
function dpdt = logisticOscilnumerical(t,p,omega,k,N0)
dpdt = N0*sin(omega*t)*p*(1-p/k);
After that save as the following codes in another file with any name and run, you can get the answer
%function RunlogOscilnumerfishfixedn0omega
omega=1;
N0=1;
k = 10;
A = 1;
p0 = 0.1;
tspan=(0:0.01:4);
[t,p] = ode45(@logisticOscilnumerical,tspan,p0,[],omega,k,N0);
figure(1)
plot(t,p)
P = @(T) interp1(t,p,T);
f = @(t) ( ( A.*( ( N0.* (sin(omega.*t)).^2 .*(1-(2.*P(t)./k))+(omega.*cos(omega.*t) ) ).^2 ) ./( (N0.^2).*((sin(omega.*t)).^4).*((P(t)-P(t).^2./k)).^2 )) ) ;
I1 = integral( f1, 0.01,1,'ArrayValued',true)./1
I2 = integral( f1, 0.01,2,'ArrayValued',true)./2
I3 = integral( f1, 0.01,3,'ArrayValued',true)./3
I4 = integral( f1, 0.01,4,'ArrayValued',true)./4
I=[I1,I2,I3,I4];
T=[1,2,3,4];
figure(2)
plot(T,I)
1;
- Badges
- Science topic