Introduction - If you have any usage issues, please Google them yourself
The following is compiled with the Fortran 90 language, the Fortran PowerStation 4.0 integrated development environment running under the ICCG method computer program, the program is also suitable for use in Visual Fortran 6.5 environment, and does not require changes. See 5.2.2,5.2.3 the algorithm.