% S-函數的編程推薦結構
function [sys,x0,str,ts] = Sfunname(t,x,u,flag) % S-函數接口
% ---根據Flag進行分支---%
switch flag,                           % 根據傳入的調用標志再調用相應的處理子函數
  case 0,                              % 如果標志要求初始化
    [sys,x0,str,ts]=mdlInitializeSizes;% 則調用mdlInitializeSizes函數做初始化
  case 1,                              % 如果要求計算連續狀態方程
    sys=mdlDerivatives(t,x,u);         % 則調用mdlDerivatives函數：微分方程計算
  case 2,                              % 如果要求計算離散狀態方程
    sys=mdlUpdate(t,x,u);              % 則調用mdlUpdate函數：差分方程遞推更新
  case 3,                              % 如果要求計算系統輸出
    sys=mdlOutputs(t,x,u);             % 則調用mdlOutputs函數
  case 4,                              % 如果要求計算下一步仿真的時刻 
    sys=mdlGetTimeOfNextVarHit(t,x,u); % 則調用mdlGetTimeOfNextVarHit函數
  case 9,                              % 如果仿真終止     
    sys=mdlTerminate(t,x,u);           % 則調用mdlTerminate函數進行最后處理
  otherwise                                    % 否則
    error(['Unhandled flag = ',num2str(flag)]);% 進行出錯處理，返回錯誤代碼
end 
% ---sfuntmpl 函數結束---
% 下面是各個子函數的具體實現代碼
funcion  [sys,x0,str,ts]=mdlInitializeSizes
(以下省略)