Monday, May 25, 2009

Listing Program Newton Raphson

Berikut ini adalah contoh listing program yang dapat digunakan untuk mencari akar persamaan dengan menggunakan program Matlab, dan metode Newton Raphson
clear;
clc;
err=0.0005
x=input('nilai X0 = ');
a=eps;
eps=1;
no=0;
clc;


fprintf('mencari akar persamaan dengn metode newton rhapson\n');
fprintf('persamaan : 2*X^3 + 4*x + 2 = 0 \n');
fprintf('taksiran awal : %5.0f\n',x);
fprintf('========================================================================\n');
fprintf(' iterasi f f` x1 selisih \n');
fprintf('========================================================================\n');

while eps>err
no=no+1;
f1=2*x^3+4*x+2;
f2=6*x^2+4;
x1=x-(f1/f2);
eps=abs(x1-x);
x=x1;
fprintf(' %3d %8.5f %8.5f %8.5f %8.5f\n', no, f1, f2, x1, eps);
end
fprintf('========================================================================\n');
fprintf('pada iterasi ke %d, selisih < %5.3f.\n',no,err);
fprintf('jadi persamaannya adalah %7.5f.\n',x1);




Hasil Program :
err =

5.0000e-004

nilai X0 = 1

mencari akar persamaan dengn metode newton rhapson
persamaan : 2*X^3 + 4*x + 2 = 0
taksiran awal : 1
========================================================================
iterasi f f` x1 selisih
========================================================================
1 8.00000 10.00000 0.20000 0.80000
2 2.81600 4.24000 -0.46415 0.66415
3 -0.05659 5.29262 -0.45346 0.01069
4 -0.00032 5.23375 -0.45340 0.00006
========================================================================

pada iterasi ke 4, selisih < 0.001.
jadi persamaannya adalah -0.45340

No comments:

Post a Comment