function [y]=linear_interp(x_vector,y_vector,desired_x_value)
%LINEAR_INTERP Função que executa a interpolação linear
%   Esta função interpola linearmente os valores do vector y_vector, em
%   função dos valores do vector x_vector, usando interpolação linear, e a
%   ordenada dada em desired x_value.
%
%   INTPUTS:
%   x_vector -> Vector de absissas;
%   y_vector -> Vector de ordenadas;
%   desired_x_value -> Valor em torno do qual se pretende fazer a
%                      interpolação;
%   OUTPUTS:
%   y -> Ordenada resultante da interpolação.
%##########################################################################
if (desired_x_value>x_vector(1))&&(desired_x_value<x_vector(size(x_vector,2)))
    for i=1:1:size(x_vector,2)
        if x_vector(i)>desired_x_value
            xb=x_vector(i);
            xa=x_vector(i-1);
            yb=y_vector(i);
            ya=y_vector(i-1);
            y=ya+(yb-ya)*((desired_x_value-xa)/(xb-xa));
            break
        end
    end
else
    y=0;
end
end