%
function dy=Josephson_Junction(t,y)
% 在程序中參數t代表文中的τ
% Bc 為阻尼參數, w 為歸一化的外加交流電流圓頻率
% Bc、ia、id、w、τ均為無量綱化時引入的變量，其中ia、id為各電流無量綱化時引入的
% y(1)= φ  約瑟夫森結中心處波函數ψ1、ψ2的相位差
% y(2)=dφ/dt 即歸一化電壓u
% 以下為約瑟夫森結RCSJ模型的無量綱化微分方程
% dφ/dt=u
% du/dt=-Bc*u-sin(φ)+id+ia*sin(w*t)
%
% ****** Author:J.F.Peng(jingyujiafu@163.com) ******
%
global Bc ia w id;
dy=zeros(3,1);
dy(1)=y(2);
dy(2)=id+ia*sin(w*t)-sin(y(1))-Bc*y(2);
dy(3)=w;