Add simple convolution with point spread function
add blur_matrix
add helper method
add mmult
add mmult and vmult
add len
allow plotting of several functions
bugfix, add matrix inverse calculation