function conv(c, n)

h = 1/n;

p1 = (-c/(2*h)) - (1/h^2);
p2 = (2/h^2);
p3 = (c/(2*h)) - (1/h^2);

M = zeros(n-1,3);
RHS = zeros(n-1,1);
RES = zeros(n-1,1);

M(:,1)=p1;
M(:,2)=p2;
M(:,3)=p3;


M(1,1) = 0;
M(n-1,3) = 0;

for K=1: n-1
    RHS(K,1)=3*(K/n)-6;
end

for K=2: n-1
    DZ = M(K,1)/M(K-1,2);
    M(K,1) = M(K,1) - DZ*M(K-1,2);
    M(K,2) = M(K,2) - DZ*M(K-1,3);
    RHS( K,1) = RHS(K,1) - DZ*RHS(K-1,1);
    RHS(K,1) = RHS(K,1)/M(K,2);
    M(K,:) = M(K,:)/M(K,2);

end

RHS(1,1) = RHS(1,1)/M(1,2);
M(1,:)=M(1,:)/M(1,2);
RES(n-1,1)=RHS(n-1,1);
for K=n-2:-1:1
    RES(K,1) = RHS(K,1)-M(K,3)*RES(K+1,1);
end

%wektor wartości od 0 do 1 z krokiem 1/n
PKT = [0:1/n:1];
PLT = zeros(1,n+1);

for K=2: n
    PLT(K) = RES(K-1);
end

plot(PKT,PLT,'rx');