function main X=sort(randn(1,10)); Y=(X-1).^9; Approxroot = zeros(3,4); for i=1:4 Approxroot(1,i)=Newton(i,X,Y); Approxroot(2,i)=Bisect(i,X,Y); Approxroot(3,i)=Muller(i,X,Y); end Approxroot Error=abs(Approxroot-1)