Нахождение всех частных производных полинома от нескольких переменных в заданной точке по схеме Горнера

public class Basic_1_2_14 {
 public static void main(String[] args) {
        //Gurner scheme
        int m = 2,n = 3;
        double[] x = new double[m];
        int[][] a = new int[m][n];
        handmadeFill(a,x);
        System.out.print("a = ");print(a, x);
        gurner(a,x);
    }
 
    private static void handmadeFill(int[][] a, double[] x) {
        double[] xt = {2,3};
        int[][] at = {
                {2,5,3},
                {3,1,2}
             };
        for (int i = 0;  i<x.length; i++) {
            x[i] = xt[i];
        }
        for (int i = 0;  i<a.length; i++) {
            for (int j = 0;  j<a[0].length; j++) {
                a[i][j] = at[i][j];
            }
        }
    }
    private static void gurner(int[][] a, double[] x){
        int m = a.length;
        int n = a[0].length-1;
        double[] d = new double[m];
        double[] dy = new double[m];
        double ma = 1,y = 1;
 
        for (int i = 0; i < m; i++) {
            ma *= x[i];
            y *= a[i][n];
        }
 
        int j = 0;
 
        while (j != n) {
            j++;
            int ca = 1;
            for (int i = 0; i < m; i++) {
                ca *= a[i][n-j];
            }
            y = y * ma + ca;
        }
 
        System.out.println("value = " + y);
    }
    private static void print(int[][] a,double[] x) {
        System.out.println();
        for (int i = 0;  i<a.length; i++) {
            for (int j = 0;  j<a[0].length; j++) {
                System.out.print((j!=0?" + ":"") + a[i][j] + "*" + x[i] + "^"+ (a[0].length-j-1)+"");
            }
            System.out.println();
        }
        System.out.println();
    } 
}

 private static void gurner(int[] a, double x){
        int n = a.length - 1, k = 0;
        double y = (double) a[n];
        double dy = 0; //this is the derevative
        //invariant: 0 <= k <= n,
        // y = a[n]*(x in degree k)+...+a(n-1)*(x in degree k-1)+...+
        //     + a(n-k]*(x in degree 0)
        while ( k != n) {
            k++;
            dy = dy*x + y;
            y = y*x + a[n-k];
        }
        System.out.println("Value in the point f(x) = " + y + " " + x);
        System.out.println("Value in the point f'(x)  = " + dy + " " + x);
  }