[x11,x22]
;  Input, Output
; Input is spatial frequency in Angstroms,
; The filter is scaled to have 1 at 100A
;  X-ray scattering power spectrum is in
; /net/saba/usr1/pawel/nenad/1D-PW/scattering8.tst
X51=1
X52=750
X53=1500
DO LB1 X51=1,1500
UD IC,X52,X31,X22
scattering8
X52=X52+1
IF(X52.gt.1500) GOTO LB2
UD IC,X52,X32,X23
scattering8

IF(X11.LE.X31) THEN
IF(X11.GE.X32)  THEN
X77=(X31-X11)**2
X78=(X32-X11)**2
IF(X78.LT.X77) X22=X23
GOTO LB2
ENDIF
ENDIF
IF(X11.LT.X32) THEN
X51=X52
X52=INT((X53-X52)/2+X52)
ELSE
X53=X52
X52=INT((X52-X51)/2+X51)
ENDIF
LB1
LB2
X22=X22/37891.
;E4
RE