RosettaCodeData/Task/Resistor-mesh/Mathematica/resistor-mesh-2.math

8 lines
229 B
Text
Raw Permalink Normal View History

2023-07-01 11:58:00 -04:00
ResistanceMatrix[g_, prec_:$MachinePrecision]:= With[{m = PseudoInverse[N[KirchhoffMatrix[g], prec]]},
Outer[Plus, Diagonal[m], Diagonal[m]] - m - Transpose[m]
]
rm = ResistanceMatrix[GridGraph[{10, 10}], 40];
rm[[12, 68]]