load fcmdata3
trn_data=fcmdata3(1:100,:);
chk_data=fcmdata3(101:200,:);
in_fismat=genfis1(trn_data);
ss=0.5;
ss_dec=0.9;
ss_inc=1.1;
epoch=300;
[fismat1,error1,stepsize,fismat2,error2]=...
	anfis(trn_data,in_fismat,[epoch nan ss ss_dec ss_inc],...
	[1 1 1 1],chk_data);

figure(1);
plot([error1,error2]);
title('RMSE curve from moisture and its diff. of two day ago.');
ylabel('RMSE');
xlabel('N');

y=[trn_data(:,5);chk_data(:,5)];
yy=evalfis([trn_data(:,1:4);chk_data(:,1:4)],fismat1);
figure(2);
plot([y yy]);
title('Actual and Predicted Rainfall (fcm4.m) ');
ylabel('Rainfall');
xlabel('N');