function xdot=stdpte(t,x)
% EDO de la planta termoelectrica, Standard Model, Usoro, MIT, 1977
% "A drum boiler-turbine power plant under emergency state control"
% eqm 8feb'19
%------------------------------------------------------------------
% equivalence
% x0=[heco, rdrs, twwm, rsso, hsso, rrho, hrho,...
%     ntr, hlho, hdew, hhho, c3md, c5ar, c5fl,...
%     c2gr, c5rh, c5sy, card, cfld, cgrd, cwryd,...
%     cxggd, csyd, c6tr, c9tr, c12tr, delta];

% dynamic equations
global dxdt
xdot=dxdt;
% xdot(1)=(wfw*(hhho-heco)+qec)/kmece;
% xdot(2)=frdrs;
% xdot(3)=(qwwgm-qwwmw)/(kmwwme*kswwm);
% xdot(4)=(wdrs+wsy-wsso)/kvsh;
% xdot(5)=(wdrs*hdrs+wsy*hfpo-wsso*hsso+qsh)/kmshe;
% xdot(6)=(whpo+wry-wrho)/kvrh;
% xdot(7)=(whpo*hhpo+wry*hfpc-wrho*hrho+qrh)/kmrhe;
% xdot(8)=(mwtro-mwgn)/(ntr*kjtre);
% xdot(9)=(hcpo+hlhd-hlho)/ktclh;
% xdot(10)=(hlho+hded-hdew)/ktcde;
% xdot(11)=(hfpo+hhhd-hhho)/ktchh;
% xdot(12)=c2md/ktc1md;
% xdot(13)=c4ar/ktc1ar;
% xdot(14)=c4fl/ktc1fl;
% xdot(15)=c1gr*kc1gr*kcgr;
% xdot(16)=c4rh/ktc1rh;
% xdot(17)=c4sy/ktc1sy;
% xdot(18)=(c6ar-card)/ktc2ar;
% xdot(19)=(c6fl-cfld)/ktc2fl;
% xdot(20)=(c2gr-cgrd)/ktc1gr;
% xdot(21)=(cwry-cwryd)/ktc2rh;
% xdot(22)=(cxgg-cxggd)/ktc3rh;
% xdot(23)=(c6sy-csyd)/ktc2sy;
% xdot(24)=c2tr/ktc1tr;
% xdot(25)=(c8tr-c9tr)/ktc2tr;
% xdot(26)=(c10tr-c12tr)/ktc3tr;
% xdot(27)=kc1gn*(ntr-kntr);
end

