To show a highly output multiplicity I want to plot a graph between output temperature and coolant flow rate (input) in open loop. CSTR have multiple steady state point but I don't know how to plot it in MATLAB. I have tried it like this. pl correct me.

Da = 0.072; gamma = 20; B = 8; eta = 0.3; Ts=0.1;

f=@(u,x)[-x(1)+ Da*(1-x(1))*exp(gamma*x(2)/(gamma+ x(2))); -x(2)+B*Da*(1-x(1))*exp(gamma*x(2)/(gamma+ x(2)))+ eta*(u-x(2))];

options = odeset('RelTol',1e-7,'AbsTol',[1e-7 1e-7]);

[u,xa] = ode45(f,[-3 3],[0.044975553686829; -0.41559367683716],options); plot(u,xa(:,2),'-.')

but plot shows only one steady state point.