load TEC %Valores de vTEC
load grilla %Todos los puntos grilla
pGrilla=182; %Cantidad de puntos grilla
vEC=[];
sd=[]; %Conjunto de las simples diferencias
dd=[]; %Conjunto de las dobles diferencias
interm=[]; %Conjunto intermedio de valores
for i=1:pGrilla %Recorro toda la grilla
k=find(TEC(:,1)==grilla(i,1) & TEC(:,2)==grilla(i,2)) %Cambio cada punto grilla
interm=TEC(k,:) %Vector que va almacenando valores
vEC(:,i)=interm(:,3) %valor de vTEC tomado
sd(:,i)=diff(interm(:,3)); %Simple diferencia
dd=diff(sd); %Doble diferencia
end