* SPICE netlist
C1 N006 0 10n
C2 N003 N006 10n
R1 N002 N006 27k
R2 N002 N003 33k
V1 x0 0 AC 1
Rkot x0 N002 1
C3 N007 0 3n
C4 N004 0 7n
R4 N001 N007 27k
R5 N001 N004 33k
Rerror x0 N001 1
B1 N008 N009 V=V(x0)*zb2
B2 N005 N008 V=V(x1)*zb1
B3 y0 N005 V=V(x2)*zb0
B4 N010 N009 V=V(y1)*za1
B5 0 N010 V=V(y2)*za0
B8 y1 0 V=V(y0) Laplace=exp(-s*ts)
B6 y2 0 V=V(y1) Laplace=exp(-s*ts)
B7 x1 0 V=V(x0) Laplace=exp(-s*ts)
B9 x2 0 V=V(x1) Laplace=exp(-s*ts)
B10 ys 0 V=V(x0) Laplace=ga*(s**2 + B1*s)/(s**2 + A1*s + A0)
.ac oct 100 10 100k
.param fs 44100
.param ts {1/fs}
;.param kz {2*pi*10k/tan(2*pi*10k*ts/2)}
.param kz {2*fs}
.param kz2 {kz*kz}
.param r1 27k
.param r2 33k
.param c1 10n
.param c2 10n
.param ga {(r1+r2)/(r1*r2)}
.param gs {1/ga}
.param z0 {1/(r2*c2)}
.param p0 {1/((r1+r2)*c2)}
.param B2 0
.param B1 {ga*p0}
.param B0 ga
.param A0 1
.param A1 {(c1*gs*z0+1)/(c1*gs)}
.param A2 {p0/(c1*gs)}
.param zb0 {(B0*kz2 + B1*kz + B2)/(A0*kz2+A1*kz+A2)}
.param zb1 {(2*B2 - 2*B0*kz2)/(A0*kz2+A1*kz+A2)}
.param zb2 {(B0*kz2 - B1*kz + B2)/(A0*kz2+A1*kz+A2)}
.param za1 {(2*A2-2*A0*kz2)/(A0*kz2+A1*kz+A2)}
.param za0 {(A0*kz2-A1*kz+A2)/(A0*kz2+A1*kz+A2)}

* LTspice schematic (.asc)
Version 4
SHEET 1 2144 2016
WIRE 128 -1008 0 -1008
WIRE 0 -896 0 -1008
WIRE 0 -704 0 -816
WIRE -608 16 -1072 16
WIRE -464 16 -528 16
WIRE -384 16 -464 16
WIRE -384 112 -384 16
WIRE 288 128 192 128
WIRE -1072 144 -1072 16
WIRE -976 144 -1072 144
WIRE -832 144 -896 144
WIRE -752 144 -832 144
WIRE -752 192 -752 144
WIRE 192 208 192 128
WIRE -832 240 -832 144
WIRE -464 240 -464 16
WIRE -752 320 -752 272
WIRE -384 320 -384 192
WIRE -1072 336 -1072 144
WIRE 192 336 192 288
WIRE 688 352 640 352
WIRE 1072 352 1024 352
WIRE -832 384 -832 320
WIRE -752 384 -832 384
WIRE -832 432 -832 384
WIRE -464 432 -464 320
WIRE 640 432 640 352
WIRE 1024 432 1024 352
WIRE 192 464 192 416
WIRE -1072 528 -1072 416
WIRE -832 528 -832 496
WIRE -832 528 -1072 528
WIRE -464 528 -464 496
WIRE -384 528 -384 384
WIRE -384 528 -464 528
WIRE -832 560 -832 528
WIRE -464 560 -464 528
WIRE 640 560 640 512
WIRE 1024 560 1024 512
WIRE 192 688 192 544
WIRE 688 736 640 736
WIRE 1072 736 1024 736
WIRE 192 816 192 768
WIRE 640 816 640 736
WIRE 1024 816 1024 736
WIRE 640 944 640 896
WIRE 1024 944 1024 896
WIRE 192 1008 192 896
FLAG -832 560 0
FLAG -464 560 0
FLAG 288 128 y0
FLAG 640 944 0
FLAG 688 736 y1
FLAG 1024 944 0
FLAG 1072 736 y2
FLAG 640 560 0
FLAG 688 352 x1
FLAG 1024 560 0
FLAG 1072 352 x2
FLAG -1072 16 x0
FLAG 192 1008 0
FLAG 0 -704 0
FLAG 128 -1008 ys
SYMBOL cap -848 432 R0
SYMATTR InstName C1
SYMATTR Value 10n
SYMBOL cap -768 320 R0
SYMATTR InstName C2
SYMATTR Value 10n
SYMBOL res -848 224 R0
SYMATTR InstName R1
SYMATTR Value 27k
SYMBOL res -768 176 R0
SYMATTR InstName R2
SYMATTR Value 33k
SYMBOL voltage -1072 320 R0
SYMATTR InstName V1
SYMATTR Value ""
SYMATTR Value2 AC 1
SYMBOL res -992 160 R270
SYMATTR InstName Rkot
SYMATTR Value 1u
SYMBOL cap -480 432 R0
SYMATTR InstName C3
SYMATTR Value 3n
SYMBOL cap -400 320 R0
SYMATTR InstName C4
SYMATTR Value 7n
SYMBOL res -480 224 R0
SYMATTR InstName R4
SYMATTR Value 27k
SYMBOL res -400 96 R0
SYMATTR InstName R5
SYMATTR Value 33k
SYMBOL res -624 32 R270
SYMATTR InstName Rerror
SYMATTR Value 1u
SYMBOL bv 192 448 R0
SYMATTR InstName B1
SYMATTR Value V=V(x0)*zb2
SYMBOL bv 192 320 R0
SYMATTR InstName B2
SYMATTR Value V=V(x1)*zb1
SYMBOL bv 192 192 R0
SYMATTR InstName B3
SYMATTR Value V=V(x2)*zb0
SYMBOL bv 192 784 R180
SYMATTR InstName B4
SYMATTR Value V=V(y1)*za1
SYMBOL bv 192 912 R180
SYMATTR InstName B5
SYMATTR Value V=V(y2)*za0
SYMBOL bv 640 800 R0
SYMATTR InstName B8
SYMATTR Value V=V(y0) Laplace=exp(-s*ts)
SYMBOL bv 1024 800 R0
SYMATTR InstName B6
SYMATTR Value V=V(y1) Laplace=exp(-s*ts)
SYMBOL bv 640 416 R0
SYMATTR InstName B7
SYMATTR Value V=V(x0) Laplace=exp(-s*ts)
SYMBOL bv 1024 416 R0
SYMATTR InstName B9
SYMATTR Value V=V(x1) Laplace=exp(-s*ts)
SYMBOL bv 0 -912 R0
SYMATTR InstName B10
SYMATTR Value V=V(x0) Laplace=ga*(s**2 + B1*s)/(s**2 + A1*s + A0)