clear all;
init_variables;
time_step_init_eqns;
rhs_eqns;           
jacobian;
updates1=jac/rhs;

j_eqns;
lu_eqns;      
       
jp=j;
j_form;
jac-jform
%V(16+1)*Bus_Ymag(118+1)*cos(-theta(16+1)+theta(30+1)+Bus_Yang(118+1))*V(30+1)