So, first off, the answer to this problem should be roughly .0028. But, somewhere in the code there is an issue and it will not compute correctly. Can you help?
clcclearv = input('Enter initial guess for the volume of methane: '); R = .518;pc = 4600; tc = 191; T = -40 + 273; P = 65000; a = .427*R^2*tc^2.5/pc; b = .0866*R*tc/pc; gv =@(v)(R*T)-(((a*(v-b)^2)/(v*(v+b)*sqrt(T))) + P * b)/P;n = 4;es = .5*10^(2-n);[r, I] = Redlich_Kwong(es, v, gv); fprintf('\nThe calculated volume is : %.6f',r); fprintf('\nIt took %1.0d iterations to converge.',I);
Function
function [r, I] = Redlich_Kwong(es, v, gv)I = 0; ea = 1;while ea > es r = gv(v); ea = 100*abs((r-v)/r); v = r; I = I +1; end end
Best Answer