function xdot=ch2example10statefun(t, x, flag, R,L,C)
% 考慮RLC串聯環路的狀態方程函數
% 輸入: t 當前計算時刻,  flag 此處不用
% x為2#1矩陣, x(1)為電感上的電流; x(2)為電容電壓
% R: 電阻值
% L: 電感值
% C: 電容值
xdot=zeros(2,1);        % 狀態變量矩陣初始化
A=[-R/L, -1/L; 1/C, 0]; % 狀態方程系數矩陣
B=[1/L; 0]; 
xdot=A*x+B*f(t);        % 狀態方程
function input=f(t)
input=(t>=0);           % 輸入信號為單位階躍