Dear Ray, Thank you. I want to calculate this Jacobin. P(i,j)=B(i,j)*(theta(i)-theta(j)) PP=[P(1,2);P(1,4);P(1,5);P(2,3);P(2,4);P(4,5)] (measurement vector) J=double(jacobian(PP,theta)) getting derivatives of power flow between line i and j in respect to all other theta J = J(:,2:end); omit first column as a slack bus.
Kind Regards, Naime