makeJac

makeJac(baseMVA, bus, branch, gen, fullJac)

makeJac() - Forms the power flow Jacobian.

J = MAKEJAC(MPC)
J = MAKEJAC(MPC, FULLJAC)
J = MAKEJAC(BASEMVA, BUS, BRANCH, GEN)
J = MAKEJAC(BASEMVA, BUS, BRANCH, GEN, FULLJAC)
[J, YBUS, YF, YT] = MAKEJAC(MPC)

Returns the power flow Jacobian and, optionally, the system admittance
matrices. Inputs can be a MATPOWER case struct or individual BASEMVA,
BUS, BRANCH and GEN values. Bus numbers must be consecutive beginning
at 1 (i.e. internal ordering). If the FULLJAC argument is present and
true, it returns the full Jacobian (sensitivities of all bus injections
w.r.t all voltage angles/magnitudes) as opposed to the reduced version
used in the Newton power flow updates. The units for all quantities are
in per unit with radians for voltage angles.

Note: This function builds the Jacobian from scratch, rebuilding the
      YBUS matrix in the process. You probably don't want to use this
      in performance critical code.