T=10u pi=3.14159 f=500 freq=f Tsample=50m set colwidth=10 printtext -n "freq" printtext -n "c2" printtext -n "s" printtext -n "corgain2" printtext -n "corgain2" printtext -n "numcycles" printtext -n "unity" printtext while freq < 25k period = 1/freq numcycles= ceil(Tsample/period) f=numcycles/Tsample w=2*pi*f s=sin(w*T*180/pi) c=cos(w*T*180/pi) gain = -2*pi*f*T cor=tan(f*pi*T*180/pi)/(pi*T*f) corgain2=ceil(cor*gain/2*2^14) c2=ceil(2*c*2^14) s =ceil(s*2^14) unity = (s*s+c2*c2/4)*2^-14 set colwidth=1 printtext -n ";" set colwidth=10 printval f set colwidth=1 printtext -n "," set colwidth=10 printval c2 set colwidth=1 printtext -n "," set colwidth=10 printval s set colwidth=1 printtext -n "," set colwidth=10 printval corgain2 printval numcycles printval unity printtext printtext -n ".int" printval c2 printtext -n "," printval s printtext -n "," printval corgain2 printtext freq=1.1*freq end