c=4
f=(50./c)*(1-exp(-9*c/5))-10;
while abs(f)>1.0e-6
f=(50./c)*(1-exp(-9*c/5))-10;
df=(90./c)*exp(-9*c/5)-(50./c^2)*(1-exp(-9*c/5));
c=c-f/df
end
c=4
f=(50./c)*(1-exp(-9*c/5))-10;
while abs(f)>1.0e-6
f=(50./c)*(1-exp(-9*c/5))-10;
df=(90./c)*exp(-9*c/5)-(50./c^2)*(1-exp(-9*c/5));
c=c-f/df
end