下面程序中为什么输出的j（a）与sin（a）/a不同，怎么该使结果相同？

program main
external j
real a
a=1.
print*,j(a),sin(a)/a
end
real function j(r)
real r
j=sin(r)/r
end function j