function ss = find_ss()

close all;

% parameters
a = 1.5;
b = 0.01;
c = 0.5;
d = 1.0;
r = 1.0;
eps = 0.009;

k = [100:10:1200];

for j = 1:length(k)
    M = [r/k(j) r/k(j) b; r/k(j)*a r/k(j) (b-eps); c*b c*(b-eps) 0];
    ss = inv(M)*[r;r;d];
    plot3(ss(1),ss(2),ss(3),'.k');
    %plot(k(j),ss(2));
    hold on;
end

grid on;
xlabel('N_1^*'); ylabel('N_2^*'); zlabel('P^*');
%xlabel('k');ylabel('N_2');