Aurora 4x v1.1

int cf, ci, dci, mi, la, bp, rp, mp, br, mr;
    double dt, t, tt, cor, dur, ven, tri, qcor, qdur, qven, qtri, dcf, dmi;

    dt = 1.0 / 365;
    t = 0;
    tt = 0;

    /*Startgebäude*/
    la = 25.0;
    ci = 5000;
    cf = 0;
    mi = 0;
    mr = 10;
    dcf = 0;
    dmi = 0;

    /*Forschung*/
    rp = 100;
    double t_f1 = 5000 / (la *  rp);
    rp = 120;
    double t_f2 = 3000 / (la *  rp);
    rp = 120;
    double t_f3 = 3000 / (la *  rp);

    /*Ressourcen*/
    dur = 1905.0, tri = 2398.0, ven = 4353.0, cor = 2398.0;
    qdur = 706888.0, qtri = 183259.0, qven = 232000.0, qcor = 155259.0;
    double  ratdur = 0.9, rattri = 0.6, ratven = 0.9, ratcor = 0.6;

    /*Bauverhältnisse*/
    double cicf = 0.75;
    double cimi = 1 - cicf;
    dci = 0;
    while (cor >= 0 && dur >= 0 && ven >= 0 && tri >= 0 && cf <= 1000) {

        if (tt >= (t_f1 + t_f2))
            br = 12;
        else
            br = 10;

        if (tt >= (t_f1 + t_f2 + t_f3))
            mr = 12;
        else
            mr = 10;

        bp = (ci - dci) + cf * br;
        mp = (ci - dci) + mi * mr;

        while (dcf <= 1 && dmi <= 1) {
            t += dt;
            dcf += (dt * cicf*bp) / 20;
            dmi += (dt * cimi*bp) / 20;


        }


        cf += (int)dcf;
        mi += (int)dmi;

        dci = mi + cf;

        dur = dur - (int)dmi * 10 - (int)dcf * 10 + mp * t * ratdur;
        tri = tri - (int)dcf * 5 + mp * t * rattri;
        ven = ven - (int)dcf * 5 + mp * t * ratven;
        cor = cor - (int)dmi * 10 + mp * t * ratcor;


        qdur = qdur - mp * t * ratdur;
        qtri = qtri - mp * t * rattri;
        qven = qven - mp * t * ratven;
        qcor = qcor - mp * t * ratcor;



        tt += t;
        t = 0;
        dcf -= (int)dcf;
        dmi -= (int)dmi;

    }

    /*Zeit bis Ressourcen ausgehen*/


    cout.precision(17);

    cout << "Ergebnis " << endl;
    cout << "Time [a] " << fixed << tt << endl;
    cout << "Constructionfactory " << fixed << cf << endl;
    cout << "Mine " << fixed << mi << endl;
    cout << "MiningPoints " << fixed << mp << endl;
    cout << "BuildingPoints " << fixed << bp << endl << endl;

    cout << "Gesamtverbrauch " << endl;
    cout << "Duranium " << fixed << cf * 10 + mi * 10 << endl;
    cout << "Tritanium " << fixed << cf * 5 << endl;
    cout << "Vendarite " << fixed << cf * 5 << endl;
    cout << "Corundium " << fixed << mi * 5 << endl << endl;

    cout << "Stockpile " << endl;
    cout << "Duranium " << fixed << dur << endl;
    cout << "Tritanium " << fixed << tri << endl;
    cout << "Vendarite " << fixed << ven << endl;
    cout << "Corundium " << fixed << cor << endl << endl;

    cout << "Ressourcen " << endl;
    cout << "Quantitiy Duranium " << fixed << qdur << endl;
    cout << "Quantitiy Tritanium " << fixed << qtri << endl;
    cout << "Quantitiy Vendarite " << fixed << qven << endl;
    cout << "Quantitiy Corundium " << fixed << qcor << endl;