using System;

​

public class Program

{

    public static void menu()

    {

        Console.WriteLine("Menu \n1. Transpose Matriks \n2. Invers Matriks \n3. Determinan Matriks");

    }

​

    public static void Main()

    {

        int i, j;

        int baris = 3;

        int kolom = 3;

        int[, ] mat = new int[3, 3];

        for (i = 0; i < baris; i++)

        {

            for (j = 0; j < kolom; j++)

            {

                Console.Write("Kolom " + (i + 1) + " Baris " + (j + 1) + "=");

                mat[i, j] = Convert.ToInt32(Console.ReadLine());

            }

        }

​

        for (i = 0; i < baris; i++) // Matriksnya

        {

            for (j = 0; j < kolom; j++)

            {

                Console.Write(mat[i, j] + " ");

            }

​

            Console.WriteLine(" ");

        }

​

        {

            menu();

            Console.Write("Input = ");

            int pilih = Convert.ToInt32(Console.ReadLine());

            {

                if (pilih == 1) // transpose matrix

                {

                    for (i = 0; i < baris; i++)

                    {

                        for (j = 0; j < kolom; j++)

                        {

                            Console.Write(mat[j, i] + " ");

                        }

​

                        Console.WriteLine(" ");

                    }

​

                    menu();

                    Console.Write("Input = ");

                    pilih = Convert.ToInt32(Console.ReadLine());

                }

                else if (pilih == 3) // determinan matriks

                {

                    double det = (((mat[0, 0] * mat[1, 1] * mat[2, 2]) + (mat[0, 1] * mat[1, 2] * mat[2, 0]) + (mat[0, 2] * mat[1, 0] * mat[2, 1])) - ((mat[2, 0] * mat[1, 1] * mat[0, 2]) + (mat[2, 1] * mat[1, 2] * mat[0, 0]) + (mat[2, 2] * mat[1, 0] * mat[0, 1])));

                    Console.WriteLine("Determinan Matriks= " + det);

                    menu();

                    Console.Write("Input = ");

                    pilih = Convert.ToInt32(Console.ReadLine());

                }

                else if (pilih == 2) // invers matriks

                {

                    double det = (((mat[0, 0] * mat[1, 1] * mat[2, 2]) + (mat[0, 1] * mat[1, 2] * mat[2, 0]) + (mat[0, 2] * mat[1, 0] * mat[2, 1])) - ((mat[2, 0] * mat[1, 1] * mat[0, 2]) + (mat[2, 1] * mat[1, 2] * mat[0, 0]) + (mat[2, 2] * mat[1, 0] * mat[0, 1])));

                    double bagian = 1 / det;

                    double[, ] adj = new double[3, 3]; // sdh termasuk ditranspose

                    //kolom 1

                    adj[0, 0] = (mat[1, 1] * mat[2, 2] - mat[2, 1] * mat[1, 2]) * bagian;

                    adj[1, 0] = -(mat[1, 0] * mat[2, 2] - mat[2, 0] * mat[1, 2]) * bagian;

                    adj[2, 0] = (mat[1, 0] * mat[2, 1] - mat[2, 0] * mat[1, 1]) * bagian;

                    //kolom 2

                    adj[0, 1] = -(mat[0, 1] * mat[2, 2] - mat[2, 1] * mat[0, 2]) * bagian;

                    adj[1, 1] = (mat[0, 0] * mat[2, 2] - mat[0, 2] * mat[2, 0]) * bagian;

                    adj[2, 1] = -(mat[0, 0] * mat[2, 1] - mat[0, 1] * mat[2, 0]) * bagian;

                    //kolom 3

                    adj[0, 2] = (mat[0, 1] * mat[1, 2] - mat[0, 2] * mat[1, 1]) * bagian;

                    adj[1, 2] = -(mat[0, 0] * mat[1, 2] - mat[1, 0] * mat[0, 2]) * bagian;

                    adj[2, 2] = (mat[0, 0] * mat[1, 1] - mat[1, 0] * mat[0, 1]) * bagian;

                    Console.WriteLine(adj[0, 0] + " " + adj[0, 1] + " " + adj[0, 2]);

                    Console.WriteLine(adj[1, 0] + " " + adj[1, 1] + " " + adj[1, 2]);

                    Console.WriteLine(adj[2, 0] + " " + adj[2, 1] + " " + adj[2, 2]);

                    menu();

                    Console.Write("Input = ");

                    pilih = Convert.ToInt32(Console.ReadLine());

                }

            }

        }

    }

}