x = [ 1.234 , 2.432 , 0.542 , -0.342 , -2.003 , - 3.678 ] ;
y =  [ 5.423 , 6.423 , -0.234 , -4.078 , - 8.762 , 3.543 ] ;

v=x';
A= [ v.^5 , v.^4 , v.^3 , v.^2 , v , ones(6,1) ] 
    
inv(A)*y.'


%Her er en enda enklere måte å lage til matrisen på.  Det er benytten en
%liten løkke.  Det er noe dere lærer mer om senere. 

V=x';
A = ones(6,1);
M=A;
for i = 1:5
  M = [ M .*V , A ];
  V=[V,x'];
end

v=inv(M)*y';

%Her har vi tatt med et plot av polynomet sammen med et plot av punktene.
syms z 
p = @(z)  v(1)*z.^5 + v(2)*z.^4 + v(3)*z.^3  + v(4)*z.^2 + v(5)*z + v(6) ;
plot(x,y, 'o r' ) 
hold on
w=-4:0.01:3;
plot(w, p(w) )
hold off

%Her er ogs\aa\ matrisen til del a) ...
C= [  16  , -8 , 4 , -2 , 1 , 3 ; 1 , -1 , 1 , -1 , 1 ,1 ; 1 , 1 , 1 , 1 , 1 , 2 ; 16  , 8 , 4 , 2 , 1  , 11 ];
    