main.m
Source code :
clc;
disp ('TUGAS 4');
disp ('METODE NEWTON RAPHSON');
disp (' ');
disp ('Nama : Desianny Amalia Fauziah');
disp ('NIM : 5302413085');
disp ('Rombel : 02');
disp (' ');
disp('PERSAMAAN : exp(x)-2-(x.^2)');
disp(' ');
x0=input('Masukkan nilai x awal = ');
f_x = inline('exp(x)-2-(x.^2)','x'); %definisi fungsi
newtonrapshon(f_x, 0.5, 0.01);
newtonraphson.m
Source code :
function newtonrapshon(fx,x0,e)
disp(' ');
i=0;
M = 9; %error awal untuk syarat perulangan
xb = 0;
disp('===============================================');
disp(' Iterasi xi f(xi) e');
disp('===============================================');
while(M>e)
syms x;
gx = diff(sym(fx));
fgx=subs(gx,x,x0);
xb = x0 - (fx(x0)/fgx);
M = abs(x0-xb); %error pendekatan
x0 = xb;
i=i+1;
fprintf('%6.0f %12.6f %12.6f %12.6f\n',i,xb,fx(x0),M);
error(i) = abs(M);
end
fprintf('Akarnya adalah %10.8f\n',xb);
plot(error)
end
Tidak ada komentar:
Posting Komentar