#include "matutil.c"
#include <math.h>

double house(int n,vector x)
{
	int i;
	double s,t;
	s=sqrt(innerproduct(n,x,x));
	if(s!=0){
		if(x[0]<0) s=-s;
		x[0]+=s;t=1/sqrt(x[0]*s);
		for(i=0;i<n;i++) x[i]*=t;
	}
	return -s;
}