The block above could benefit from a  few comments. E.g. setting up data matrix + results vector + sat something about the shape required etc.