%
% 分岔圖(Fork)
% 
% ****** Author:J.F.Peng(jingyujiafu@163.com) ******
%
clear
clc
%
% tic、toc用來計算程序運行時間，分別表示開始和結束計時
%
tic
global Bc ia w id; 
id=0.00; 
w=0.20; 
ia=1.50;
for Bc=0.00:0.002:2
    [T,Y]=ode45('Josephson_Junction',[0,200],[0;0;0]); 
    data=Y(:,2); 
    n=length(data); 
    % 去掉一些前面的點，防止因計算機配置問題而畫不出圖，不影響結果
    m=round(n/6); 
    a=data(m-2); 
    b=data(m-1); 
    for i=m:n 
        if b>=a&b>=data(i); 
           plot(Bc,b,'k'); 
           hold on; 
        end 
        a=b; 
        b=data(i); 
    end 
end 
%
% 根據需要選擇圖像標注語句
%
% ia-u 分岔圖（fork）
%
%xlabel('Aternative Current ia')
%ylabel('Voltage u=dφ/dt')
%desc={'Bc=0.50,w=0.66,id=0.00'};
%text(0.5,3,desc)
%
% id-u 分岔圖 (fork)
%
%xlabel('Directive Current id')
%ylabel('Voltage u=dφ/dt')
%desc={'Bc=1.00,w=0.50,ia=0.65'};
%text(1,3,desc)
%
% Bc-u 分岔圖 (fork)
%
xlabel('Damping Coefficient \betac')
ylabel('Voltage u=dφ/dt')
desc={'ia=1.50,w=0.20,id=0.00'};
text(0.75,8,desc)
toc
t=toc