%Calculate v and find column
close all;
clc;
clear all;

load('M.mat');
t = M(:,1);
T2 = M(:,2);

v = ((T2(i)-T2(i+1))/(t(i+1)-t(i)));

k = find(T2 > 0 & v > 0.1);
columnT2 = k(1);