Version 4
SHEET 1 1420 1180
WIRE 96 240 80 240
WIRE 208 240 192 240
WIRE 80 256 80 240
WIRE 192 256 192 240
WIRE 544 272 544 256
WIRE 544 384 544 352
WIRE 608 384 544 384
WIRE 544 416 544 384
WIRE 320 448 304 448
WIRE 304 464 304 448
WIRE 208 496 192 496
WIRE 192 512 192 496
WIRE 544 528 544 496
WIRE 96 544 80 544
WIRE 80 560 80 544
FLAG 80 336 0
FLAG 80 240 t
FLAG 192 336 0
FLAG 192 240 u
FLAG 608 384 1
FLAG 80 640 0
FLAG 80 544 3
FLAG 544 528 0
FLAG 544 256 0
FLAG 192 592 0
FLAG 192 496 4
FLAG 304 544 0
FLAG 304 448 5
SYMBOL res 528 400 R0
WINDOW 3 36 68 Left 0
SYMATTR Value R=V(u)
SYMATTR InstName R1
SYMBOL bv 192 240 R0
WINDOW 3 -33 136 Left 0
SYMATTR Value V=if(V(t) < 0, Rn(V(t)), Rp(V(t)))
SYMATTR InstName B1
SYMBOL voltage 80 240 R0
SYMATTR InstName V1
SYMATTR Value 0
SYMBOL bv 80 544 R0
WINDOW 3 -33 136 Left 0
SYMATTR Value V=(sqrt((A**2 - 4*B)*R0**2 + 4*B*R0*V(1)/1m) - A*R0)/(2*B*R0)
SYMATTR InstName B2
SYMBOL current 544 272 R0
SYMATTR InstName I1
SYMATTR Value 1m
SYMBOL bv 192 496 R0
WINDOW 3 -33 136 Left 0
SYMATTR Value V=if(V(3) < 0, fv4(V(3)), 0)
SYMATTR InstName B3
SYMBOL bv 304 448 R0
WINDOW 3 -33 136 Left 0
SYMATTR Value V=V(3) - V(4)
SYMATTR InstName B4
TEXT 48 720 Left 0 !.option numdgt=16 measdgt=8\n.param R0=100\n.param A=3.9083e-3 B=-5.775e-7 C=-4.183e-12\n.func Rn(t) { R0*(1 + A*t + B*t**2 + C*(t-100)*t**3) }\n.func Rp(t) { R0*(1 + A*t + B*t**2) }\n.dc V1 -100 100 0.01
TEXT 48 976 Left 0 ;.func fv4(x) { c4*x**4 + c3*x**3 + c2*x**2 + c1*x }\n.param c4 = -9.0548722256116693E-10 ,c3 = 1.2002237069770213E-7\n+         c2 = 4.6535300827000332E-7 c1 = 5.0191938484069306E-6
TEXT 48 904 Left 0 !.func fv4(x) { c2*x**2 + c1*x }\n.param c2 = -3.3462046194357326E-5 c1 = -0.0012784081351104
RECTANGLE Normal 656 512 512 400 2