figure(1); hold on;
figure(2); hold on;
%UU=[0.5:0.05:1.05];
UU=1;
h=1;

for l=1:length(UU)
U=UU(l);
om=[0.01:0.005:2];
%om=0.25;
a=om;
b=om;
b(49)=0.25-0.12*i
a(49)=0.25-0.12*i;
for j=50:length(om) a(j)=fsolve(@(k) jetvari(k,om(j),U,h), a(j-1)); end;
for j=49:-1:1  a(j)=fsolve(@(k) jetvari(k,om(j),U,h), a(j+1)); end;
for j=50:length(om) b(j)=fsolve(@(k) jetsinu(k,om(j),U,h), b(j-1)); end;
for j=49:-1:1  b(j)=fsolve(@(k) jetsinu(k,om(j),U,h), b(j+1)); end;

figure(1);plot(om,-imag(a));
hold all;
plot(om,-imag(b));

xlabel('\omega');
ylabel('-k_i');
legend('varicose','sinuous');

figure(2);plot(om,real(a));
hold all;
plot(om,real(b));
xlabel('\omega');
ylabel('k_r');
legend('varicose','sinuous');


end;
