a1=1.2236,a2=1.4642,a3=0.2074,a4=1.2537,a5=0.4394,a6=0.6242,a7=1.4191,a8=0.3836,a9=1.9154,a10=0.9443,a11=1.4965,a12=2.2248,
%a1=0.5812,a2=1.5760,a3=0.2203,a4=1.9863,a5=1.5416,a6=0.8068,a7=0.7682,a8=0.6435,a9=1.6297,a10=0.9486,a11=2.2162,a12=1.3504,
%a1=0.5383,a2=1.1635,a3=0.1623,a4=1.4599,a5=2.0658,a6=0.6715,a7=0.7918,a8=1.4105,a9=0.2826,a10=1.2193,a11=1.7092,a12=2.2914,
%a1=1.4750,a2=1.4814,a3=0.2482,a4=1.2429,a5=0.9723,a6=0.2805,a7=1.6941,a8=1.6533,a9=1.6920,a10=0.7767,a11=1.2623,a12=2.0830,
%a1=1.7134,a2=1.5179,a3=0.1236,a4=1.2816,a5=0.5361,a6=1.3547,a7=1.4814,a8=1.4256,a9=0.9916,a10=0.7510,a11=2.0658,a12=2.2355,
%a1=1.4728,a2=2.1603,a3=0.3900,a4=2.2935,a5=1.4428,a6=1.9218,a7=1.8810,a8=1.2795,a9=2.0185,a10=0.2547,a11=1.8123,a12=1.5351,
%a1=2.0379,a2=0.6307,a3=0.1623,a4=0.4137,a5=0.4051,a6=0.7209,a7=1.9412,a8=1.8982,a9=0.9615,a10=0.2955,a11=0.9271,a12=1.5974,
%a1=1.7134,a2=1.5179,a3=0.1236,a4=1.2816,a5=0.5361,a6=1.3547,a7=1.4814,a8=1.4256,a9=0.9916,a10=0.7510,a11=2.0658,a12=2.2355,
%a1=0.5383,a2=1.1635,a3=0.1623,a4=1.4599,a5=2.0658,a6=0.6715,a7=0.7918,a8=1.4105,a9=0.2826,a10=1.2193,a11=1.7092,a12=2.2914,
%a1=0.1623,a2=0.7230,a3=0.3148,a4=0.7896,a5=1.6941,a6=1.8638,a7=0.8305,a8=0.2783,a9=0.5340,a10=0.3449,a11=1.1398,a12=1.3998,
%a1=1.4728,a2=2.1603,a3=0.3900,a4=2.2935,a5=1.4428,a6=1.9218,a7=1.8810,a8=1.2795,a9=2.0185,a10=0.2547,a11=1.8123,a12=1.5351,
%a1=0.6586,a2=1.6297,a3=0.4609,a4=0.7359,a5=1.2945,a6=1.7006,a7=2.0722,a8=1.2924,a9=2.1066,a10=0.7531,a11=1.2494,a12=2.1131,
%a1=0.5,a2=0.5,a3=0.5,a4=0.5,a5=0.5,a6=0.5,a7=0.5,a8=0.5,a9=0.5,a10=0.5,a11=0.5,a12=0.5,
a=newfis('nm');
a=addvar(a,'input','e',[-4,4]);
a=addmf(a,'input',1,'NB','gaussmf',[a1,-4]);
a=addmf(a,'input',1,'NS','gauss2mf',[a1,-2,a2,-2]);
a=addmf(a,'input',1,'ZO','gauss2mf',[a2,0,a3,0]);
a=addmf(a,'input',1,'PS','gauss2mf',[a3,2,a4,2]);
a=addmf(a,'input',1,'PB','gaussmf',[a4,4]);
a=addvar(a,'input','ec',[-4,4]);
a=addmf(a,'input',2,'NB','gaussmf',[a5,-4]);
a=addmf(a,'input',2,'NS','gauss2mf',[a5,-2,a6,-2]);
a=addmf(a,'input',2,'ZO','gauss2mf',[a6,0,a7,0]);
a=addmf(a,'input',2,'PS','gauss2mf',[a7,2,a8,2]);
a=addmf(a,'input',2,'PB','gaussmf',[a8,4]);
a=addvar(a,'output','u',[-4,4]);
a=addmf(a,'output',1,'NB','gaussmf',[a9,-4]);
a=addmf(a,'output',1,'NS','gauss2mf',[a9,-2,a10,-2]);
a=addmf(a,'output',1,'ZO','gauss2mf',[a10,0,a11,0]);
a=addmf(a,'output',1,'PS','gauss2mf',[a11,2,a12,2]);
a=addmf(a,'output',1,'PB','gaussmf',[a12,4]);
f1=[1 1 1 1 1];
f2=[1 2 1 1 1];
f3=[1 3 2 1 1];
f4=[1 4 2 1 1];
f5=[1 5 3 1 1];
f6=[2 1 1 1 1];
f7=[2 2 2 1 1];
f8=[2 3 2 1 1];
f9=[2 4 3 1 1];
f10=[2 5 3 1 1];
f11=[3 1 2 1 1];
f12=[3 2 2 1 1];
f13=[3 3 3 1 1];
f14=[3 4 3 1 1];
f15=[3 5 4 1 1];
f16=[4 1 2 1 1];
f17=[4 2 3 1 1];
f18=[4 3 3 1 1];
f19=[4 4 4 1 1];
f20=[4 5 4 1 1];
f21=[5 1 3 1 1];
f22=[5 2 3 1 1];
f23=[5 3 4 1 1];
f24=[5 4 4 1 1];
f25=[5 5 5 1 1];
rulelist=[f1;f2;f3;f4;f5;f6;f7;f8;f9;f10;f11;f12;f13;f14;f15;f16;f17;f18;f19;f20;f21;f22;f23;f24;f25];
a=addrule(a,rulelist);
a=setfis(a,'DefuzzMethod','mom');
    sim('nmcc');
    y=ScopeData1.signals.values;
    t=ScopeData1.time;
    e=abs(50-y);
    eval=trapz(t,e),