%微分方程函數
function xdot=elecsys(t,x)  %狀態導數
Vs=1;
R=1.4;
L=2;
C=0.32;
xdot=[x(2)/C;1/L*(Vs-x(1)-R*x(2))];

%格式
%  function xdot=filename(t,x)
%     xdot=[表達式1；表達式2；表達式3；...；表達式n-1]
%  表達式1 對應  x1'=x2
%  表達式2 對應  x2'=x3
%  表達式3 對應  x3'=x4
%  ...
%  表達式n-1 對應  xn-1'=xn
%	如例exp3_1.m
%	x(1)=Vo   x(2)=iL
% x(1)'=x(2)/C
% x(2)'=(Vs-x(1)-R*x(2))/L