%微分方程函數(shù)
function xdot=elecsys(t,x)  %狀態(tài)導(dǎo)數(shù)
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=[表達(dá)式1；表達(dá)式2；表達(dá)式3；...；表達(dá)式n-1]
%  表達(dá)式1 對(duì)應(yīng)  x1'=x2
%  表達(dá)式2 對(duì)應(yīng)  x2'=x3
%  表達(dá)式3 對(duì)應(yīng)  x3'=x4
%  ...
%  表達(dá)式n-1 對(duì)應(yīng)  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