Matlab code used:

% Given:

% dX/dL = eta*RNH3*A/2*FN2

% Known parameters

eta = 0.6527

RNH3 = 56791.1128

A = 0

FN2 = 135.43

% Function used

func = @(A) eta*RNH3*A/2*FN2

X = func(A)

while abs(x-y)>0.00001

x=y

n=n+1

y=func(A)

end

I can't seem to iterate this to 97% conversion to determine the area of reactor (A) required using Matlab. Pls correct the Matlab code too.