load Isaac_bright_oh.dat
x=Isaac_bright_oh(:,1)
y=Isaac_bright_oh(:,2)
load atmosphere.dat
a=atmosphere(:,1)
b=atmosphere(:,2)
[AX,H1,H2] = plotyy(a,b,x,y,'plot','stem')
axes(AX(1))
axis('autoy')
axes(AX(2))
axis('autoy')
set(H2,'Marker','none')
axis([1,2,0,1])
