Multisim Simulation of Strain Gage bridge compensation

* model for the stain gage
.SUBCKT STRAINGAGESIMF__STRAIN_GAGE__1__1 1 2 8 9 
.param R0=1000
.param TCGF=-2.0E-3
.param TCR=1e-3
.param GF0=100
EBUFF 3 0 1 2 1
R1 1 0 10G
R2 2 0 10G
R3 3 0 10G
EOPA 5 0 3 4 10G
RIN 4 0 1
RF 5 4 {GF0-1}
eout 7 69 poly(2) (5,0) (6,0) 0 0 0 0 1
fcopy 0 6 vsense 1
riz 5 0 10G
vsense 69 9 0
R4 7 8 {R0*(1+TCR*(TEMP-25))}
R5 6 0  {R0*(1+(TCR+TCGF)*(TEMP-25))}
.ENDS STRAINGAGESIMF

*model for the NTC thermistor
.SUBCKT NTCLE101E3C90173__DEF_1__1 1 2 Tin
R1 Tin 0 1T
Vsense 1 N002 0V
B3 N002 2 V= I(Vsense)*Rntc(V(Ttot),beta)
B1 N003 0 V= V(Tin)
B2 Ttot N003 V= +abs(V(vcth))
R2 Ttot 0 1T
I2 Beta 0 -1
R4 Vcth N001 1
C1 0 Vcth {Cth/Gth}
B4 N001 0 V= +abs(V(1,2)*I(Vsense))/Gth
.FUNC Rntc(x,y)
+ {if(x<25,+Rref*exp((y-Beta)*(1/(x+273.15)-1/(273.15+Tref)))*exp(B*(pwr(1/(273.15+x),1)-pwr(1/(273.15+Tref),1))+C*(pwr(1/(273.15+x),2)-pwr(1/(273.15+Tref),2))+D*(pwr(1/(273.15+x),3)-pwr(1/(273.15+Tref),3))),
+ +Rref*exp((y-Beta)*(1/(x+273.15)-1/(273.15+Tref)))*exp(B1*(pwr(1/(273.15+x),1)-pwr(1/(273.15+Tref),1))+C1*(pwr(1/(273.15+x),2)-pwr(1/(273.15+Tref),2))+D1*(pwr(1/(273.15+x),3)-pwr(1/(273.15+Tref),3))))}
.PARAM A=-12.89228328
.PARAM B=4245.148
.PARAM C=-87493
.PARAM D=-9588114
.PARAM a1=-12.89228328
.PARAM b1=4245.148
.PARAM c1=-87493
.PARAM d1=-9588114
.PARAM Gth=0.0068
.PARAM Rref=10000
.PARAM tolB=0.01
.PARAM tolRref=0.01
.PARAM Tref=25
.PARAM Cth=0.0476
.PARAM beta=3435
.IC V(vcth)=0
.ENDS NTCLE101E3C90173