Crate reimsum.m

Takes a row vector v representing the coefficients of a (n-1)st degree polynomial and three scalars a, b, n. It returns c1, which is the midpoint Reimann sum from a to b with n iterations and c2, which is the actual integral of the row vector v from a to b