%RSCJ equations to do iv, with non linear subgap resistance
%
% v is in units of IcRn
% t is in units of t_J
% so that this is the Method 1 formulation 
% rnonl is the r(v), the non linear resistancs.
% from jjnlin units
%
function xdot = jjnlin(t,x)
global  betac Ib Rs Rsg;
xdot(1) = (-sin(x(2))./(betac)) - x(1)./(betac .* rnonl(x(1),Rsg))+ Ib/(betac); 
xdot(2) = x(1);

