private static void ROT(double g ,double h ,double s ,double tau, double[,] a ,int i ,int j ,int k ,int l) { g=a[i,j]; h=a[k,l]; a[i,j]=g-s*(h+g*tau); a[k,l]=h+s*(g-h*tau); }