% 編寫求解起動過程微分方程的腳本函數
% 將該腳本函數定義為sh_motor_vary_so(shunt_motor_self_excited_solver)
[t,y]=ode23(@sh_motor_vary_ode,[0:.01:10],[20.8332996141; 872.6049312471;   63.8736361288])
% 首先繪制if=f(t)曲線
subplot(3,1,1),
plot(t,y(:,1))
xlabel('Time[T]')
ylabel('If[A]')
grid on
% 下面繪制ia=f(t)曲線
subplot(3,1,2),
plot(t,y(:,2)),
xlabel('Time[T]')
ylabel('Ia[A]')
grid on
% 下面繪制w=f(t)曲線
subplot(3,1,3),
plot(t,y(:,3))
xlabel('Time[T]')
ylabel('w[s/s]')
grid on