double xgaus[]={-0.3333333, -0.0597158717,-0.0597158717, -0.8805682564, -0.7974269853, -0.7974269853, 0.5948539707};
double ygaus[]={-0.3333333, -0.0597158717,-0.8805682564, -0.0597158717, -0.7974269853,  0.5948539707,-0.7974269853};
double wgaus[]={0.45,0.2647883055, 0.2647883055, 0.2647883055, 0.251878361, 0.251878361, 0.251878361};