5 * Revision 1.1.1.1 1996/04/01 15:02:25 mclareni
10 FUNCTION TRISUM(X,A,N,B,M,IOP)
12 C PURPOSE = TO COMPUTE THE SUM OF THE TRIGONOMETRIC SERIES
14 C TRISUM = A(1)/2+A(2)*COS(X)+....................A(N)*COS((N-1)*X)
15 C +B(1)*SIN(X)+....................B(M)*SIN(M*X)
19 C X = FUNCTION ARGUMENT
20 C A = ARRAY DIMENSIONED A(N),KEEPING THE COEFFICIENTS FOR THE
22 C N = THE NUMBER OF COEFFICIENTS IN THE COSINE SERIE
23 C B = ARRAY DIMENSIONED B(M),KEEPING THE COEFFICIENTS FOR THE
25 C M = THE NUMBER OF COEFFICIENTS IN THE SINE SERIE
26 C IOP = OPTION NUMBER , IOP = 1 , THE GENERAL CASE
27 C IOP = 2 , THE EVEN CASE
28 C IOP = 3 , THE ODD CASE
32 C*UL 1000 COF=2.*COS(X)
36 IF (IOP.EQ.3) GO TO 1030
38 C COMPUTE SUM OF COSINE TERMS
44 IF (N.EQ.1) GO TO 1020
54 1020 TRISUM=.5*(A0-A2)
56 1030 IF (IOP.EQ.2) RETURN
58 C COMPUTE SUM OF SINUS TERMS
64 IF (M.EQ.1) GO TO 1050
74 1050 TRISUM=TRISUM+A0*SIN(X)