|
Source Code:
|
|
Hasil Eksekusi:
|
|
Fungsi : log(x)-1
|
|
Nilai terendah : 1
|
|
Nilai tertinggi : 5
|
|
Batas toleransi : 0,001
|
|
|
|
Listing program
|
|
a=input('Enter
Function:','s');
|
|
f=inline(a);
|
|
xl=input('Enter lower
guess:') ;
|
|
xu=input('Enter upper
guess:');
|
|
tol=input('Enter
tolerance(recommended 0.001):');
|
|
for i=2:100
|
|
if f(xu)>f(xl)
|
|
xr=xu-((f(xu)*(xu-xl))/(f(xu)-f(xl)));
|
|
else
|
|
xr=xl-((f(xl)*(xl-xu))/(f(xl)-f(xu)));
|
|
end
|
|
if f(xu)*f(xr)<0
|
|
xl=xr;
|
|
else
|
|
xu=xr;
|
|
end
|
|
if f(xl)*f(xr)<0
|
|
xu=xr;
|
|
else
|
|
xl=xr;
|
|
end
|
|
xnew(1)=0;
|
|
xnew(i)=xr;
|
|
if
abs((xnew(i)-xnew(i-1))/xnew(i))<tol,break,end
|
|
end
|
|
str = ['Root: ',
num2str(xr), '']
|
|
Hasil eksekusi
|
|
Inputx0, xTol
|
|
iters = 1
|
|
dx = -f(0)/fDeriv(x0) (*
fcns f and fDeriv *)
|
|
root = x0 + dx
|
|
While (Abs(dx) > xTol)
|
|
dx =
-f(root)/fDeriv(root)
|
|
root = root + dx
|
|
iters = iters + 1
|
|
Endofwhile
|
|
Return root, iters
|
 Contoh
|
|
y=x^3 + 2*x^2 + 10*x –
20;
|
|
dalam selang [1, 1.5]
dengan epsilon = 10^-6= 0.000001:
|
