% 編寫求解空載建立電壓微分方程的函數(sy_ge_damp_noload_sol)
function sy_ge_damp_noload_sol 
%設定時間和初始值
tspan = [0 10];
y0 = [0;0;0;0;0];
% 方程求解
[t,y]=ode113('tongbufadianji',tspan,y0);
%ic=cos(t+2*pi/3).*y(:,1)-sin(t+2*pi/3).*y(:,2);
plot(t,y(:,3))

% 將該腳本函數定義為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)曲線