clear all; clc; x1=input('x1 = '); e=input('e = '); g1=-3/x1^2; if(abs(g1)<1) x2=1+3/x1; while(abs(x2-x1)>e) x1=x2; x2=1+3/x1 end else disp('is failed') end