fcall.m

这是一个很适合研究和学习用的电力系统仿真软件

function fcall(a)

global DAE

if ~a.n, return, end

x = DAE.x(a.x);
V1 = DAE.V(a.bus);
iTf = 1./a.con(:,5);

DAE.f(a.x) = -(V1.*iTf+x).*iTf;