//x
        J_sec[0][0] = fx1;
        J_sec[1][0] = fx2;
        J_sec[2][0] = fx3;

        for (j = 0; j < 3; j++) {
            for (n = 1; n < 3; n++) {
                J_sec[j][n] = J[j][n];
                printf("k = %d, n = %d\nJ_sec[%d][%d] = J[%d][%d]\n", k,n,j,n,j,n);
            }
        }
        w_x = wyzn(J_sec);

        //y

        J_sec[0][1] = fx1;
        J_sec[1][1] = fx2;
        J_sec[2][1] = fx3;

        for (j = 0; j < 3; j++) {
            for (n = 0; n < 3; n++) {
                if (n != 1) {
                    J_sec[j][n] = J[j][n];
                    printf("k = %d, n = %d\nJ_sec[%d][%d] = J[%d][%d]\n", k,n,j,n,j,n);
                }
            }
        }
        w_y = wyzn(J_sec);

        //z
        J_sec[0][2] = fx1;
        J_sec[1][2] = fx2;
        J_sec[2][2] = fx3;

        for (j = 0; j < 3; j++) {
            for (n = 0; n < 2; n++) {
                if (n != 1) {
                    J_sec[j][n] = J[j][n];
                    printf("k = %d, n = %d\nJ_sec[%d][%d] = J[%d][%d]\n", k,n,j,n,j,n);
                }
            }
        }
        w_z = wyzn(J_sec);