1        program test
2        double precision a,b,c
3        data a,b/1.0d-46,1.0d0/
4        c=fun(a,b)
5        print*,'in main: fun=',c
6        end
7        double precision function fun(a,b)
8        double precision a,b
9        print*,'in sub: a,b=',a,b
10        fun=a*b
11        print*,'in sub: fun=',fun
12        return
13        end
14