home *** CD-ROM | disk | FTP | other *** search
/ c't freeware shareware 1999 February / CT_SW9902.ISO / pc / software / wissen / daten / gnuplot.hqx / gnuplot.2.0b4 / Scripts / demo / reflect.fnc < prev    next >
Text File  |  1997-05-03  |  897b  |  27 lines
  1. #
  2. #   Model function for Reflectivity evaluation
  3. #
  4.  
  5. mu = 1.130469005513490E-001                     # (cm-1) @ 17.479 keV
  6. t0 = 0.18                                       # cm
  7. tb = 11.417823202820120 * 0.01745329251994      # thetaB (radians)
  8. A = mu * t0 / cos(tb)
  9. P = (1 + (cos(2.*tb))**2) / 2
  10. Fhkl = sqrt(3.536346308456155**2 + (4.58815426260982e-4)**2) * 0.968
  11. r0 = 2.81794092e-13                             # classical electron radius
  12. lambda = 7.09338062818239e-9                    # Mo K in cm
  13. V = 1.62253546981499e-23
  14. P = (1. + (cos(2.*tb))**2) / 2.
  15.  
  16. #
  17.  
  18. W(x) = 1./(sqrt(2.*pi)*eta) * exp( -1. * x**2 / (2.*eta**2) )
  19. Y(tc) = tc/sin(tb) * (Fhkl*r0*lambda) / V
  20. f(tc)= (tanh(Y(tc)) + abs(cos(2.*tb)) * tanh(abs(Y(tc)*cos(2.*tb)))) / (Y(tc)*(1.+(cos(2.*tb))**2))
  21. Q(tc) = (r0*Fhkl/V)**2 * (lambda**3/sin(2.*tb)) * P * f(tc)
  22. a(x) = W(x) * Q(tc) / mu
  23.  
  24. #
  25.  
  26. R(x) = sinh(A*a(x)) * exp(-1.*A*(1.+a(x)))
  27.