problema con scilab!!

[phytona90]

Tengo que hace 1 programa con SCILAB que calcule la variación de Tº de la llama maxima durante la combustión de gas naturral en un proceso adiabatico en funcion de la concentración de metano,sabiendo que la concentración de este varía entorno al 95-63% del gas y etano 3-33%.Y HACER UNA GRAFICA que los confronte.
He plantEado lo siguiente pero me da error en x=fsolve y no se resolverlo.espero q alguien me eche 1 cable y veis algun error de planteamiento y demas comentarmelo por favor.gracias

Código scilab:

Qcomb=(-212.798);
Qcomb2=(-372.820);
CH4=95; a1=3.381; b1=18.044*10^(-3);c1=-4.3*10^(-6);// Cantidad de entrada y constantes de metano
C2H6=3 ;a2=2.247; b2=38.201*10^(-3);c2=-11.049^10^(-6);// Cantidad de entrada y constantes de etano
N2=160;a3=6.524;b3=10.396*10^(-3);c3=-0.001*10^(-6);//Cantidad de entrada/salida de N2
O2e=42;a4=6.148;b4=3.102*10^(-3);c4=-0.923*10^(-6);// Cantidad de entrada y constantes de OXIGENO
CO2=98;a5=3.214;b5=10.396*10^(-3);c5=-3.545*10^(-6);// Cantidad de salida y constantes de CO2
H2O=199;a6=7.256;b6=2.298*10^(-3);c6=0.283*10^(-6);// Cantidad de salida y constantes de H2O
 
 
t0=298;
i=0
 
scilab.org/product/dic-mat-sci/M2SCI_doc.htm">while i<=32//Condicion impuesta porque CH4 VARÍA DEL 95% AL 63%
    CH4=CH4-i;
    C2H6=C2H6+i
    CO2=(CH4)+(2*C2H6)
    H20=(2*CH4)+(3*C2H6)
    O2s=42-((CH4*2)+(C2H6*(7/2)))
    CONC=[]//Inicializamos vectores de % de ch4 y tº vacios para recoger una lista de valores y luego ilustrarlos
    T_LLAMA=[]
    //Calculo de la entalpia de entrada 
    Cp1=(a1+(b1*t0)+(c1*(t0^2)))*10^(-3); Q1=CH4*Cp1// CP del metano y calor 
    Cp2=(a2+(b2*t0)+(c2*(t0^2)))*10^(-3);Q2=C2H6*Cp2// CP del etano y calor
    Cp3=(a3+(b3*t0)+(c3*(t0^2)))*10^(-3);Q3=Cp3*O2e// CP del O2 y calor
    Cp4=(a4+(b4*t0)+(c4*(t0^2)))*10^(-3);Q4=Cp4*N2// CP del N2 y calor
    He=Cp1+Cp2+Cp3+Cp4
    //Calculo de la entalpia de salida
    Cp5=(a3+(b3*t)+(c3*(t^2)))*10^(-3); Q5=N2*Cp5// CP del N2 y calor
    Cp6=(a4+(b4*t)+(c4*(t^2)))*10^(-3); Q6=O2s*Cp6// CP del O2 y calor
    Cp7=(a5+(b5*t)+(c5*(t^2)))*10^(-3); Q1=CO2*Cp7// CP del CO2 y calor
    Cp8=(a6+(b6*t0)+(c6*(t^2)))*10^(-3); Q1=H2O*Cp8// CP del H2O y calor
    //calculo de la entalpia de reacción
    AHr=((CH4*Qcomb)+(C2H6*Qcomb2))
    deff("fc=llama(t)",'fc=(Hs-He+AHr)'))
    x=fsolve(298,llama)
    printf(("La respuesta es :%2.3f"),x)
    disp("Para este % de CH4:",CH4)
    T_LLAMA=T_LLAMA+x
    CONC=CONC+CH4
    i=i+1
scilab.org/product/dic-mat-sci/M2SCI_doc.htm">end
plot(T_LLAMA,CONC)
xtitle("Variación De La Concentración Frente a la Tº","Tº de Llama","% De CH4")
 

[BrujoNic]

Cerrado.

Son tus deberes NO de los demás.