DEFDBL A-Z SCREEN 9 COLOR 8, 7 WINDOW (0, -.5)-(3, 4) startomega = 0 endomega = 3 omegasteps = 2 ^ 9 deltaomega = (endomega - startomega) / omegasteps pi = 4 * ATN(1) twopi = 2 * pi DEF fnf (t) = 3 * SIN(7 * t - pi / 3) a = 0 b = 10 numberofsteps = 2 ^ 6 deltat = (b - a) / numberofsteps omega = startomega oldomega = omega oldpower = 0 FOR j = 1 TO omegasteps t = a + deltat / 2 accumS = 0 accumC = 0 power = 0 FOR k = 1 TO numberofsteps deltaS = 2 * fnf(t) * SIN(twopi * omega * t) * deltat / (b - a) accumS = accumS + deltaS deltaC = 2 * fnf(t) * COS(twopi * omega * t) * deltat / (b - a) accumC = accumC + deltaC t = t + deltat NEXT k power = SQR(accumS ^ 2 + accumC ^ 2) omega = omega + deltaomega LINE (oldomega, oldpower)-(omega, power) oldomega = omega oldpower = power NEXT j LINE (7 / twopi, 0)-(7 / twopi, 3), 6 'this marks the location of the 'peak in the power spectrum