dx = [j11(x_old),j12(x_old),j13(x_old);j21(x_old),j22(x_old),j23(x_old);j31(x_old),j32(x_old),j33(x_old)]\[f1(x_old);f2(x_old);f3(x_old)]; ...