simulate non linear system

how can I simulate (48)? usually I used the function lsim, but in this case I have only matrices A and B.

[darova]

for k = 1:n-1
    dx(k+1) = A*dx(k) + B*du(k) + d(k);
end