%%Compositions  A, B, C, in product streams D, S, R, D1, R1 after optimization
%display('D, S, B')
%display('RL, RV')
%[x_nom(im_xt(1)),x_nom(im_xt(2)),1-x_nom(im_xt(1))-x_nom(im_xt(2));...
    %x_nom(im_xS1(1)),x_nom(im_xS1(2)),1-x_nom(im_xS1(1))-x_nom(im_xS1(2));...
    %x_nom(im_xr(1)),x_nom(im_xr(2)),1-x_nom(im_xr(1))-x_nom(im_xr(2));...
        %x_nom(im_x1(1)),x_nom(im_x1(coldata.Nt(1)+1)),1-x_nom(im_x1(1))-x_nom(im_x1(1+coldata.Nt(2)));...
    %x_nom(im_x2(coldata.Nt(2))),x_nom(im_x2(2*coldata.Nt(2))),1-x_nom(im_x2(coldata.Nt(2)))-x_nom(im_x2(2*coldata.Nt(2)))]
%display('L, VB, RL, RV, S1')
%UU([226 227 228 229 230])
%%Compositions  A, B, C, in product streams D, S, R, D1, R1 after
%%simulation with ode solver
[x(end,im_xt(1)),x(end,im_xt(2)),1-x(end,im_xt(1))-x(end,im_xt(2));...
    x(end,im_xS1(1)),x(end,im_xS1(2)),1-x(end,im_xS1(1))-x(end,im_xS1(2));...
    x(end,im_xr(1)),x(end,im_xr(2)),1-x(end,im_xr(1))-x(end,im_xr(2));...
        x(end,im_x1(1)),x(end,im_x1(coldata.Nt(1)+1)),1-x(end,im_x1(1))-x(end,im_x1(1+coldata.Nt(2)));...
   x(end,im_x2(coldata.Nt(2))),x(end,im_x2(2*coldata.Nt(2))),1-x(end,im_x2(coldata.Nt(2)))-x(end,im_x2(2*coldata.Nt(2)))]
