mainSimul function

simulate many times the first model and calculate the risk constant