format long e x = [1 2 0.03;5 6 0.07]; y = [4;8]; dx1 = [1 0 0.03;5 0 -0.07]; dx2 = [0 2 0.03;0 6 0.07]; j = -dx1/dx2 j.*y