Funcion spline

function f=ec(X)
epsi=X(1);
Ts=X(2);
%Temperatura en kelvin(K)
To=425+273.15;
%delta H en cal/mol
dH=-23490;
%presion en (atm)
P=1
%moles que Entran
ESO2=8;
EO2=11;
EN2=82.7619;nE=[8 11 82.7619 0];
%moles que Salen
SSO2=8-epsi;
SO2=11-0.5*epsi;
SN2=82.7619;
SSO3=epsi;

syms T
nS=[8-epsi, 11-0.5*epsi, 82.7619, epsi];
n=sum(nS);

%IntegralesCpO2=6.148+3.102*10^-3*T-0.923*10^-6*T^2;
CpN2=6.524+1.25*10^-3*T-0.001*10^-6*T^2;
CpSO2=7.116+9.512*10^-3*T+3.511*10^-6*T^2;
CpSO3=6.077+23.537*10^-3*T-0.687*10^-6*T^2;
CpT=[CpSO2 CpO2 CpN2 CpSO3];Cp=sum(nE.*INT(Cp,T,To,298));
F=subs(Cp+epsi*dH+sum(nS.*INT(CpT,T,298,Ts)));
EXP(-2.07*log(Ts)+3.1389*10^-3*Ts-0.3134*10^6*Ts^2+11467.1546/Ts+0.6589)-(SSO3/SSO2)*SQRT(n/(SO2*P));fprintf('TEMPERATURA A LA SALIDA =%f K\n',Ts)
fprintf('EPSILON =%f moles\n',epsi)
end



x =
6.1272 861.8324

TEMPERATURA A LA SALIDA = 870.000000 EPSILON=7.000000 molesTEMPERATURA A LA SALIDA = 870.000000 EPSILON=7.000000 moles TEMPERATURA A LA SALIDA =870.000000 EPSILON=7.000000 moles TEMPERATURA A LA SALIDA = 870.000000EPSILON=7.000000 moles



>> [X]= fsolve(@sistema,[7,870])
TEMPERATURA A LA SALIDA = 870.000000 EPSILON=7.000000 molesTEMPERATURA A LA SALIDA = 870.000000 EPSILON=7.000000 moles TEMPERATURA A LA SALIDA = 870.000000 EPSILON=7.000000 molesTEMPERATURA A LA SALIDA = 870.000000 EPSILON=7.000000 moles...