% Emacs mode header: -*-Octave-*- source. %

tw_nonlin_4
yy=abs(rh);
F_int_convex
F_lambda
recc=1000./(3.5*omega.*lambda(1:8));
latcum=cumtrapz(-x(ilist(1:8)),recc);
hold on
grid on
plot(-x,sqrt(rh.*rh+ih.*ih)/40,'g')
plot(-x,rh/40,'b')
plot(-x(ilist(1:8)),3.5*omega.*lambda(1:8)/50,'r')
plot(-x(ilist(1:8)),recc,'m')
plot(-x(ilist(1:8)),recc,'mo')
plot(-x(ilist),list,'ko')
plot(-x(ilist(1:8)),latcum,'k')
title(' Numerical calculation of latency and velocity')
