%  This function solves project [P5.1], part B
%  -------------------------------------------
%
echo on
%
%  Stability domain of Houbolt
%
[x1,y1] = stabdom4(2,+1);
%
plot(x1,y1)
axis equal
grid on;
title('ttt')
xlabel('xxx')
ylabel('yyy')
%
print -deps nsds_pr5_1b.eps
%
echo off
return
