function Ws = shaft_work(x, aux)
sc = aux.sc;
i_com = length(sc.com.x);
i_gco = length(sc.gco.x);
i_val = length(sc.val.x);
i_vap = length(sc.vap.x);
i_ihx = length(sc.ihx.x);
com = x(1:i_com);
gco = x(i_com+1:i_com+i_gco);
val = x(i_com+i_gco+1:i_com+i_gco+i_val);
vap = x(i_com+i_gco+i_val+1:i_com+i_gco+i_val+i_vap);
ihx = x(i_com+i_gco+i_val+i_vap+1:end);

Ws = com(1)*val(2);
