fx=sz_fx(nu) Spectral distortion function f(x) due to SZ. x is dimensionless frequency x=(h*nu)/(k*T) e.g. nu=[1:1:300]*1e9; plot(nu,sz_fx(nu)); see also sz_relcor.m
0001 function fx=sz_fx(nu) 0002 % fx=sz_fx(nu) 0003 % 0004 % Spectral distortion function f(x) due to SZ. 0005 % 0006 % x is dimensionless frequency x=(h*nu)/(k*T) 0007 % 0008 % e.g. 0009 % 0010 % nu=[1:1:300]*1e9; 0011 % plot(nu,sz_fx(nu)); 0012 % 0013 % see also sz_relcor.m 0014 0015 h=6.626e-34; 0016 k=1.38e-23; 0017 T=2.726; 0018 0019 x=(h*nu)/(k*T); 0020 0021 fx=x.*(exp(x)+1)./(exp(x)-1)-4; 0022 0023 return