public void takeOff(DifferentialMeanings d, IntegralMeanings S, Forses F, Moments M, Controls del, Statics s)
        {
            d.Vx   = (1 / s.m) * F.x;
            d.BETA = (1 / (s.m * S.Vx)) * F.z + S.Wy;
            d.Wy   = 1 / s.Iyy * M.y;
            d.X    = S.Vx * cos(S.BETA);
            d.Z    = -S.Vx * sin(S.BETA);

            S.Y = s.Cyo * S.q * s.S;
            S.Q = s.Cxo * S.q * s.S;

            F.x = 20 * del.Tr - S.Q - (s.G - S.Y) * (s.froll + s.fslip * (s.InhibitionL + s.InhibitionR));
            F.z = F.zBETA * S.BETA + F.zDELn * del.n + F.zt;

            F.zBETA = s.Cz_BETA * s.S * S.q;
            F.zDELn = s.Cz_del_n * s.S * S.q;
            F.zt    = (s.G - S.Y) * s.Czk_BETA * (S.BETA + s.fW * del.fW);

            M.y = M.yz + M.yWy + M.yp + M.yt;

            M.yz  = s.my_BETA * S.q * s.S * s.L * S.BETA + s.my_DELn * del.n;
            M.yWy = s.my_Wy * S.q * s.S * s.L * d.Wy;
            M.yp  = 0; // = (P1 - P2) * s.Zengine
            M.yt  = (s.fW * s.XfrontWheel * del.fW + s.fslip * (s.G - S.Y) * (s.InhibitionL + s.InhibitionR) * (s.Zchassis / 2));

            S.q = (s.PO_h * S.Vx * S.Vx / 2);

            d.Xg = S.Vx * cos(S.BETA);
        }
        /// <summary>
        /// Set initial conditions
        /// </summary>
        public void SetInitialConditions()
        {
            d = new DifferentialMeanings();
            S = new IntegralMeanings();
            F = new Forses();
            M = new Moments();
            C = new Coefficient();
            c = new LinearizedCoefficient();
            N = new Overload();
            del = new Controls();
            s = new Statics();

            SetInitialValues();

            if (stForm == null)
            { stForm = new StaticsForm(this, s); }
            else if (stForm != null)
            {
                stForm.Close();
                stForm = new StaticsForm(this, s);
            }

            if (fpForm == null)
            { fpForm = new FlyParamForm(this); }
            else if (fpForm != null)
            {
                fpForm.Close();
                fpForm = new FlyParamForm(this);
            }

            if (glForm == null)
            { glForm = new GraphicsListForm(this); }
            else if (glForm != null)
            {
                glForm.Close();
                glForm = new GraphicsListForm(this);
            }

            if (dgForm == null)
            { dgForm = new DataGridForm(this); }
            else if (dgForm != null)
            {
                dgForm.Close();
                dgForm = new DataGridForm(this);
            }

            if (avForm == null)
            { avForm = new AviatingForm(this); }
            else if (avForm != null)
            {
                avForm.Close();
                avForm = new AviatingForm(this);
            }

            if (pvForm == null)
            { pvForm = new PlaneVisualisationForm(this.toolStripButton6_Visualisation, S); }
            else if (pvForm != null)
            {
                pvForm.Close();
                pvForm = new PlaneVisualisationForm(this.toolStripButton6_Visualisation, S);
            }

            if (auForm == null)
            { auForm = new AutopilotForm(this); }
            else if (auForm != null)
            {
                auForm.Close();
                auForm = new AutopilotForm(this);
            }

            if (scrForm == null)
            { scrForm = new ScriptingAutopilotForm(this); }
            else if (scrForm != null)
            {
                scrForm.Close();
                scrForm = new ScriptingAutopilotForm(this);
            }
        }
        public void takeOff(DifferentialMeanings d, IntegralMeanings S, Forses F, Moments M, Controls del, Statics s)
        {
            d.Vx = (1 / s.m) * F.x;
            d.BETA = (1 / (s.m * S.Vx)) * F.z + S.Wy;
            d.Wy = 1 / s.Iyy * M.y;
            d.X = S.Vx * cos(S.BETA);
            d.Z = -S.Vx * sin(S.BETA);

            S.Y = s.Cyo * S.q * s.S;
            S.Q = s.Cxo * S.q * s.S;

            F.x = 20 * del.Tr - S.Q - (s.G - S.Y) * (s.froll + s.fslip * (s.InhibitionL + s.InhibitionR));
            F.z = F.zBETA * S.BETA + F.zDELn * del.n + F.zt;

            F.zBETA = s.Cz_BETA * s.S * S.q;
            F.zDELn = s.Cz_del_n * s.S * S.q;
            F.zt = (s.G - S.Y) * s.Czk_BETA * (S.BETA + s.fW * del.fW);

            M.y = M.yz + M.yWy + M.yp + M.yt;

            M.yz = s.my_BETA * S.q * s.S * s.L * S.BETA + s.my_DELn * del.n;
            M.yWy = s.my_Wy * S.q * s.S * s.L * d.Wy;
            M.yp = 0; // = (P1 - P2) * s.Zengine
            M.yt = (s.fW * s.XfrontWheel * del.fW + s.fslip * (s.G - S.Y) * (s.InhibitionL + s.InhibitionR) * (s.Zchassis / 2));

            S.q = (s.PO_h * S.Vx * S.Vx / 2);

            d.Xg = S.Vx * cos(S.BETA);
        }