global dyndata

indmat


[Nz,n] = size([dyndata.x_nom;dyndata.u_nom]);


Ns=6;
As=zeros(Ns,Nz);bs=zeros(Ns,1);

%Constraints 1-5 to fix feed conditions
As(1,7)=1;
As(2,8)=1;
As(3,9)=1;
As(4,10)=1;
As(5,11)=1;
% As(6,11)=1;
% As(7,236+6)=1;

bs(1:5)=dyndata.u_nom([7:11]);

%Constant boilup
As(6,2)=1;

bs(6)=dyndata.u_nom(2);

%Non-negativity of states, helps the optimizer
Ain=-eye(Nz,Nz);bin=zeros(Nz,1);