function fangzhen 
%RIGIDODE  Euler equations of a rigid body without external forces
tspan = [0  1];
y0 = [17.8876/0.13;0;0;0];
% Solve the problem using ode45
%options = odeset('MaxStep',1e-6)

[t,y]=ode23tb('tongbufadianji',tspan,y0);
%ic=cos(t+2*pi/3).*y(:,1)-sin(t+2*pi/3).*y(:,2);
subplot(3,1,1)
plot(t,y(:,1))
hold on 
subplot(3,1,2)
plot(t,y(:,2))
hold on
subplot(3,1,3)
plot(t,y(:,4))
hold on
%ib=y(:,1).*sin(-120)+y(:,2).*cos(120);
%plot(t,ib)
% 編寫求解自激過程微分方程的腳本函數
% 將該腳本函數定義為sh_ge_se_ex_so(shunt_generator_self_excited_solver)
%[t,iff]=ode23(@sh_ge_se_ex_ode,[0 10],[0])
% 指定全局變量（注意全局變量必須在同時使用的多個函數中同時指定）
% 繪圖（這里采用多子圖法同時繪制ea=f(if)、uf=f(if)、ea=f(t)、if=f(t)等多條曲線）
% 首先繪制ea=f(if)和uf=f(if)曲線
%subplot(3,1,1),
% 下面指定縱橫軸標簽
xlabel('Time[s]')
ylabel('If[A]')
% 下面繪制if=f(t)曲線