a = newfis('hw3ans');
a = addvar(a,'input','psi',[-45 45]); 
a = addmf(a,'input',1,'LN','trapmf', [-135 -45 -30 -14]);
a = addmf(a,'input',1,    ,        , [                ]);
a = addmf(a,'input',1,    ,        , [                ]);
a = addmf(a,'input',1,    ,        , [                ]);
a = addmf(a,'input',1,    ,        , [                ]);
a = addvar(a,'input','yA',[-5 5]);
a = addmf(a,'input',2,'LN','trapmf', [-10 -9 -1 -0.5]);
a = addmf(a,'input',2,    ,        , [              ]);
a = addmf(a,'input',2,    ,        , [              ]);
a = addmf(a,'input',2,    ,        , [              ]);
a = addmf(a,'input',2,    ,        , [              ]);
a = addvar(a,'output','beta',[-70 70]);
a = addmf(a,'output',1,'LN','trimf',[-70 -45 -20]);
a = addmf(a,'output',1,    ,       , [           ]);
a = addmf(a,'output',1,    ,       , [           ]);
a = addmf(a,'output',1,    ,       , [           ]);
a = addmf(a,'output',1,    ,       , [           ]);

ruleList = [1 1 5 1 1
            1 2 5 1 1  % etc.
            
                ];
a = addrule(a,ruleList);
tune_rules.impMethod = 'prod';