;$Id: pfit.pro,v 1.1 2024/05/24 15:33:33 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=12,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
!p.charsize=1.7
!x.margin=[7.8,.5]
!y.margin=[3.2,.5]
!x.title='!8t!6'
!y.title='!7n!6!dM!n'
;
plot_oo,dt,1./kI,yr=[8e-3,1]
fit1=.13*2.80E-04^.33333*dt^.66667
fit2=2.3e-2*dt^(4./9.)
;oplot,dt,.13*2.80E-04^.33333*dt^.66667,col=122
;oplot,dt,2.3e-2*dt^(4./9.),col=155            
;oplot,dt,2.3e-2*dt^(4./9.),col=155            
;oplot,dt,.13*2.80E-04^.33333*dt^.66667+.12*1e-10^(1./9.)*dt^(4./9.),col=55,th=3
nexp=4.
;oplot,dt,((.13*2.80E-04^.33333*dt^.66667)^nexp+(.12*1e-10^(1./9.)*dt^(4./9.))^nexp)^(1./nexp),col=55,th=3
oplot,dt,fit1,col=122
oplot,dt,fit2,col=155
oplot,dt,fit1+fit2,col=55
oplot,dt,(fit1^nexp+fit2^nexp)^(1./nexp),col=55,th=3,li=2
;
print,"$mv idl.ps ~/Overleaf/Axel/Relax-Hosking/fig/pfit.eps"
END