function [y] = out_eq3(xc,xd,t)
%
%  Output model of [H9.7]
%
RLoad = 1;
RSh = 10;
LLoad = 0.01;
V0 = 750;
p = 16 + 2/3;
tp = 1/p;
%
iL = xc(1);
m0 = xd(1);
m1 = xd(2);
ms = xd(3);
%
vLine = V0*sin(2*pi*t/tp);
inpt = [iL;vLine];
a11 = RSh*(1 - m0) + m0;
a12 = 1;
a21 = m0*RLoad;
a22 = RLoad + RSh;
A = [a11,a12;a21,a22];
b11 = RSh;
b12 = 1;
b21 = RLoad*RSh;
b22 = RLoad;
B = [b11,b12;b21,b22];
tear = A\B*inpt;
s = tear(1);
vRLoad = tear(2);
vTh = m0*s;
vLoad = vLine - vTh;
vRSh = vLoad - vRLoad;
iRSh = vRSh/RSh;
iLoad = iL + iRSh;
diL = vRSh/LLoad;
%
y = zeros(2,1);
y(1) = vLoad;
y(2) = iLoad;
%
return