function U=my_solution(a,b,N);
x=uniform(a,b,N); % Omioimorfos diamerismos
h=x(2)-x(1); %yplogismos h
F=zeros(N,1);
for i=1:N
    F(i)=h^2*f(x(i+1));
end
A=my_triag(N); %triafiagonios pinakas
y=A\F;
U=zeros(N+2,1);
for i=1:N;
    U(i+1)=y(i);
end

