Added basic version of linear solution

NB! Still using inaccurate right boundary condition (du/dx = 0)!