restore, './dif20000/disp_dif.sav'
restore, './kdv20000/disp_kdv.sav'
n=20
fm1=dblarr(n+1)
fm2=dblarr(n+1)

for h=0,n do begin
fm1(h)=max(fph(*,h))
fm2(h)=max(fan(*,h))
endfor

loadct,39

set_plot,'ps'
device,bits_per_pixel=8,/encapsulated,/color,filename='fig7.eps'

!P.MULTI = [0, 2, 2]
 plot,xph,fphd(*,0),xrange=[min(xph)/10.0,max(xph)/10.0],xstyle=1,yrange=[0,1],ystyle=1,charsize=1.2,xtitle='!4n!X',ytitle='A!dz!n'
oplot,xph,fphd(*,n/20),color=60;blue
oplot,xph,fphd(*,n),color=240 ;red

plot,t,fm,/xlog,/ylog,xrange=[1,12.5],xstyle=1,yrange=[0.1,1],ystyle=1,charsize=1.2,xtitle='s',ytitle='A!dz!n!U*!N'
oplot,t,exp(r[0])*t^r[1], psym=5

 plot,xph,fph(*,0),xrange=[min(xph)/10.0,max(xph)/10.0],xstyle=1,yrange=[-0.4,1],ystyle=1,charsize=1.2,xtitle='!4n!X',ytitle='A!dz!n'
oplot,xph,fph(*,n/20),line=0,color=60 ;blue
oplot,xph,fan(*,n/20),nsum=20,color=60,psym=4 ;blue
oplot,xph,fph(*,n),line=0,color=240 ;red
oplot,xph,fan(*,n),nsum=50,color=240,psym=4 ;red

plot,t,fm1,/xlog,/ylog,xrange=[1,12.5],xstyle=1,yrange=[0.25,1],ystyle=1,charsize=1.2,xtitle='s',ytitle='A!dz!n!U*!N'
oplot,t,fm2,line=2
oplot,t,exp(r1[0])*t^r1[1], psym=5
oplot,t,exp(r2[0])*t^r2[1], psym=1
oplot,t,t^(-1.0/3.0)*sqrt(2.0*!pi)/4.06235,line=0,color=240 ;red
oplot,t,t^(-1.0/3.0)*sqrt(!pi)/4.06235,line=0,color=60 ;blue

!P.MULTI = 0

device,/close

print,r
print,r1
print,r2

end