function konvergenceMKD1d()

figure;
for i=1:7
    h = 10^(-i);
    [A,b] = mkd1d('sin',h);
    condest(A)
    u = A\b;
    
    hold on;
    n = length(u),
    plot(linspace(0,1,n+2),[0;u;0]');
    pause;
end

end

