Su-34

Люди добрые, помогите. Надо вот нелинейное уравнение решить методом ПИ в Фортране. Помогите с расположением кодов.

Эту функцию решил Методом половинного деления, верно?:
subroutine MPD(F,A,B,EPS,C,N,IER)
IF((F(A)*F(B)>0).OR.(A>B)) then
IER=1
return
end if
A1=A;B1=B;N=0;IER=0
do
Y=F(A1)
C=(A1+B1)/2
N=N+1
IF(B1-A1<EPS) exit
IF(Y*F(C)<0) then
B1=C
else
A1=C
end if
end do
end
program RGZ
external F
A=0;B=1
EPS=0.1
call MPD(F,A,B,EPS,C,N,IER)
print*,'x=',C,'IER=',IER
print*,'N=',N
end
function F(x)
f=3*x-cos(x)-1
end