파이(원주율)을 구하는중에 문제가..ㅠㅠ[fortran]

program pi
implicit double precision(a-h,o-z)
implicit integer*4(i-n)
double precision a(2)
common/seed/iseed

open(unit=10,file='pi.txt')
iseed=irtc()

do i=0, 10000, 500

carea=0.

do j=1, i
a(1)=ranf()*10-5
a(2)=ranf()*10-5

write(*,*) a(1),a(2)
len=length(a(1),a(2))
write(*,*) len
if(len.lt.5) carea=carea+1
enddo

write(10,*) carea
write(10,*) 'pi=', carea*4/i
enddo

stop
end

real function length(a,b)
double precision a,b

length=sqrt(a**2+b**2)

return
end

이렇게 만들어 봤는데요..function을 사용하니까 갑자기 값이 나오지를 않네요..function 사용하지 않고 식을 main 에 바로 쓰면 잘 되는데...왜 이런거죠??