Sabtu, 24 Desember 2011

program Matlab Newton

%program newton dengan matlab simbolik
%created by masduki
%input program
f=input ('masukkan fungsi (dalam x)= ');
x0=input('masukkan nilai awal (x0) = ');
n=input ('masukkan banyak interasi = ');
exact=input ('masukkan salah satu solusi eksak = ');

%membuat header table newton
disp (' ')
disp ('           METODE NEWTON             ')
disp ('=====================================')
disp (' i        akar            error             ')
disp ('=====================================')

%perulangan pada metode newton
for i= 1:n
   fx0=subs(f,'x',x0);                %menghitung nilai f(xi)
   fderx0=subs(diff(f),'x',x0);    %menghitung nilai f'(xi)
   xr=x0-fx0/fderx0;                    %menghitung akar dengan newton
   eror=abs(exact-xr);                %menghitung error absolut
   fprintf ('%2.0f    %5.5f        %2.5f\n',i,xr,eror)        %mencetak bentuk tabel
   x0=xr ;                                %mengganti harga xi dengan xr
end

%membuat footer table newton
disp ('=============================')



output dari program diatas sebagai berikut :

Masukkan fungsi (dalam x)         =......................
Masukkan nilai awal (x0)              =...............................
Masukkan banyak iterasi               = ..............................
Masukkan salah satu solusi eksak =....................................

METODE NEWTON
================================================
                                    i                                   akar                                  error
================================================

program Matlab fixed point

%program fixed point dengan matlab simbolik
%input program
g=input ('masukkan fungsi (dalam x)= ');
x0=input('masukkan nilai awal (x0) = ');
n=input ('masukkan banyak interasi = ');
exact=input ('masukkan salah satu solusi eksak = ');

%membuat header table fixed point
disp (' ')
disp ('           METODE FIXED POINT        ')
disp ('=====================================')
disp (' i        akar             error            ')
disp ('=====================================')

%perulangan pada metode fixed point
for i= 1:n
   gx0=subs(g,'x',x0);                %menghitung nilai g(x0)
   xr=gx0;                                %menghitung akar dengan fixed
   eror=abs(exact-xr);                %menghitung error absolut
   fprintf ('%2.0f    %5.5f        %2.5f\n',i,xr,eror)        %mencetak bentuk tabel
   x0=xr ;                                %mengganti nilai x0 dengan xr
end

%membuat footer table fixed
disp ('=============================')

output dari program diatas sebagai berikut
Masukkan fungsi g(dalam x)         =......................
Masukkan nilai awal (x0)              =...............................
Masukkan banyak iterasi               = ..............................
Masukkan salah satu solusi eksak =....................................
METODE FIXED POINT
================================================
                                    i                                   akar                                  error
================================================