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