Tim,
Thanks for your response.
X is my ODE solution vector for molar flow rates (Fi's), pressure (P) and
temperature (T). L is my integration range. I don't know if you had a chance
to check, but I have the code uploaded.
I will try your suggestion.
to evaluate the external function, would you
First of all, sorry if this is a silly question. I am getting practice in
Scilab
I have a set of 8 coupled ODEs and a set of 28 "external equations".
One example of external equation is my volumetric flow rate:
function v=F18(FE,FY,FH,FM,FP,FA,FL,FB,P,T)
v=(v0*ro0)/F17()
endfunction;
F17 is