Tampilkan postingan dengan label Secant. Tampilkan semua postingan
Tampilkan postingan dengan label Secant. 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