子例程仅适用于某些值
PROGRAM po
IMPLICIT NONE
INTEGER :: i, j, nc, nd,ok,iter
REAL :: alph, bet, chi
REAL, DIMENSION(:), ALLOCATABLE :: uexact, x,u,up2, xd , uexactd
REAL ...