Diplomarbeit Quantitative Analyse des Ausscheidungs- verhaltens ...

;************** 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
;window,2
;plot,q,av1+int2+int3,/xlog,/ylog
;oplot,q,av1,linestyle=1
;oplot,q,int2,linestyle=2
;oplot,q,int3,linestyle=3
intc=total(av1*q^2*dq)
intg=total(int2*q^2*dq)
intq4=total(int3*q^2*dq)
print,'max q=',q(num-1)
pc=av1(num-1)
pg=int2(num-1)
pq4=int3(num-1)
rpc=intc/pc
rpg=intg/pg
rpq4=intq4/pq4
intpar=[intc,pc,rpc,intg,pg,rpg,intq4,pq4,rpq4]
return,intpar
end
pro nimonifit,q,vector,av
;OP Nov. 2001
;Procedure called by SANSFIT.pro
;
;spherical average of Intensity of
;a rotation-ellipsoidal particle including
;Log-Normal size distribution with fixed width
;
;& model function for g'
;Factor for Carbides
;factor*i0=(N*deltarho)^2
factor=1e-12
;Fittig parameters
i0=vector(0) ;Amplitude
a0=vector(1) ;avg rotational half axis = fitting Parameter
n=vector(2) ;aspect ratio (a>1 Needles, a