%
% HPastar.m
%
p=pwd

[ap,Fname]=adhtproc;
cd(p)

adhbproc
cd(p)

ast=ap-ad;

%display
figure
hold on
for i=1:n
  plot(350:2:800,ap(:,i),'r')
end
hold off  
title('ab total')

%display
figure
hold on
for i=1:n
  plot(350:2:800,ad(:,i),'b')
end
hold off  
title('ad, Bleached')

%display

figure
hold on
for i=1:n
 plot(350:2:800,ast(:,i),'g')
end
hold off
title('astar')