#include <oxstd.h>

MyOls(const mY, const mX)
{
    decl b;
    olsc(mY, mX, &b);
    return b;
}
main()
{
    decl y  = <0.58; 1.1; 1.2;  1.3;  1.95;
               2.55; 2.6; 2.9;  3.45; 3.5;
               3.6;  4.1; 4.35; 4.4;  4.5>;
    decl x1 = <1 : 15>';    // note transpose!
    decl mx = 1 ~ x1 ~ x1 .^ 2;  // regressors
    print(y ~ mx);           // print all data

// 2do: use MyOls to regress y on mx & print results
}
