Написать программу циклами

using System;
using System.IO;
class Data
{
    long z = 1000000;
    double[,] Yp;
    double[,] y;
    double[] time;
    double[] m = new double[4];
    double[] n = new double[4];
    double[] l = new double[4];
    double[] p = new double[4];
    double dt;
    double k1;
    double KU;
    double S;
    double T1;
    double k2;
    double TK;
    double C;
    double TM;
    double II;
    long k;
    int Qin;
    double a;
    double b;
    double y1;
    double y2;
    double y3;
    double y4;
    double e;
    double d;
    public Data()
    {
        Yp = new double[4, z];
        y = new double[4, z];
        time = new double[z];
        KU = 1;
        TM = 0.02;
        TK = 0.01;
        T1 = 0.002;
        C = 20;
        II = 400;
        k1 = 5;
        k2 = 2;
        S = 60;
        Qin = 1;
        dt = 0.001;
        e = 0.0001;
        a = 0;
        b = 0.8;
        time[0] = a;
    }
    public double F1(double t, double y1, double y2, double y3, double y4)
    {
        return (k1 * KU * S * (Qin - y2) - y1) / T1;
    }
    public double F2(double t, double y1, double y2, double y3, double y4)
    {
        return y4 / II;
    }
    public double F3(double t, double y1, double y2, double y3, double y4)
    {
        return (k2 * y1 - y3) / TK;
    }
    public double F4(double t, double y1, double y2, double y3, double y4)
    {
        return (C * y3 - y4) / TM;
    }
    public void Pohidni()
    {
        Yp[0, k] = F1(time[k], y[0, k], y[1, k], y[2, k], y[3, k]);
        Yp[1, k] = F2(time[k], y[0, k], y[1, k], y[2, k], y[3, k]);
        Yp[2, k] = F3(time[k], y[0, k], y[1, k], y[2, k], y[3, k]);
        Yp[3, k] = F4(time[k], y[0, k], y[1, k], y[2, k], y[3, k]);
    }
    public void Znach()
    {
        y1 = y[0, k] + dt * F1(time[k], y[0, k], y[1, k], y[2, k], y[3, k]) / 2.0;
        y2 = y[1, k] + dt * F2(time[k], y[0, k], y[1, k], y[2, k], y[3, k]) / 2.0;
        y3 = y[2, k] + dt * F3(time[k], y[0, k], y[1, k], y[2, k], y[3, k]) / 2.0;
        y4 = y[3, k] + dt * F4(time[k], y[0, k], y[1, k], y[2, k], y[3, k]) / 2.0;
        y[0, k + 1] = y[0, k] + dt * F1(time[k] + dt / 2.0, y1, y2, y3, y4);
        y[1, k + 1] = y[1, k] + dt * F2(time[k] + dt / 2.0, y1, y2, y3, y4);
        y[2, k + 1] = y[2, k] + dt * F3(time[k] + dt / 2.0, y1, y2, y3, y4);
        y[3, k + 1] = y[3, k] + dt * F4(time[k] + dt / 2.0, y1, y2, y3, y4);
    }
    public void Provirka()
    {
        do
        {
            Znach();
            d = l[0] / 150.0 - l[2] * 3.0 / 100.0 + l[3] * 16.0 / 75.0 + l[4] / 20.0 - l[5] * 6.0 / 25.0;
            if (Math.Abs(d) > e)
                dt = dt / 2.0;
            if (Math.Abs(d) < e / 30.0)
                dt = 2.0 * dt;
        } while (Math.Abs(d) > e || Math.Abs(d) < e / 30.0);
    }
    public void Prod()
    {
        do
        {
            Pohidni();
            Znach();
            time[k + 1] = time[k] + dt;
            k++;
        } while (time[k] < b);
        Pohidni();
    }
    public void vuvid()
    {
        Random rand = new Random(); 
        StreamWriter log_out;
        log_out = new StreamWriter("logfile.txt");
        Console.SetOut(log_out);
        for (int j = 0; j < k; j += rand.Next(1, 5)) Console.WriteLine("{0:0.####}\t{1:0.###}", time[j], y[1, j]);
        log_out.Close();
    }
    public void vuv()
    {
        Console.WriteLine("t\ty1\ty1'");
        for (int j = 0; j < k; j += 2)
        {
            Console.WriteLine("{0:0.####}\t{1:0.###}", time[j], y[1, j]);
        }
    }
}
class Program
{
    static void Main(string[] args)
    {
        Data d = new Data();
        d.Prod();
        d.vuv();
        d.vuvid();
        Console.ReadLine();
    }
}