%Determina el indice Esfuerzo de control CE

ws=input('Introduce variable a calcular incremento \n','s');
ws1=input('Introduce columna \n','s');
ws1=str2num(ws1);
ws=eval(ws);
b=zeros(7212,1);
t=zeros(7212,1);
CE=0;
eU=zeros(7212,1);
for i=1:7211;
    b(i+1,1)=ws(i,ws1)-ws(i+1,ws1);
end;
ws2=input('Introduce valor estacionario valvula 1 \n','s');
ws2=str2num(ws2);
ws3=input('Introduce valor estacionario valvula 2 \n','s');
ur1=ws2;
ws3=str2num(ws3);
ur2=ws3;
eU=abs([b(1:3609,1)/ur1;b(3610:7212,1)/(ur2-ur1)]);
t=[Tsm(1:3609,1);Tsm(3610:7212,1)-3600];
CE=t'*eU
