1 INTEGER FUNCTION mltmod(a,s,m)
31 INTEGER a0,a1,k,p,q,qh,rh
38 IF (.NOT. (a.LE.0.OR.a.GE.m.OR.s.LE.0.OR.s.GE.m)) go to 10
39 WRITE (*,*)
' A, M, S out of order in MLTMOD - ABORT!'
40 WRITE (*,*)
' A = ',a,
' S = ',s,
' M = ',m
41 WRITE (*,*)
' MLTMOD requires: 0 < A < M; 0 < S < M'
42 CALL xstopx(
' A, M, S out of order in MLTMOD - ABORT!')
44 10
IF (.NOT. (a.LT.h)) go to 20
53 IF (.NOT. (a1.GE.h)) go to 50
56 p = h* (s-k*qh) - k*rh
57 30
IF (.NOT. (p.LT.0)) go to 40
67 60
IF (.NOT. (a1.NE.0)) go to 90
73 70
IF (.NOT. (p.LT.0)) go to 80
82 p = h* (p-k*qh) - k*rh
83 100
IF (.NOT. (p.LT.0)) go to 110
88 120
IF (.NOT. (a0.NE.0)) go to 150
97 130
IF (.NOT. (p.LT.0)) go to 140