StI=load('RepStor.txt');
It=load('RepInf.txt');
Rwt=load('RangeWatDem.txt');
Rit=load('RangeIrrDem.txt');
EnDem=load('RangeHyDem.txt');

%for j=1:12
    for u=1:30
        for k=1:10
            for n=1:6
                for m=1:11
                    for p=1:11
                  %StF(j,u,k,n,m,p)=StI(u)+It(k,j)-Rwt(n)-EnDem(m)-Rit(p);
		StF1(u,k,n,m,p)=StI(u)+It(k,1)-Rwt(n)-EnDem(m)-Rit(p);
		StF2(u,k,n,m,p)=StI(u)+It(k,2)-Rwt(n)-EnDem(m)-Rit(p);
		....
		StF12(u,k,n,m,p)=StI(u)+It(k,12)-Rwt(n)-EnDem(m)-Rit(p);
                    end
                end
            end
        end
    end
%end