function h=hfn(x)
% EE362G Smart Grids
% Nonlinear Measurement Function for State Estimation on 2-bus network 

h=[x(1); x(2); 15*x(1)*x(2)*sin(x(3)); ... 
    15*x(2)^2-15*x(1)*x(2)*cos(x(3))];
