Tampilkan postingan dengan label Regula Falsi. Tampilkan semua postingan
Tampilkan postingan dengan label Regula Falsi. Tampilkan semua postingan

Jumat, 14 Agustus 2015

Modul praktikum analisi numerik dengan menggunakan matlab Penyelesaian persamaan linier/tak linier



DAFTAR ISI

1.        Penyelesaian persamaan linier/tak linier 
A.    Metode Biseksi (Bagi Dua) ..................................................  1
B.     Metode Regula Falsi ............................................................  1
C.    Metode Newton Rhapson ....................................................  2
D.    Metode Secant .......................................................................  3
2.        Penyelesaian Sistem Persamaan Linier
A.    Jacobi .......................................................................................  3
3.        Interpolasi
A.  Lagrange .................................................................................  4
4.        Integral
A.  Integral Trapesium ...............................................................  5
B.   Integral Simson  ...............................................................  5
C.  Integral Simson  ...............................................................  6
5.        Persamaan Differensial Biasa
A.    Meode Euler ..........................................................................  6
B.     Metode Runge Kutta Orde IV ............................................  6
6.        Lampiran (Dasar-dasar Matlab)............................................... 8
 


MODUL PRAKTIKUM METODE NUMERIK
1.    Penyelesaian Persamaan Linier/Tak Linier
A.    Metode Biseksi (Bagi Dua)
clc, clear
e1=1/100000000000;
x1=-4;x2=3;xt=(x1+x2)/2;
iterasi=100;
y=inline('x^3','x')
ezplot(y)
axis ([-5 5 -5 5])
grid on
for i=1:iterasi
    if abs(y(x1))<e1
        Pada_iterasi_ke=i;Pada_iterasi_ke, nilaiX=x1; nilaiX,         break
    elseif abs(y(x2))<e1
        Pada_iterasi_ke=i;Pada_iterasi_ke, nilaiX=x2; nilaiX,         break
    elseif abs(y(xt))<e1
        Pada_iterasi_ke=i;Pada_iterasi_ke, nilaiX=xt; nilaiX,         break
    elseif y(x1)*y(x2)<0
        xt=(x1+x2)/2;
        if y(x1)*y(xt)<0
            x1=x1;            x2=xt;
        elseif y(x2)*y(xt)<0
            x1=xt;            x2=x2;
        end
    else
        disp 'tidak ada akar dalam rentang ini',         break
    end
end
B.     Metode Regula Falsi
clc, clear
e1=0.001;

x1=-2;x2=3;
iterasi=100;

y=inline('x^3','x')
ezplot(y)
axis ([-5 5 -5 5]);grid on

for i=1:iterasi
    xc=x2-((y(x2)*(x2-x1))/(y(x2)-y(x1)))
    if abs(y(x1))<e1
        Pada_iterasi_ke=i;Pada_iterasi_ke, nilaiX=x1; nilaiX,         break
    elseif abs(y(x2))<e1
        Pada_iterasi_ke=i;Pada_iterasi_ke, nilaiX=x2; nilaiX,         break
    elseif abs(y(xc))<e1
        Pada_iterasi_ke=i;Pada_iterasi_ke, nilaiX=xc; nilaiX,         break
    elseif y(x1)*y(x2)<0
        if y(x1)*y(xc)<0
            x1=x1;            x2=xc;
        elseif y(x2)*y(xc)<0
            x1=xc;            x2=x2;
        end
    else
        disp 'tidak ada akar dalam rentang ini'
        break
    end
end
C.     Metode Newton Rhapson
clc, clear
e1=0.00000000000000001;
x=3;
iterasi=100
y=inline('x^3','x')
dy=diff('x^3','x');
dy=inline(dy)
ezplot(y)
axis ([-5 5 -5 5])
grid on

for i=1:iterasi
    x=x-y(x)/dy(x)
    if abs(y(x))<e1
        Pada_iterasi_ke=i;Pada_iterasi_ke, nilaiX=x; nilaiX,         break
    end
end
D.    Metode Secant
clc, clear
e1=0.001;
x1=3; x2=2;
iterasi=100
y=inline('x^3','x')
ezplot(y)
axis ([-5 5 -5 5]); grid on

for i=1:iterasi
    x3=x2-(y(x2)*((x2-x1)/(y(x2)-y(x1))))
    if abs(x3-x2)<e1
        Pada_iterasi_ke=i;Pada_iterasi_ke, nilaiX=x3; nilaiX,        break
    end
    x1=x2;    x2=x3;
end