NEWXYZ
 

In the subroutine NEWXYZ, I used four functions to calculate Bessel function. But when I write the result from the function and the subroutine, they are different. I just give the value of the function to a variable in the subroutine. They should give me the value that the function calculated. But why not?