MatN_Lab2_ex1.3

clear
close
clc

n=input('Hello, input number: ')
h=1/(n-1);
D=-1*eye(n);

for i=1:n-1
    D(i,i+1)=1;
end

D

b=zeros(n,1);

hold on
for k=(1:n)
    b(k)=h*cos(2*pi*(k-1)/(n-1));
    x=D\b;
end

plot(1:n,x,'b-');

xlabel('k');
ylabel('x_k');
title('LES solution');