dens_ini
amp = 0.1*us;
bias = us;
T = [0:0.5:200]';
n = size(T,1);
u1 = amp(1)*prbs(n,10,100,0.8,1)+bias(1);
u2 = amp(2)*prbs(n,10,100,0.8,1)+bias(2);
u3 = amp(3)*prbs(n,10,100,0.8,1)+bias(3);
u4 = amp(4)*prbs(n,10,100,0.8,1)+bias(4);

u = [u1 u2 u3 u4];

op=simset('solver','ode15s','InitialState',x0);

[t,x,y]=sim('dens_io', T, op, [T u]);

plot(t,y);
