void QR_dec(real **,int,int,real **);
void QR_back(real **,real **,int,int,real *);