>I need a good function or routine (or algoritm) to make a deconvolution of
>Kalpha2, from x,y points or from function expression.
Alpha 1 and alpha 2 components are added, so that this will not be
a deconvolution, but a desummation. You will find the basic algorithm
in the Fortran source code of the AXRLEC program, subroutine STRIP2 :
http://sdpd.univ-lemans.fr/also.html
SUBROUTINE STRIP2(BUF,RA,IL,NBMES,TH0,STP)
DIMENSION DLAM(4),BUF(40000)
DATA PIF,DLAM/8.726646E-3,2.48484E-3,2.1716E-3,2.034E-3,1.71E-3/
DTH=STP*PIF
DLAM2=DLAM(IL)*.5
IF(DTH.LE.0.)RETURN
DO 10 J=1,NBMES
TH=(TH0+(J-1)*STP)*PIF
TG=SIN(TH)/COS(TH)
DX=TG*DLAM(IL)*(1.+DLAM2*TG**2)/DTH
INC=DX
I=J+INC
IF(I.GT.NBMES)RETURN
REST=INC-DX+1.
DY=BUF(J)*RA
BUF(I)=AMAX1(0.,BUF(I)-DY*REST)
IF(I.GE.NBMES)RETURN
I=I+1
10 BUF(I)=AMAX1(0.,BUF(I)-DY*(1.-REST))
END
Best,
Armel Le Bail
http://sdpd.univ-lemans.fr/course/
