циклы — вводный C ++ и метод Монте-Карло

Так что я только знакомлюсь с программированием, и я должен сказать, что это единственная самая полезная, но разочаровывающая вещь, которую я когда-либо делал. Я беру проекты возрастающей сложности, последний из которых связан с использованием метода Монте-Карло и множеством циклов. Ниже приведен код, который уже завершен:

#include <iostream>
#include <cstdlib>
#include <cmath>
#include <string>
#include <ctime>

using namespace std;int main ()
{
srand (time(0));
string operation;
cout << "Using the letters 'o', or 'q', please indicate if you would like to simulate once, or quit the program: " << endl;
cin >> operation;
string reservoir_name; // Creating variables for reservoir
double reservoir_capacity;
double outflow;
double inflow_min;
double inflow_max;

if (operation == "q")
{
cout << "Exiting program." << endl;
system ("pause");
return 0;
}

while (operation == "o") // Choose one or multiple simulations.
{

string reservoir_name; // Creating variables for reservoir function
double reservoir_capacity;

double inflow_min = 0;
double inflow_max = 0;
double inflow_range = inflow_min + inflow_max;
double inflow_difference = inflow_max - inflow_min;
double inflow_threshold = .9 * inflow_range/2; // Math for acceptable flow threshold.cout << "What is the name of the reservoir?" << endl;
cin.ignore ();
getline (cin,reservoir_name); // Grab whole string for reservoir name.
cout << "What is the capacity of the reservoir in MAF (Millions of Acre Feet)?" << endl;
cin >> reservoir_capacity;
cout << "What is the minimum inflow?" << endl;
cin >> inflow_min;
cout << "What is the maximum inflow?" << endl;
cin >> inflow_max;
cout << "What is the required outflow?" << endl;
cin >> outflow;
inflow_range = inflow_min + inflow_max;
inflow_threshold = .9 * inflow_range/2;
cin.ignore ();

if (outflow > inflow_threshold) // Check for unacceptable outflow levels.
{
cout << "Warning! The outflow is over 90% of the average inflow. Simulation aborted. Returning to main menu." << endl;
}
else
{
const int number_simulations = 10;
double fill_level = 0;
int years = 1;
cout << "Running simulation." << endl;
for (int i = 1; i < number_simulations; i++) // Each year
{

for (years; fill_level < reservoir_capacity; years++ )
{
double r = rand() * 1.0 / RAND_MAX;
double x = inflow_min + inflow_range * r;// SHOULD be between minimum inflow and maximum inflow.
if (fill_level < 0)
{
fill_level = 0;
}
}   // Simulate the change of water level.
cout << years << endl;
}

}
cout << "What would you like to do now?" << endl; // Saving for later. The menu re-prompt message and code.
cout << "Using the letters 'o', or 'q', please indicate if you would like to simulate once, or quit the program: " << endl;
cin >> operation;
}system ("pause");
return 0;
}

Итак, я полагаю, что мой главный вопрос заключается в том, что я сталкиваюсь со стеной, касающейся настройки циклов for под «Запуск симуляции», где мне нужно настроить первый цикл for для запуска внутреннего цикла for 10 раз с каждой из этих 10 итераций. внутреннего цикла for, предлагающего случайные числа для диапазона приемлемых результатов запроса для случайного значения. Мне сказали, что идея состоит в том, чтобы использовать метод Монте-Карло, т.е.

double r = rand() * 1.0 / RAND_MAX;
double x = inflow_min + inflow_range * r;// SHOULD be between minimum inflow and maximum inflow.

поэтому программа создаст случайное значение для притока. Идея состоит в том, что внутренний цикл for будет продолжать работать до тех пор, пока fill_level резервуара, который начинается с 0, не достигнетservoir_capacity. Процесс имитации того, сколько лет (каждая итерация внутреннего цикла for, представляющего год) должен повторяться 10 раз родительским циклом для имитации fill_level для цикла.

Когда я пытаюсь запустить программу, как вы видите ее здесь, она будет проходить до самого «Запуск симуляции», а затем не будет дальше. Кто-то более опытный, чем я, понимает, что я говорю, и знает, что происходит?