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