метод ньютона и метод простых итераций

[Igor1997]

здравствуйте!
program uravnenie
implicit none
real :: F
real:: eps1, eps2, m, R
integer:: k
real:: b, a
print*, "выберите способ решения"
print*, "1 - метод Ньютона"
print*, "2 - метод простых итераций"
print*, "3 - метод дихотомии"
read(*,*) k
print*,"введите левую границу"
read(*,*) a
print*, "введите правую границу"
read(*,*) b
m=-8.0
select case(k)
case(1)
print*, "выберите точность eps1"
read(*,*) eps1
call newton(b, eps1, m, r)

print*, "корень=", R
print*, "F(r)=", f(r)


case(2)

print*, "выберите точность eps1"
read(*,*) eps1

call iteracii(b, eps1, m, r)
print*, "корень=", R
print*, "F(r)=", f(r)




case(3)
print*,"выберите точность eps1"
read(*,*) eps1
print*,"выберите точность eps2"
read(*,*) eps2

call dihotomy(a, b, eps1, eps2, m, R)
print*, "корень=", R
print*, "F(r)=", f(r)





case default
print*,"Ошибка при выборе метода"
go to 1
end select



open(unit=1, file='C:\Users\data.out')
write(1,*) R, F(r)
stop

end program
real function f(x)
real :: x
f=x**3-3*x**2+5 !моя функция
end function f
real function df(x)
real:: x
df=((f(x+epsilon(x))-f(x))/epsilon(x))
end function df
subroutine newton (b, eps1, m, r)






return
end subroutine newton
subroutine iteracii (b, eps1, m, r)


subroutine dihotomy(a, b, eps1, eps2, m, c)
real:: b, a, eps1, eps2, m, c
integer:: n
n=0
do while (abs(a-b)>eps.or.N<1000)
c=(a+b)/2
if (f(a)*f(c)<0)then
b=c
else
a=c
end if
N=N+1
end do
c=(a+b)/2

end subroutine dihotomy
как организовать метод ньютона и простых итераций?