header <napatool.hdr>

title  "Intermodulation, a theoretical study"

fs     1.0e9

dvar   k3   -0.01       // in physical circuit, signal is clipped, k3 must therefore be negative
dvar   k2    0.1
dvar   k1    1.0

dvar   f1   10.0e6
dvar   f2   f1 + 1.0e6


ivar   npt   POWEROF2(18)
ivar   nstp  51

node   a0   duser    synchro_linsweep   -25.0 25.0  nstp  (up)
node   a1   dalgebra DB2LIN(a0, 1.0)

export a0

#*  DUT

node  in1   cosine   0.0  a1 f1 0.0
node  in2   cosine   0.0  a1 f2 0.0
node  in    sum   in1  in2

node   s0   poly  k3 k2 k1 0.0  in


tool  im3  "~/Output/im3a.out"  s0 1.0  f1 f2  npt

terminate  nstp <= TOOL_INDEX