function [U,x]=my_solution_gen(a,ua,b,ub,N,p,q);
x=uniform(a,b,N); % Omioimorfos diamerismos
h=x(2)-x(1); %yplogismos h
x_inside=x(2:end);
F=zeros(N,1);
for i=1:N
    F(i)=h^2*f(x_inside(i));
end
F(1)=F(1)+ua*(1-p(x(1))*h/2);
F(end)=F(end)+ub*(1+p(x(end))*h/2);
A=my_triag3(N,x_inside,q,p); %triafiagonios pinakas
y=A\F;  %epilisi grammikou sistimatos
U=zeros(N+2,1);
for i=1:N;
    U(i+1)=y(i);
end
U(1)=ua;U(end)=ub;

