quad3d.f90

visual fortran常用数值算法集—原代码

SUBROUTINE quad3d(x1,x2,ss)
REAL ss,x1,x2,h
EXTERNAL h
!USES h,qgausx
call qgausx(h,x1,x2,ss)
END SUBROUTINE quad3d
FUNCTION f(zz)
REAL f,zz,func,x,y,z
COMMON /xyz/ x,y,z
!USES func
z=zz
f=func(x,y,z)
END FUNCTION f
FUNCTION g(yy)
REAL g,yy,f,z1,z2,x,y,z
EXTERNAL f
COMMON /xyz/ x,y,z
!USES f,qgausz,z1,z2
REAL ss
y=yy
call qgausz(f,z1(x,y),z2(x,y),ss)
g=ss
END FUNCTION g
FUNCTION h(xx)
REAL h,xx,g,y1,y2,x,y,z
EXTERNAL g
COMMON /xyz/ x,y,z
!USES g,qgausy,y1,y2
REAL ss
x=xx
call qgausy(g,y1(x),y2(x),ss)
h=ss
END FUNCTION h