You need to rewrite your code so that averything ‘myfun’ needs to calculate that are functions of ‘x’ are calculated within it, using the ‘x’ vector supplied to it by fsolve.
The statements must be in the correct order, so that the values needed later in the code are calculated first.
For example:
pC=yC*P;
appears under ‘%outlet gas partial pressure’ however
yC=NC/(NC+NH+G2-NW+NCH);
is calculated later in ‘%outlet gas composition’
I am certain you can get your code to run successfully. However you need to write ‘myfun’ so that it calculates its variables in the correct order. Then use fsolve to optimise ‘myfun’, returning the optimal values for ‘x’.
