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

%load act_data.dat %x_act rh_act
load pas_data.dat %x_pas ah_pas

subplot(2,1,1), plot(-x_pas,ah_pas), grid on
title('Frank-s sources, passive and active wave')

% x_act=x;
% rh_act=rh;

subplot(2,1,2), plot(-x_act,rh_act), grid on

% envelope % plot(-x,sqrt(rh.*rh+ih.*ih),'k.')