fortran程序,subroutine solvereal::i,x,y,mreal::h(i),d(i),a(i),