disp("Wspolrzedne wezlow: ");
w1=[-5,5]
w2=[-3,1]
w3=[-2,4]
w4=[1,2]
w5=[2,-1]
w6=[4,-1]
A=[w1(1,1)^0,w1(1,1)^1,w1(1,1)^2,w1(1,1)^3,w1(1,1)^4,w1(1,1)^5;
w2(1,1)^0,w2(1,1)^1,w2(1,1)^2,w2(1,1)^3,w2(1,1)^4,w2(1,1)^5;
w3(1,1)^0,w3(1,1)^1,w3(1,1)^2,w3(1,1)^3,w3(1,1)^4,w3(1,1)^5;
w4(1,1)^0,w4(1,1)^1,w4(1,1)^2,w4(1,1)^3,w4(1,1)^4,w4(1,1)^5;
w5(1,1)^0,w5(1,1)^1,w5(1,1)^2,w5(1,1)^3,w5(1,1)^4,w5(1,1)^5;
w6(1,1)^0,w6(1,1)^1,w6(1,1)^2,w6(1,1)^3,w6(1,1)^4,w6(1,1)^5];
disp("Macierz A:");
disp(A);
B=[w1(1,2);w2(1,2);w3(1,2);w4(1,2);w5(1,2);w6(1,2)];
disp("Macierz B:");
disp(B);
C=inv(A)*B;
disp("Macierz C:");
disp(C);
fprintf(" W(x)=%f*x^5 + %f*x^4 + %f*x^3 + %f*x^2 + %f*x + %f \n",
C(6,1),C(5,1),C(4,1),C(3,1),C(2,1),C(1,1));
x=[-6:1:6];
y=zeros(1,13);
 for i=1:13

y(i)=C(1,1)+C(2,1)*x(i)+C(3,1)*x(i)^2+C(4,1)*x(i)^3+C(5,1)*x(i)^4+C(6,1)*x(i)^5;

 end
plot(x,y);
grid on;
hold on;
vector_x=[w1(1,1), w2(1,1), w3(1,1), w4(1,1), w5(1,1), w6(1,1)];
vector_y=zeros(1,6);
 for i=1:6
 vector_y(1,1)=vector_y(1,1) + C(i,1) * w1(1,1)^(i-1);
 vector_y(1,2)=vector_y(1,2) + C(i,1) * w2(1,1)^(i-1);
 vector_y(1,3)=vector_y(1,3) + C(i,1) * w3(1,1)^(i-1);
 vector_y(1,4)=vector_y(1,4) + C(i,1) * w4(1,1)^(i-1);
 vector_y(1,5)=vector_y(1,5) + C(i,1) * w5(1,1)^(i-1);
 vector_y(1,6)=vector_y(1,6) + C(i,1) * w6(1,1)^(i-1);
 end
plot(vector_x,vector_y,'bs');
plot(w1(1,1),w1(1,2),'-r*');
plot(w2(1,1),w2(1,2),'-r*');
plot(w3(1,1),w3(1,2),'-r*');
plot(w4(1,1),w4(1,2),'-r*');
plot(w5(1,1),w5(1,2),'-r*');
plot(w6(1,1),w6(1,2),'-r*');
title("Interpolacja wielomianowa");
ylabel("Os y");
xlabel("Os x");