*
simulator lang = spectre

ahdl_include "varistor.va"
x1 (top 0) varistor T=1 C=700pF L=13.0nH B1=2.7233755 B2=0.0258453 
+  B3=-0.0005746 B4=0.0046033

i1 (0 top) isource type=sine dc=1e-5 sinedc=1e3 freq=1e3 ampl=1e7


* to verify the linear interpolation
*dc0 dc dev=i1 start=-0.2e-6 stop=0.2e-6


* these three sweeps replicate Figure 2: V-I characteristic
* with 10% tolerance bands
dc1 dc dev=i1 start=1e-5 stop=1e4 dec=10

alter1 alter dev=x1 param=T value=0.9
dc2 dc dev=i1 start=1e-5 stop=1e4 dec=10

alter2 alter dev=x1 param=T value=1.1
dc3 dc dev=i1 start=1e-5 stop=1e4 dec=10


* can you replicate the transient pulse of Figure 3?