8e5578ea |
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 |