#!/usr/bin/octave
% Octave file: a figure for the file StabMiso.tex

gset output "kdv-spectrum.tex"
% Clear labels and arrows:
gset nolabel
%gset noarrow

gset term pstex norotate
gset nokey
gset zeroaxis
gset nogrid
gset noborder
gset ticscale 1
gset noxtics
gset noytics
gset size 0.66,0.66
#gset linestyle 1 linewidth 3

k=(0.:0.1:0.7)';
#gset label '$\\sigma\\sb{ess}(e^{-\\mu x}J\\mathcal{H}\\sb{c}e^{\\mu x})$' at 0.4,0.6;
#gset label '$\\mu=\\frac 1{10}$' at -0.25,-0.55
#gset label '$\\mu=\\frac 1 4$' at -0.45,-0.55

gset label '${\\rm Re}\\,\\lambda$' at 0.85,0.1
gset label '${\\rm Im}\\,\\lambda$' at 0.05,0.95

k=(-2:0.02:2)';
xaxis=[k,0.*k];
yaxis=[0.*k,k];

mu=0.1;
#FIG1=[-mu.^3-3.*mu.*k.^2,k.^3+3.*mu.^2.*k];
FIG1=[mu.^3-3.*mu.*k.^2-mu,k.^3-3.*mu.^2.*k+k];
mu=0.2;
FIG2=[mu.^3-3.*mu.*k.^2-mu,k.^3-3.*mu.^2.*k+k];
mu=0.3;
FIG3=[mu.^3-3.*mu.*k.^2-mu,k.^3-3.*mu.^2.*k+k];
mu=0.4;
FIG4=[mu.^3-3.*mu.*k.^2-mu,k.^3-3.*mu.^2.*k+k];
mu=0.4;
FIG4=[mu.^3-3.*mu.*k.^2-mu,k.^3-3.*mu.^2.*k+k];
mu=0.5;
#FIG3=[-mu.^3-3.*mu.*k.^2,k.^3+3.*mu.^2.*k];
FIG5=[mu.^3-3.*mu.*k.^2-mu,k.^3-3.*mu.^2.*k+k];
mu=0.6;
#FIG3=[-mu.^3-3.*mu.*k.^2,k.^3+3.*mu.^2.*k];
FIG6=[mu.^3-3.*mu.*k.^2-mu,k.^3-3.*mu.^2.*k+k];
mu=0.65;
#FIG3=[-mu.^3-3.*mu.*k.^2,k.^3+3.*mu.^2.*k];
FIG666=[mu.^3-3.*mu.*k.^2-mu,k.^3-3.*mu.^2.*k+k];
mu=0.7;
FIG7=[mu.^3-3.*mu.*k.^2-mu,k.^3-3.*mu.^2.*k+k];
mu=0.8;
FIG8=[mu.^3-3.*mu.*k.^2-mu,k.^3-3.*mu.^2.*k+k];

PTS=[0,0];

gplot [-1.2:1.2] [-1.2:1.2] FIG1 with lines 1, \
 FIG666 with lines 3, \
    xaxis with lines 1, \
    yaxis with lines 1, \
    PTS using 1:2 with points 2;