{ z = 0.0; p = 1.00*(2*3.14159); pwrong = p*1.02; pbad = p*1.43; m0 = 0.0; a = 3.0; e = 0.00; dt = 0.29; for (t = 0.0; t < 100; t += dt) { phase = (t - z)/p; iphase = int(phase); phase = phase - iphase; wrongphase = (t - z)/pwrong; iwrongphase = int(wrongphase); wrongphase = wrongphase - iwrongphase; badphase = (t - z)/pbad; ibadphase = int(badphase); badphase = badphase - ibadphase; arg = phase*(2*3.14159); mag = m0 + a*sin(t) + e*rand(); printf "%9.4f %6.3f %9.3f %6.3f %6.3f \n", t, phase, mag, wrongphase, badphase; } }