Merged in feature/dirac_discr (pull request #17) Add multi-d dirac discretization with tests Approved-by: Vidar Stiernström <vidar.stiernstrom@it.uu.se> Approved-by: Martin Almquist <malmquist@stanford.edu> Approved-by: Jonatan Werpers <jonatan.werpers@it.uu.se>

% [x,h] = get_grid(a,b,m)
function [x,h] = get_grid(a,b,m)
    % TODO: allow the interval to be a vector
    x = linspace(a,b,m)';
    h = (b-a)/(m-1);
end