function [ I ] = LinSpline( x,y )
%Calculates coefficients of spline

n = length(x);

plot(x,y, 'or');
hold on

for i = 1:n-1
    A = [1 x(i); 1 x(i+1)];
    b = [y(i); y(i+1)];
    a = A \ b;
    
    xrange = linspace(x(i), x(i+1),10);
    line = a(2)*xrange + a(1);
    plot(xrange, line);
    hold onp
end
    
end