%——————————————
% Program Latihan
%——————————————
clear all;
clc;

fname1=input(‘nama file output ‘);
foutput=fopen(fname1,’w');
fprintf(foutput,’%10s %10s %12s\n’,'iterasi’,'xx’,'yy’);
fprintf(foutput,’——————————\n’);

v=[-3:0.05:3];
syms x;
f=(x^4-x+1);
c=diff(f);
z=diff(f,2);
d=subs(f,v);
plot(v,d);
y=3;
p=1;
n=1;
while p>0;
cc=subs(c,y);
zz=subs(z,y);
hold on;
k=subs(f,y);
plot(y,k);
ww=cc/zz;
w=y-ww;
fprintf(foutput,’%5f %10f %13f\n’,n,y,p);
p=y-w;
q=y;
y=w;
end;
fclose(foutput);

gw puas ama nich coding, lumayan buat pemula.