%% Johdattelua interpolaatioon
% L2demo1 
clc
close all
xdata=[-2 -1 3]
ydata=[1 -2 5];
plot(xdata,ydata,'*')
%%
 A=[ones(3,1) [-2;-1;3] [-2;-1;3].^2]
 c=A\ydata'
 
%%
axis([-2.5 3.5 -5 10 ])
grid on
hold on
%%
x=linspace(-2.1,3.1,100);
p=c(1)+c(2)*x+c(3)*x.^2;
length(p)
%%
plot(x,p)
%%
plot(x,p,'r')   % Eri varilla.
%% Vektorien lineaarikombinaatio "riviajattelulla"
plot(x,p,'r')
%%
V=[ones(100,1) x' (x.^2)'];
y=V*c;
plot(x,y,'g')
title('Interpolaatiopoly 3:n pisteen dataan')
