lab3

clc
a = 0;
b = 1;
N = 5;

x = linspace(a, b, N+1);
h = x(2)-x(1);


A = zeros(N+1);
B = zeros(N+1, 1);

for i = 1:N+1
    for j = 1:N+1
        delta = 0;
        if i == j
            delta = 1;
        end
        c = h;
        if j == 1 || j == N+1
            c = h/2;
        end
        A(i, j) = delta - c*k(x(i), x(j));
    end
    B(i) = f(x(i));
end
A

% c=h*ones(N+1,1);
% c(1)=h/2;
% c(N+1)=h/2;
% for i = 1:N+1
%     for j = 1:N+1
%           A1(i, j) = - c(j)*k(x(i), x(j));
%     end
%     B(i) = f(x(i));
%     A1(i,i)=1+A1(i,i);
% end
% A1
% A-A1
%

V = A\B
plot(x,V)