Pc(2): out of bound 1 error

[Bala Vignesh]

I am trying to solve ODE equation inside a for loop for a given mdot, Ab and
time step. I am getting error .. can anyone help me please

Code:
rhop=1765;
a_burn=2.6292e-5;
n=0.33;
a=1577.1;
Ath=0.008167;
a=1577.1;
b=9e-7;
gamma=0.6465;
Vc=0.0943;
pc0=0;
%
%Read excel files
A=xlsread('Ignition transient.xlsx');
m=A(:,7);
tspan=A(:,4);
Area=A(:,6);
for i=2:40
  mdot=m(i)
  Ab=Area(i)
[t,Pc]= ode45(@(t,Pc)
(mdot+(rhop*Ab*a*Pc^n)-(Pc(i)*Ath/(a+(b*Pc))))*(gamma^2*(a+(b*Pc))^3/(Vc*(a-(b*Pc)))),
tspan, 0);
end

Error message:
error: Pc(2): out of bound 1
error: called from
    ignitionTransient>@<anonymous> at line 29 column 69
    starting_stepsize at line 46 column 5
    ode45 at line 192 column 25
    ignitionTransient at line 29 column 7



--
Sent from: https://octave.1599824.n4.nabble.com/Octave-General-f1599825.html