p17 oxtut2e long.ox

时间序列分析中著名的OxMetrics软件包

#include <oxstd.h>

MyOls(const mY, const mX)
{ decl b;
b= (mX'mX)^-1 * (mX'mY);
print("in MyOls(): b=", b);
return b;
} 

main()
{ decl b;
// mY, mX both just random
b = MyOls( rann(10, 1), ranu(10, 2) );
// b now holds the result
print("in main(): b=", b);
}