function[x,u] = uexakt(a,b,N)

h = (b-a)/N;

for i = 1:N+1
    x(i) = a + (i-1)*h;
    u(i) = 1/(log(x(i)*x(i)+1)+0.5);
end