function B = extract(T, A, tArray)
%  extract(T, A, tArray);
%  get values of function A(T)
%  corresponing to argument values tArray

L = length(T);
iA = 1;
B = [];

for i = 1 : length(tArray)
    tt = tArray(i);
    while T(iA) < tt && iA < L
        iA = iA + 1;
    end;
    B = [B, A(iA)];
end;
