using System;

public class Program

{

    public static double D(double x)

    {

        return (x - (int)x);

    }

    public static double[] LPTAU(int[,] NR, int i,int n, double[] q)

    {

        int a = i;

        double m = 1+(int)(Math.Log(a)/0.693147);

        for (int j=1;j<=n;j++)

        {

            double s = 0;

            for (int k=1;k<=m;k++)

            {

                double ns = 0;

                for (int l=k;l<=m;l++)

                {

                 int b = NR[j-1,l-1];

                    ns = ns + (int)(2*D(Math.Pow(a/2,l)))*(int)(2*D(Math.Pow(b/2,l+1-k)));

                }

                s = s + Math.Pow(D(0.5*ns)/2,k-1);

            }

            q[j]=s;

        }

        return q;

    }

    public static void Main()

    {

        int[,] NR = { 

{1, 1, 1, 1, 1, 1, 1, 1, 1, 1}, 

{1, 3, 5, 15, 17, 51, 85, 255, 257, 771}, 

{1, 1, 7, 11, 13, 61, 67, 79, 465, 721}, 

{1, 3, 7, 5, 7, 43, 49, 147, 439, 1013}, 

{1, 1, 5, 3, 15, 51, 125, 141, 177, 759}, 

{1, 3, 1, 1, 9, 59, 25, 89, 321, 835}, 

{1, 1, 3, 7, 31, 47, 109, 173, 181, 949}, 

{1, 3, 3, 9, 9, 57, 43, 43, 225, 113} 

};

        int n=2;

        double[] q = new double[n+1];

        for (int i=0;i<100;i++)

        {

            double[] qq = new double[n];

            q = LPTAU(NR,i,n,q);

​

            Console.WriteLine(q[1]+" "+q[2]);

        }

    }

}