function xdot=rcstateequation(t, x, flag, par)
% RC串聯電路的狀態方程函數
% 輸入: t 當前計算時刻, x 當前狀態－電容電壓, flag 此處不用
% par 為2#1矩陣, par(1)為電阻值參數, par(2)為電容值參數
R=par(1); % 電路參數：電阻值
C=par(2); % 電路參數：電容量
xdot=(-1./(R*C)).*x+(1/(R*C)).*f(t); % 狀態方程實現
%---以下是用來產生輸入信號f(t)的子函數---
function inputsignal=f(t) 
inputsignal=(t>=0);% 產生的輸入信號為單位階躍