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

