Note

In the theory the data points are counted through . Inside the program you will need to count the data points from 1 to m. This can be in principle avoided by specifying the subscript range for the vectors.