Diplomarbeit Quantitative Analyse des Ausscheidungs- verhaltens ...

print,vec
print,sigma
plot,q,i,/ylog,/xlog,psym=6,title=ti
oplot,q,fit
;************** Calculate functions with fitted Parameters *********
;definitions
amp2=vec(3)
qmax=vec(4)
;Call fixed values from Nimonival
f=[0,0,0,0,0,0,0,0]
fp=nimonival(f)
pi=fp(3)
aa=fp(4)
bb=fp(5)
cc=fp(6)
dd=fp(7)
;-------------------------------------------------------
;Model function for g'
x=q/qmax
int2=amp2*(aa*x^4/(x^4+cc))*bb/(bb+(x^2-1+dd)^2)
oplot,q,int2,linestyle=3
;q-4 contribution
int4=vec(5)/q^4
oplot,q,int4,linestyle=3
;Scattering from carbides
int3=fit-int2-int4
oplot,q,int3,linestyle=3
;calculate integral parameters
intpar=nimonicalc(vec)
print,intpar
;write results to file
OPENW,unt,output2,/GET_LUN
PRINTF,unt,'Ic(1e36 cm-1),Rc,nc,Ig,qmax,Aq4'
PRINTF,unt,vec(0)
PRINTF,unt,vec(1)
PRINTF,unt,vec(2)
PRINTF,unt,vec(3)
PRINTF,unt,vec(4)
PRINTF,unt,vec(5)
PRINTF,unt,'Qc,Pc,Rpc,Qg,Pg,Rpg,Qq4,Pq4,Rpq4'
PRINTF,unt,intpar(0)
PRINTF,unt,intpar(1)
PRINTF,unt,intpar(2)
PRINTF,unt,intpar(3)
PRINTF,unt,intpar(4)
PRINTF,unt,intpar(5)
PRINTF,unt,intpar(6)
PRINTF,unt,intpar(7)
PRINTF,unt,intpar(8)
FREE_LUN,unt
OPENW,unt,output1,/GET_LUN
FOR iii=0,N_elements(q)-1 DO BEGIN
PRINTF,unt,FORMAT='(f12.6,1x,f12.6,1x,f12.6,1x,f12.6,1x,f12.6,1x,f12.6)',$
q(iii),i(iii),fit(iii),int2(iii),int3(iii),int4(iii)
endfor
FREE_LUN,unt
end
95