Lake of Fire

Kaotic Klique – What It Is

BILL $ABER - MONSTER*

Wasted Wednesday

Keep Believing

intro track no. 1 from the upcoming album.
Credits to "Jhonny Charles" for the instrumental.

Lux - One (Abusive)

Care

MARI - Broken Heart (Prod. DevilPrado)

Gone for the Moment

Gone for the Moment (Produced by Chunksta D)

Vice city (feat 29Dimes)

In peace

Whats's up bro?
My name is Kirill, i'm beginner beatmaker and rap artist. Make
using System;

class Program
{
    // Напишем класс вектора из 3 координат
    class Vector
    {
        public double a, b, c;
        public Vector()
        {
            
            a = 5.1;
            b = 5.1;
            c = 2.9;
        }
        public Vector(double a, double b, double c)
        {
            this.a = a;
            this.b = b;
            this.c = c;
        }
        public static Vector operator +(Vector v1, Vector v2)
        {
            return new Vector { a = v1.a + v2.a, b = v1.b + v2.b, c = v1.c + v2.c };
        }
        public static Vector operator -(Vector v1, Vector v2)
        {
            return new Vector { a = v1.a - v2.a, b = v1.b - v2.b, c = v1.c - v2.c };
        }
        public static Vector operator *(double number, Vector v1)
        {
            return new Vector { a = v1.a * number, b = v1.b * number, c = v1.c * number };
        }

        // Перегружаем оператор <
        // Принимает вектор и число, если каждая координата вектора, взятая по модулю, меньше числа - 
        // то возвращает true, иначе false
        public static bool operator <(Vector v1, double number)
        {
            return Math.Abs(v1.a) < number && Math.Abs(v1.b) < number && Math.Abs(v1.c) < number;
        }

        // Перегружаем оператор >
        // Принимает вектор и число, если каждая координата вектора, взятая по модулю, больше числа - 
        // то возвращает true, иначе false
        public static bool operator >(Vector v1, double number)
        {
            return Math.Abs(v1.a) > number && Math.Abs(v1.b) > number && Math.Abs(v1.c) > number;
        }
    }
    static Vector MultMatrixOnVector(Vector v1)
    {
        Vector NewVector = new Vector
        {
            a = A[0, 0] * v1.a + A[0, 1] * v1.b + A[0, 2] * v1.c,
            b = A[1, 0] * v1.a + A[1, 1] * v1.b + A[1, 2] * v1.c,
            c = A[2, 0] * v1.a + A[2, 1] * v1.b + A[2, 2] * v1.c
        };
        return NewVector;
    }
    static double ScalarMult(Vector v1, Vector v2)
    {
        return v1.a * v2.a + v1.b * v2.b + v1.c * v2.c;
    }

    // Матрица A
    static double[,] A;

    // Вектор b
    static Vector b = new Vector();

    static void Main(string[] args)
    {
        // Заносим в матрицу нужные данные
        A = new double [3, 3]
        {
            { 6,  0.6,  0.7 },
            {  0.6, 3, -0.3},
            {  0.7, -0.3,  4 }
        };

        // Переменная k будет отвечать за максимальное число шагов,
        // чтобы программа не зациклилась, если будет расхождение
        Console.Write("Введите максимальное число шагов для алгоритма: ");
        int k = int.Parse(Console.ReadLine());

        // Введем точность вычисления
        double eps = 0.0001;

        // Начальные данные для алгоритма
        // Так как матрица A диагонально преобладающая, то в качестве X0 берем bi / Aii
        Vector Xk = new Vector(b.a / A[0, 0], b.b / A[1, 1], b.c / A[2, 2]);
        Vector Rk = b - MultMatrixOnVector(Xk);
        Vector Zk = new Vector(Rk.a, Rk.b, Rk.c);

        // В цикле идем до тех пор, пока не достигнем нужной точности (Rk > eps),
        // либо пока не истратим лимит шагов k
        int i = 0;
        for ( i = 0; Rk > eps && i < k; ++i)
        {
            double ak = ScalarMult(Rk, Rk) / ScalarMult(MultMatrixOnVector(Zk), Zk);
            Xk = Xk + ak * Zk;
            Vector Rk1 = new Vector(Rk.a, Rk.b, Rk.c);
            Rk = Rk1 - ak * MultMatrixOnVector(Zk);
            double bk = ScalarMult(Rk, Rk) / ScalarMult(Rk1, Rk1);
            Zk = Rk + bk * Zk;
        }

        // Выводим на экран ответ
        Console.WriteLine("\nX1 = " + Xk.a + "\nX2 = " + Xk.b + "\nX3 = " + Xk.c);
        // Выводим на экран ответ
        Console.WriteLine("\nВектор невязки  \nX= " + Math.Round(Rk.a, 14) + "\nY = " + Math.Round(Rk.b, 7) + "\nZ = " + Math.Round(Rk.c, 7));
    }
}

Vice city (feat 29Dimes)

DevilPrado ✪

Related tracks