% msf7.m
%clear

load fnnwb7
%load plantv
t=0.05;
e=-2:0.05:2;
T0=0.005:0.005:20;
in=0.6104*sin(0.1*pi*T0);  %20000/32768=0.6104 is plant can track
T=in(1:100:4000);  

[n,N]=size(e);
%k3=50;
ec=e;
P=[e;ec];
NWB=wbg7;

i=1;
 W1=[NWB(i,1) 0;NWB(i,2) 0;NWB(i,3) 0;NWB(i,4) 0;NWB(i,5) 0;...
	0 NWB(i,6);0 NWB(i,7);0 NWB(i,8);0 NWB(i,9);0 NWB(i,10);...
	0 NWB(i,11);0 NWB(i,12);0 NWB(i,13);0 NWB(i,14)];
  B1=[NWB(i,15);NWB(i,16);NWB(i,17);NWB(i,18);NWB(i,19);NWB(i,20);...
	NWB(i,21);NWB(i,22);NWB(i,23);NWB(i,24);NWB(i,25);...
	NWB(i,26);NWB(i,27);NWB(i,28)];
  W2=[NWB(i,29) NWB(i,30) NWB(i,31) NWB(i,32) NWB(i,33) NWB(i,34)...
       NWB(i,35) NWB(i,36) NWB(i,37) NWB(i,38) NWB(i,39) NWB(i,40)...
       NWB(i,41) NWB(i,42) NWB(i,43) NWB(i,44) NWB(i,45) NWB(i,46)...
       NWB(i,47) NWB(i,48) NWB(i,49) NWB(i,50) NWB(i,51) NWB(i,52)...
       NWB(i,53) NWB(i,54) NWB(i,55) NWB(i,56) NWB(i,57) NWB(i,58)...
       NWB(i,59) NWB(i,60) NWB(i,61) NWB(i,62) NWB(i,63) NWB(i,64)...
       NWB(i,65) NWB(i,66) NWB(i,67) NWB(i,68) NWB(i,69) NWB(i,70)...
       NWB(i,71) NWB(i,72) NWB(i,73) NWB(i,74) NWB(i,75) NWB(i,76)...
       NWB(i,77)];

%Calculate the fitness values

% NETWORK PARAMETERS
 %W1 = [W10(:,1) zeros(5,1); zeros(5,1) W10(:,2)] ;

lfun=gause(W1*P,B1);
xp=multil(lfun);
cm=purelin(W2*xp);


plot(cm),ylabel('control value'),
pause,

%plot([e'*327.68 ec'*327.68]),
%ylabel('error'),

%plot([e(5:30)'*655.36 ec(5:30)'*655.36]),
%ylabel('error'),pause
plot(e,lfun(1,:)),
hold on
plot(e,lfun(2,:)),
plot(e,lfun(3,:)),
plot(e,lfun(4,:)),
plot(e,lfun(5,:)),
plot(e,lfun(6,:)),
plot(e,lfun(7,:)),

hold off,pause,

plot(e,lfun(8,:)),
hold on
plot(e,lfun(9,:)),
plot(e,lfun(10,:)),
plot(e,lfun(11,:)),
plot(e,lfun(12,:)),
plot(e,lfun(13,:)),
plot(e,lfun(14,:)),

hold