false_position_a.f90

	program false_position_a 

	implicit none
	real*8, parameter :: toleransi=1.d-5
	integer :: langkah
	real*8 :: a,b,x1,x2,kesrel
	
	
	do

 		write(*,*) 'Masukkan batas kiri, batas kanan:'
 		read (*,*) a,b

 		if (f(a)*f(b) < 0.d0) then
 			exit
		else
			write(*,*) 'Batas kiri-kanan tidak mengapit akar fungsi.'
			write(*,*) 'Coba lagi dengan nilai-nilai yang lain.'
		end if
	end do


	langkah=0
	
	do 
	
		x1=(a*f(b)-b*f(a))/(f(b)-f(a))

 		kesrel=abs((b-x1)/x1)

 		langkah=langkah+1
 		
 		if (kesrel < toleransi) exit

 		b=x1

	end do

	write(*,*) 'Pencarian akar konvergen pada langkah ke-', langkah
	write(*,*) 'Akar    = ',x1
	write(*,*) 'f(akar) = ',f(x1) 
	write(*,*) 'Kesalahan relatif =',kesrel

	stop

	contains

	function f(x) result(y)

	implicit none
	real*8, intent(in) :: x
	real*8 :: y

	y=cos(x)-x

	return

	end function f

	end program false_position_a