Mathematical Physics - Volume II - Numerical Methods

Chapter 3. Comparison of finite element method and finite difference method

92

b = (155.0-sq15)/2400.0 wt(1) = 9.0/80.0

wt(2) = a wt(3) = a wt(4) = a wt(5) = b wt(6) = b wt(7) = b

a = (6.0+sq15)/21.0 b = (6.0-sq15)/21.0

c = (9.0+2.0*sq15)/21.0 d = (9.0-2.0*sq15)/21.0

xit(1,1) = 1.0/3.0 xit(2,1) = 1.0/3.0

xit(1,2) = a xit(2,2) = d xit(1,3) = d xit(2,3) = a xit(1,4) = a xit(2,4) = a xit(1,5) = b xit(2,5) = c xit(1,6) = c xit(2,6) = b xit(1,7) = b xit(2,7) = b c tri point rule for line segment [-1,+1] etaq(1) = -sq35 etaq(2) = 0.

etaq(3) = sq35 ww(1) = 5./9. ww(2) = 8./9. ww(3) = 5./9. RETURN END

SUBROUTINE prep CHARACTER title*80 COMMON /ctitl / title READ (5,*) READ (5,99001) title IF ( title.EQ.’end’ .OR. title.EQ.’END’ ) THEN

CLOSE (6) CLOSE (5) STOP

ENDIF PRINT * PRINT * , ’reading problem data ...’ WRITE (6,*)