% rload.m
% 實現一個阻性負載

function[sys,xo]=rload(t,x,u,flag,R)
% 由于不是離散系統，輸出參數沒有包括ts
% 此外需要用戶輸入1個新的參數
% 負載電阻R

if flag==0
   sys=[0 0 2 1 1 1];
   % 0個狀態
   % 2個輸出，sys(1)為負載電流；sys(2)為輸出電壓
   % 1個輸入，u為交流輸入電壓
   % 初始狀態為空
   xo=[];
elseif flag==1
   sys=0;
elseif flag==3
   sys(1)=u(1)/R;
   sys(2)=u(1);
   else
   sys=[];
   end