A=[1 0 0 0 0; -0.5 1 -0.5 0 0; 0 -0.5 1 -0.5 0; ...
    0 0 -0.5 1 -0.5; 0 0 0 -1 1];
D = [0 -0.1 -0.1 -0.1 0];
U = inv(A)*D'
plot(U)
xlabel('x')
ylabel('u')

    
    




