A method for determining sensitivity coefficients of an electric power network comprising a set of nodes and a set of branches, the electric power network further being provided with a monitoring infrastructure comprising metering units at each one of a plurality of nodes (called measuring nodes (N1,..., N4)) in the set of nodes, each of one of said metering units being arranged to measure a nodal voltage of one of said measuring nodes, branch currents flowing into or out of said one of said measuring nodes, and respective phase differences between the branch currents and the nodal voltage, said branch currents either flowing through branches of the network that are incident on said one of said measuring nodes or being associated with power injections at said one of said measuring nodes, the monitoring infrastructure further comprising at least one processing unit and a communication infrastructure arranged for allowing communication between the metering units and said at least one processing unit, the method determines current sensitivity coefficients of a plurality of selected branches (A,..., D), with respect to the plurality of measuring nodes (N1,..., N4) of the electric power network, by carrying out the following steps: I. Have the metering units measure at the same time, at each one of said measuring nodes (N1,..., N4), repeatedly over a time window (τ), sets of data comprising the nodal voltage V n ( t) and values of branch currents I b ( t) flowing into or out of said one of said measuring nodes, timestamp t ∈ { t 1,..., t m} the measured sets of data, and compute a timestamped nodal active power P n ( t) and a timestamped nodal reactive power Q n ( t) from each set of measured data; II. for each one of a selection of branches (called the selected branches), compute variations (Δ Ĩ b ( t)) of the branch current measured in step I (box 01) flowing through said selected branch, by subtracting from the measured value of the branch current flowing through said selected branch, the precedent measured value of the current through the same selected branch, and for each one of said measuring nodes, compute concomitant variations (Δ Ĩ b ( t),Δ P̃ n ( t);Δ Q̃ n ( t)) of the current measured in step I (box "01") and of the nodal active and reactive powers computed in step I (box "01"), by subtracting from each set of concomitant values of the current, the active power and the reactive power respectively, the precedent value of the same variables; III. compile chronologically ordered tables of the variations of the current (Δ Ĩ b ( t)) through each selected branch (A,..., D) in relation to concomitant variations of the active (Δ P̃ 1 ( t),..., Δ P̌ N ( t)) and reactive powers (Δ Q̃ 1 ( t),...,Δ Q̃ N ( t)) at all measuring nodes (N1,..., N4); IV. perform a Maximum Likelihood Estimation of current sensitivity coefficients ( KIP bn, KIQ bn) linking variations of the current through the selected branches (Δ Ĩ b ( t)), as compiled during step III (box "02"), to the nodal active and reactive power variations (Δ P n ( t),Δ Q n ( t)), while taking into account serial correlation between error terms corresponding to discrepancies between the actual variations and the variations predicted by the Maximum Likelihood Estimation, and obtain from the determined current sensitivity coefficients ( KlP bn, KIQ bn) current sensitivity coefficient matrices. I. Have the metering units measure at the same time, at each one of said measuring nodes (N1,..., N4), repeatedly over a time window (τ), sets of data comprising the nodal voltage V n ( t) and values of branch currents I b ( t) flowing into or out of said one of said measuring nodes, timestamp t ∈ { t 1,..., t m} the measured sets of data, and compute a timestamped nodal active power P n ( t) and a timestamped nodal reactive power Q n ( t) from each set of measured data; II. for each one of a selection of branches (called the selected branches), compute variations (Δ Ĩ b ( t)) of the branch current measured in step I (box 01) flowing through said selected branch, by subtracting from the measured value of the branch current flowing through said selected branch, the precedent measured value of the current through the same selected branch, and for each one of said measuring nodes, compute concomitant variations (Δ Ĩ b ( t),Δ P̃ n ( t);Δ Q̃ n ( t)) of the current measured in step I (box "01") and of the nodal active and reactive powers computed in step I (box "01"), by subtracting from each set of concomitant values of the current, the active power and the reactive power respectively, the precedent value of the same variables; III. compile chronologically ordered tables of the variations of the current (Δ Ĩ b ( t)) through each selected branch (A,..., D) in relation to concomitant variations of the active (Δ P̃ 1 ( t),..., Δ P̌ N ( t)) and reactive powers (Δ Q̃ 1 ( t),...,Δ Q̃ N ( t)) at all measuring nodes (N1,..., N4); IV. perform a Maximum Likelihood Estimation of current sensitivity coefficients ( KIP bn, KIQ bn) linking variations of the current through the selected branches (Δ Ĩ b ( t)), as compiled during step III (box "02"), to the nodal active and reactive power variations (Δ P n ( t),Δ Q n ( t)), while taking into account serial correlation between error terms corresponding to discrepancies between the actual variations and the variations predicted by the Maximum Likelihood Estimation, and obtain from the determined current sensitivity coefficients ( KlP bn, KIQ bn) current sensitivity coefficient matrices.