file proc qr.prc
file proc hh.prc 
* solve a least squares problem by HH 
gen1 N_=3
gen1 K_=2
read X_ / rows=N_ cols=K_ list
3 -1
0 0
4 7
read y / rows=N_ cols=1 list
0
18
25
set nodoecho
EXEC QR
stop