Doubt in the Impulse integral

I'm not get to resolve this equation, I don't know what is wrong...

I = @(Tc,t)(mdot.*(sqrt(((2.*y.*R_gas.*Tc)./(y - 1)).*(1 - ((Pe)/(Patm)).^((y - 1)/(y))))));

Q = dblquad(I,0,Tc,0,t);

[Ameer Hamza]

What is the error? What are the values of constants?

[Maria Angélica Pagane]

I don't know what is wrong this code. The values of constants fallow below:

y = 1.1332; % y = gama (Cp/Cv)

R = 8.314e+9; % Pa*mm3/mol*K

MM = 0.037192; % Kg/mol (propelant molar mass)

MM_gas = 0.60809; % Kg/mol (gas molar mass)

m = 0.2650; % Kg (propelant mass)

Patm = 0.10135; % MPa

Pe = 0.10135; % MPa (exhaust pressure)

Pc = [0,2000,6.895]; % MPa (combustion pressure)

Tc = [0,100,1750]; % K (combustion temperature)

t = [0,0.1,10]; % s (time)

Being,

R_gas = R / (MM_gas)

r = a*(Pc).^(n);

mdot = m./r ;

Answer 1

[Ameer Hamza]

Run the following code

y = 1.1332; % y = gama (Cp/Cv)
R = 8.314e+9; % Pa*mm3/mol*K 
MM = 0.037192; % Kg/mol (propelant molar mass)
MM_gas = 0.60809; % Kg/mol (gas molar mass)
m = 0.2650; % Kg (propelant mass)
Patm = 0.10135; % MPa 
Pe = 0.10135; % MPa (exhaust pressure)
Pc = [0,2000,6.895]; % MPa (combustion pressure)
Tc = [0,100,1750]; % K (combustion temperature)
t = [0,0.1,10]; % s (time)
R_gas = R / (MM_gas)
a = 1; % <---- undefined in your code
r = a*(Pc).^(m);
mdot = m./R ;

I = @(Tc,t)(mdot.*(sqrt(((2.*y.*R_gas.*Tc)./(y - 1)).*(1 - ((Pe)/(Patm)).^((y - 1)/(y))))));
Q = dblquad(I,0,Tc(end),0,t(end));

Also, use integral2 instead of dblquad

Q = integral2(I,0,Tc(end),0,t(end));