vector<double> minPar;
    vector<double> maxPar;
    minPar.push_back(0.1); maxPar.push_back(7.);
    minPar.push_back(1.); maxPar.push_back(7.);
    minPar.push_back(0.001); maxPar.push_back(5.);
    minPar.push_back(0.001); maxPar.push_back(5.);
    minPar.push_back(0.0001); maxPar.push_back(7.);
    minPar.push_back(1e-10); maxPar.push_back(0.1);

    minPar.push_back(0.001); maxPar.push_back(5.);
    minPar.push_back(5000); maxPar.push_back(50000);
    minPar.push_back(0.001); maxPar.push_back(50.);
    minPar.push_back(-10.); maxPar.push_back(-0.00001);
    vector<double> pars = PS.params;
    for (int i = 0; i < 1; i++)
    {
        for (int j = 0; j < 10; j++)
        {
            ofstream grfile;
            grfile.open("pargraph\\par" + to_string(i * 12 + j) + ".dt");
            for (int k = 0; k < 100; k++)
            {

                pars[i * 12 + j] = minPar[j] + (double)k*(maxPar[j] - minPar[j]) / 100.0;
                PS.calcProps(pars);
                double fmin = PS.getFunMin();
                cout << (i * 12 + j) * 100 + k << "\t" << fixed << setprecision(15) << pars[i * 12 + j] << "\t" << fmin << endl;

                grfile << pars[i * 12 + j] << ", " << fixed << setprecision(15) << fmin << endl;
            }
            grfile.close();
            pars[i * 12 + j] = PS.params[i * 12 + j];
        }
    }