!copyright (C) 2001 MSC-RPN COMM %%%RPNPHY%%% ***FUNCTION TRAPEZ2 *REAL FUNCTION TRAPEZ2(DEL,F,N,NM) 5 #include "impnone.cdk"
INTEGER N, NM, I REAL F(N),DEL(NM) * *Author * L.Garand (1989) * *Revision * 001 G.Pellerin(Mar90)Standard documentation * *Object * to calculate integrals * *Arguments * * - Input - * DEL thickness of layers * F variable to be integrated * N dimension of F * NM dimension of DEL * ** * TRAPEZ2=0. DO 10 I=1,NM 10 TRAPEZ2=TRAPEZ2+(F(I)+F(I+1))/2.*DEL(I) RETURN END