10 DIM Y(32),A(15),B(15)
15 PRINT "INPUT N": INPUT N 
20 PRINT "INPUT VALUES OF FUNCTION (4/LINE)"  
25 FOR I=0 TO INT((N+2)/2)-1 
30 INPUT Y(4*I+1),Y(4*I+2),Y(4*I+3),Y(4*I+4): NEXT I 
35 PRINT "INPUT MAXIMUM HARMONIC TO BE CALCULATED (<=N)": INPUT M 
40 PRINT : PRINT "HARMONIC        COSINE COEFF.  SINE COEFF." 
45 C3=2/(2*N+1): K1=3.1415927*C3 
50 S1=SIN(K1): C1=COS(K1)
55 J,C=1: S=0: W=Y(1) 
60 U1,U2=0: I=2*N+1
65 U3=Y(I)+2*C*U1-U2: U2=U1: U1=U3 
70 I=I-1: IF I>1 THEN 65 
75 A(J)=C3*(W+C*U1-U2): IF J<>1 THEN 85 
80 A(1)=A(1)*.5
85 B(J)=C3*S*U1: PRINT J-1,A(J),B(J)  
87 IF J>=M+1 THEN 100 
90 Q=C1*C-S1*S: S=C1*S+S1*C: C=Q 
95 J=J+1: GOTO 60 
100 PRINT 
110 PRINT "INPUT X-VALUE TO BE EVALUATED. TO END PROGRAM INPUT 99999" 
120 INPUT W: IF W=99999 THEN 999 
130 I=0 
140 FOR J=1 TO M+1: I=I+A(J)*COS((J-1)*W)+B(J)*SIN((J-1)*W):NEXT J 
150 PRINT "Y=";I: PRINT : PRINT "INPUT X": GOTO 120  
999 END 
