Протабулировать функцию у на заданном отрезке

using System;
 
 
namespace ConsoleApp5
{
    class Program
    {
        static double func (double x)
        {
           return 3*Math.Cos(5*x) + 2*Math.Sin(x);
        }
 
        static void Printarr(double[] arr)
        {
            for (int j = 0; j < arr.Length; j++)
                Console.Write(" {0:0.0000}", arr[j]);
            Console.WriteLine();    
        }
 
        static void Main(string[] args)
        {
           
            double[] a = new double[25];
            double[] b = new double[25];
            double[] wrk = new double[25];
            double h = Math.PI / 10;
            int g = 0;
            double x1 = - Math.PI;
            Console.WriteLine("Табуляция функции:");
            Console.WriteLine("    X     |     Y  ");
            while (x1 <= Math.PI)
            {
                wrk[g] = func(x1);
                Console.WriteLine(" {0,8:0.0000} | {1,8:0.0000}", x1, wrk[g]);
               
                g = g + 1;
                x1 = x1 + h;
          
            }
            double maxelem = double.MinValue;
            double minelem = double.MaxValue;
            int kmax = 0;
            int kmin = 0;
            for (int k = 0; k < g; k++)
            {
                if (wrk[k] < minelem)
                {
                    minelem = wrk[k];
                    kmin = k;
                }
 
                if (wrk[k] > maxelem)
                {
                    maxelem = wrk[k];
                    kmax = k;
                }
            }
            Console.WriteLine("Элемент с максимальным значением {0:0.0000}", maxelem);
            Console.WriteLine("Элемент с минимальным значением {0:0.0000}", minelem);
            for (int j = 0; j < kmin; j++)
                a[j] = wrk[j];
            int j1 = 0;
            for (int j = 0; j < kmax; j++)
                if (wrk[j] < 0)
                {
                    b[j1] = wrk[j];
                    j1 = j1 + 1;
                }
            Array.Resize(ref a, kmin);
            Array.Resize(ref b, j1);
            Console.WriteLine("Массив A:");
            Printarr(a);
            Console.WriteLine("Массив B:");
            Printarr(b);
            Array.Reverse(a);
            Console.WriteLine("Измененный массив A:");
            Printarr(a);
            Array.Sort(b);
            Array.Reverse(b);
            Console.WriteLine("Отсортированный массив B:");
            Printarr(b);
            int n = Math.Min(a.Length, b.Length);
            Console.WriteLine("Наименьшая размерность массива: {0}",n);
            double[,] e = new double[n, n];
 
            for (int i = 0; i < n; i++)
                for (int j = 0; j < n; j++)
                {
                    if (i == j)
                        e[i, j] = 1;
                    else
                        e[i, j] = a[i] * b[j] - a[j] * b[i];  // Формула для заполнения боковых диагоналей матрицы.
 
                }
            Console.WriteLine("Матрица:");
            for (int i = 0; i < n; i++)
            {
                for (int j = 0; j < n; j++)
                    Console.Write(" {0,7:0.0000}", e[i,j]);
                Console.WriteLine();
            }
            int cou = 0;
            for (int i = 0; i < n; i++)
                for (int j = 0; j < n; j++)
                    if (e[i,j] < 0)
                        cou = cou + 1;
            int Tret = Convert.ToInt32(n * n / 3);
            Console.WriteLine("Количество отрицательных элементов {0} ",cou);
            Console.WriteLine("Одна третья элементов {0} (округлено до целого) ", Tret);
            if (cou > Tret)
            {
                Console.WriteLine("Количество отрицательных элементов матрицы превышает треть количества элементов");
                double tmp;
                for (int i = 0; i < n; i++)
                {
                    for (int j = 0; j < i; j++)
                    {
                        tmp = e[i,j];
                        e[i,j] = e[j,i];
                        e[j,i] = tmp;
                    }
                }
                Console.WriteLine("Транспонированная матрица:");
                for (int i = 0; i < n; i++)
                {
                    for (int j = 0; j < n; j++)
                        Console.Write(" {0,7:0.0000}", e[i, j]);
                    Console.WriteLine();
                }
            }
            else
                Console.WriteLine("Количество отрицательных элементов меньше  одной трети элементов матрицы");
 
 
            Console.ReadLine();
        }
    }
}