%Parser V1.0
%Input Spice Deck:  C:\Users\sudhirk\Downloads\Intern_2016\ckt_inout\2T2Rsameside.txt
%Equation Type:     NODAL
%Datum Node:        2

%*****************************************************************************
%                      Component Values:
q1IS = 1e-016;
q1BF = 100;
q1BR = 1;
q1N = -1.25262e-058;
q2IS = 1e-016;
q2BF = 100;
q2BR = 1;
q2N = -1.25262e-058;
r1 = 1;
r2 = 2;

%*****************************************************************************

%                      Circuit Equations: 

F(1) =  (q1IS )*(exp(-q1N*(-X(3))) -1) + (-q1IS /q1BR )*(1+q1BR )*(exp(-q1N*(X(2)-X(3))) -1) + (q2IS )*(exp(q2N*(X(1)-X(2))) -1) + (-q2IS /q2BR )*(1+q2BR )*(exp(q2N*(X(3)-X(2))) -1) + (-q2IS /q2BF )*(1+q2BF )*(exp(q2N*(X(1)-X(2))) -1) - (-q2IS )*(exp(q2N*(X(3)-X(2))) -1) + (X(2)-X(1))/r1 + (X(2))/r2 ;

F(2) =  (-q1IS )*(exp(-q1N*(-X(3))) -1) + (q1IS /q1BR )*(1+q1BR )*(exp(-q1N*(X(2)-X(3))) -1) + (q1IS /q1BF )*(1+q1BF )*(exp(-q1N*(-X(3))) -1) - ( q1IS )*(exp(-q1N*(X(2)-X(3))) -1) + (-q2IS )*(exp(q2N*(X(1)-X(2))) -1) + (q2IS /q2BR )*(1+q2BR )*(exp(q2N*(X(3)-X(2))) -1) ;

F(3) =  (q2IS /q2BF )*(1+q2BF )*(exp(q2N*(X(1)-X(2))) -1) - ( q2IS )*(exp(q2N*(X(3)-X(2))) -1) + (X(1)-X(2))/r1 ;

F(4) =  (-12) ;

%*****************************************************************************

%                      Jacobians: 

JAC(1, 2) =  (q1IS * q1N / q1BR)*(1 + q1BR)*(exp( - q1N*(X(2) - X(3)))) + ( - q2IS * q2N)*(exp(q2N*(X(1) - X(2)))) + (q2IS * q2N / q2BR )*(1 + q2BR )*(exp(q2N*(X(3)-X(2)))) + (q2IS * q2N / q2BF )*(1 + q2BF )*(exp(q2N*(X(1) - X(2)))) - (q2IS * q2N)*(exp(q2N*(X(3) - X(2)))) + 1/r1 + 1/r2 ;

JAC(1, 3) =  (q1IS * q1N)*(exp( - q1N*( - X(3))))  - (q1IS * q1N / q1BR)*(1 + q1BR)*(exp( - q1N*(X(2)-X(3)))) + ( - q2IS * q2N / q2BR)*(1 + q2BR)*(exp(q2N*(X(3) - X(2)))) + (q2IS * q2N)*(exp(q2N*(X(3) - X(2)))) + 0 + 0 ;

JAC(1, 1) =  0 + (q2IS * q2N)*(exp(q2N*(X(1) - X(2)))) - (q2IS * q2N / q2BF)*(1 + q2BF)*(exp(q2N*(X(1) - X(2)))) + (-1/r1) + 0 ;

JAC(1, 4) =  0 + 0 + 0 + 0 ;

JAC(2, 2) =  ( - q1IS * q1N / q1BR)*(1 + q1BR)*(exp(-q1N*(X(2) - X(3)))) + (q1IS * q1N)*(exp(-q1N*(X(2) - X(3)))) + (q2IS * q2N)*(exp( q2N*(X(1) - X(2))))  - (q2IS * q2N / q2BR)*(1 + q2BR)*(exp( q2N*(X(3)-X(2)))) ;

JAC(2, 3) =  ( - q1IS * q1N)*(exp(-q1N*( - X(3)))) + (q1IS * q1N / q1BR )*(1 + q1BR )*(exp(-q1N*(X(2)-X(3)))) + (q1IS * q1N / q1BF )*(1 + q1BF )*(exp(-q1N*( - X(3)))) - (q1IS * q1N)*(exp(-q1N*(X(2) - X(3)))) + (q2IS * q2N / q2BR)*(1 + q2BR)*(exp( q2N*(X(3) - X(2)))) ;

JAC(2, 1) =  0 + ( - q2IS * q2N)*(exp( q2N*(X(1) - X(2)))) ;

JAC(2, 4) =  0 + 0 ;

JAC(3, 2) =  ( - q2IS * q2N / q2BF)*(1 + q2BF)*(exp(q2N*(X(1)-X(2)))) + (q2IS * q2N)*(exp(q2N*(X(3) - X(2)))) + (-1/r1) ;

JAC(3, 3) =  ( -  q2IS * q2N)*(exp(q2N*(X(3) - X(2)))) + 0 ;

JAC(3, 1) =  (q2IS * q2N / q2BF)*(1 + q2BF)*(exp(q2N*(X(1) - X(2)))) + 1/r1 ;

JAC(3, 4) =  0 + 0 ;

JAC(4, 2) =  0 ;

JAC(4, 3) =  0 ;

JAC(4, 1) =  0 ;

JAC(4, 4) =  0 ;
