Функция для системы нелинейных уравнений методом Ньютона и синтаксическая ошибка - C#

double Newton(double x, double y)
        {
            double f1, f2, x_sled, y_sled;
            double df1dx, df1dy, df2dx, df2dy, J;
            do{
            i++;
            //описание производных
            df1dx = cos(x + 1.8);
            df1dy = -2;
            df2dx = 0.9;
            df2dy = -sin(y + 1.6);
            //Якобиан
            J = df1dx * df2dy - df1dy * df2dx;
            //вычисление очередного значения х
            x_sled = x - ((f1 * df2dy - f2 * df1dy) / J);
            y_sled = y + ((f1 * df2dx - f2 * df2dx) / J);
            if (Math.Abs(x_sled - x) < eps && Math.Abs(y_sled - y) < eps)
                goto M;
            x = x_sled;
            y = y_sled;
            }
            while (1);
            M: return(x_sled,y_sled);
        }