10 REM Fourier- 20 REM Analyse 50 INPUT "ZAHL DER HARMONISCHEN=";NH 60 INPUT "ZAHL DER PUNKTE=";NP 70 IF NP<2*NH+1 THEN 50 80 RETURN 100 REM BEISPIEL SAEGEZAHN 1 105 GOSUB 50 110 DIM F(NP) 115 F(0)=-1:DY=2/NP 120 FOR I=1 TO NP-1 125 F(I)=F(I-1)+DY:NEXTI 128 F(0)=0 130 GOSUB 500 140 STOP 200 REM BEISPIEL SAEGEZAHN 2 205 GOSUB 50 210 DIMF(NP) 215 F(0)=0:DY=2/NP 220 FOR I=1 TO NP/2-1 225 F(I)=F(I-1)+DY:NEXTI 230 F(NP/2)=0.5 240 FOR I=NP/2+1 TO NP-1 245 F(I)=0:NEXTI 250 GOSUB 500 260 STOP 300 REM BEISPIEL RECHTECK 302 GOSUB 50 305 DIMF(NP) 310 F(0)=0:F(NP/2)=0 315 FOR I=1 TO NP/2-1:F(I)=1:NEXTI 320 FORI=NP/2+1 TO NP-1:F(I)=-1:NEXTI 330 GOSUB 500 340 STOP 500 PI=3.141592653 502 PRINT 503 PRINT" N";TAB(7);"COS";TAB(13);"SIN";:PRINT:PRINT 505 PRINT 510 Z=2/NP 520 C1=COS(Z*PI):S1=SIN(Z*PI) 530 C=1:S=0:P=0 540 U1=0:U2=0 550 KLN=NP-1 560 U0=F(KLN)+2*C*U1-U2 570 U2=U1:U1=U0:KLN=KLN-1 580 IF KLN=0 GOTO 600 590 GOTO 560 600 A=Z*(F(0)+C*U1-U2) 605 B=Z*S*U1 610 EPS=1E-6 620 IF ABS(A)