%Parser V1.0
%Input Spice Deck:  C:\Users\sudhirk\Desktop\2T2Raltside.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 ;

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) + (X(3))/r2 ;

F(3) =  (X(1)) -12;

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

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

%                      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 ;

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 ;

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) ;

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)))) + 0 ;

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)))) + 1/r2 ;

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

JAC(3, 2) = 0;

JAC(3, 3) = 0;

JAC(3, 1) = 1;
