Exemplo n.º 1
0
 static double[][] Outer(int g, int r, int x)
 {
     Rk4 solver = new Rk4(
       (double t, double[] V, double[] res) =>
       {
     // Console.WriteLine(t);
     res[0] = rate(t) * V[0] - GmFemale(x+t) * (Middle(g,r,x+t,t) - V[0]);
       },
       1);
     double[][] reserves = solver.SolveAll(120-x, 0, steps, new double[] { 0.0 });
     return reserves;
 }