Diplomarbeit Quantitative Analyse des Ausscheidungs- verhaltens ...

f=[0,0,0,0,0,0,0,0]
fp=nimonival(f)
s=fp(0) ; width of logNormal SD
N0=fp(1) ; Number density
num2=fp(2) ;Number of datapoints for sizedist
pi=fp(3)
a=findgen(num2)
a=a*a0/(num2/2)+a0*s
da=a(1)-a(0)
fa=N0/((sqrt(2*pi)*s*a0))*exp(-0.5*(alog(a/a0)/s)^2)
teta=findgen(num)
teta=teta/(num-1)*pi/2
deltateta=teta(1)-teta(0)
int=fltarr(num,num)
intges=fltarr(num,num)
for l=0,num2-1 do begin
for i=0,num-1 do begin
for j=0,num-1 do begin
qr=q(j)*sin(teta(i))
qz=q(j)*cos(teta(i))
u=sqrt((qr*a(l))^2+(n*a(l)*qz)^2)
int(i,j)=factor*fa(l)*((4*pi/3)*(a(l)^3*n))^2*da*((sin(u)u*cos(u))/u^3/3)^2
endfor
endfor
intges=intges+int
endfor
intges=intges/total(fa)
av1=fltarr(num)
for i=0,num-1 do begin
av1(i)=i0*total(intges(*,i)*sin(teta))*deltateta
endfor
;************** Calculate Modelfunction for g' *********
aa=fp(4)
bb=fp(5)
cc=fp(6)
dd=fp(7)
int2=fltarr(num)
for i=0,num-1 do begin
int2(i)=amp2*(aa*(q(i)/qmax)^4/((q(i)/qmax)^4+cc))*bb/(bb+((q(i)/qmax)^2-1+dd)^2)
endfor
;************** Calculate Total intensity *********
;q-4 Contribution
int3=amp3/q^4
av=fltarr(num)
for i=0,num-1 do begin
av(i)=av1(i)+int2(i)+int3(i)
endfor
end
function nimoniparc,q,i,vector
;OP Nov. 2001
;Procedure called by SANSFIT.pro
;
98