function [u,j,f]=oofcmf(r,e,ce,k)%r為給定的值，e＝r，ce＝0,k=0
if(ce >= 1)                   %遞歸出口
    {
        f=1/1+j;
        return f ; 
    } 
end    
x=0;                        %設置規則：if x is Ai and y is Bi Then u is Ui
y=0;
a=0;
j=0;
u=[1,2];                    %Ui 用一個矩陣來存儲  
for i=1:2                   %設置三角形隸屬函數的寬度
    b=i;
    for h=2:2:4
        c=h;
        x+=trig(e,a,b,c)*trig(ce,a,b,c)*u(i);           %利用計算公式：U(k)=""
        y+=trig(e,a,b,c)*trig(e,a,b,c);
    end
end
u=x/y;                     %第k時刻u(k)的值
j+=u;                      %累加j      
ts=0.001;                  %采樣時間
y_1=0;                     %有疑問
y_2=1;
u_1=0;
u_2=1;
dsys=c3d(sys,ts,'z');
[num,den]=tfdata(dys,'v');
yout(k)=-den(2)*y_1-den(3)*y_2+num(2)*u_1+num(3)*u_2;   %如何知道u_1,u_2,y_1.y_2
y_2=y_1;
y_1=y;
u_2=u_1;
u_1=u;
x=e;
e=r-y;
ce=e-x;
k=k+1;                  
oofcmf(r,e,ce,k);%遞歸調用主函數，計算第 k+1時刻