ошыбка в параметре з ref, наверное...

[Natok]

вот есть у меня реализацыя метода Рунге-Кутта для решения системи диф.уравнений 1 порядка.

и есть система:
y[1]'=F(i, x, y);
y[2]' =F(i, x, y);

вот так выглядят мои начальные условия:
(0, y[1], y[2])

n=2 в моем случае.

using System;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Drawing;
using System.Linq;
using System.Text;
using System.Windows.Forms;

namespace WindowsFormsApplication9
{
    public partial class Form1 : Form
    {
        public Form1()
        {
            InitializeComponent();
        }
        static double f(int i,
        double x,
        double[] y)
        {
            if (i == 1)
            {
                double  k1, k2;
                k1 = 0.2;
                k2 = 0.5;

                return k1*y[1]-k2*y[2]*y[1];
            }

            else if (i == 2)
            {
                double k3, k4;

                k3 = 0.4;
                k4 = 0.3;
                return -k3 * y[2] + k4 * y[1] * y[2];
            }
            else return 0;


        }/*************************************************************************
    Метод Рунге-Кутта четвертого порядка для решения
    системы дифферециальных уравнений.

    
    Алгоритм решает систему диффуров y[i]'=F(i,x,y) для i=1..n
    методом Рунге-Кутта 4 порядка.

    Начальная точка имеет кординаты (x,y[1], ..., y[n])

    До конечной точки мы добираемся через n промежуточных
    с постоянным шагом h=(x1-x)/m

    Результат помещается в переменную y
    *************************************************************************/
        public static void solvesystemrungekutta(double x,
        double x1,
        int m,
        int n,
        ref double[] y)
        {
            //double h = 0;
            int i = 0;

            for (i = 0; i <= m - 1; i++)
            {
                solvesystemrungekuttastep(x + i * (x1 - x) / m, (x1 - x) / m, n, ref y);
            }
        }
/*************************************************************************
    Один шаг метода Рунге-Кутта четвертого порядка для решения
    системы дифферециальных уравнений.

    
    Алгоритм совершает один шаг метода для системы
    диффуров y[i]'=F(i,x,y) для i=1..n

    Начальная точка имеет кординаты (x,y[1], ..., y[n])

    После выполнения алгоритма в переменной y содержится состояние
    системы в точке x+h
    *************************************************************************/

         public static void solvesystemrungekuttastep(double x,
        double h,
        int n,
        ref double[] y)
    {
        int i = 0;
        double[] yt = new double[0];
        double[] k1 = new double[0];
        double[] k2 = new double[0];
        double[] k3 = new double[0];
        double[] k4 = new double[0];

        yt = new double[n+1];
        k1 = new double[n+1];
        k2 = new double[n+1];
        k3 = new double[n+1];
        k4 = new double[n+1];
        for(i=1; i<=n; i++)
        {
            k1[i] = h*f(i, x, y);
        }
        for(i=1; i<=n; i++)
        {
            yt[i] = y[i]+0.5*k1[i];
        }
        for(i=1; i<=n; i++)
        {
            k2[i] = h*f(i, x+h*0.5, yt);
        }
        for(i=1; i<=n; i++)
        {
            yt[i] = y[i]+0.5*k2[i];
        }
        for(i=1; i<=n; i++)
        {
            k3[i] = h*f(i, x+h*0.5, yt);
        }
        for(i=1; i<=n; i++)
        {
            yt[i] = y[i]+k3[i];
        }
        for(i=1; i<=n; i++)
        {
            k4[i] = h*f(i, x+h, yt);
        }
        for(i=1; i<=n; i++)
        {
            y[i] = y[i]+(k1[i]+2.0*k2[i]+2.0*k3[i]+k4[i])/6;
        }
    }

        private void button1_Click(object sender, EventArgs e)
        {
             double[] y1;
            double[] yrez;
            y1=new double [3];

            y1[1]=20;
            y1[2]=10;
            
            solvesystemrungekutta(1, 5, 20 , 2,   ref y1);//<---------вызывает ошыбку. 
        
        }
}



    }

[Natok]

ага, понятно ,вызывать надо так:

solvesystemrungekutta(1, 5, 20, 2, ref y1);

написано очень класно здесь:

http://msdn.microsoft.com/en-us/library/14akc2c7.aspx

но еще есть проблеми с индексами...выходит за границы в исчислениях

[Natok]

поправила, уже работает...