share/contrib/fresnel/fg.mac

   1 matchdeclare(mv,true)$
   2 defrule(fg1,cfresnel(mv),(1/2+f(mv)*sin((%pi/2)*mv^2)-g(mv)*cos((%pi/2)*mv^2)))$
   3 defrule(fg2,sfresnel(mv),(1/2-f(mv)*cos((%pi/2)*mv^2)-g(mv)*sin((%pi/2)*mv^2)))$
   4 defrule(ga,g(mv),1/((%pi*mv)*(%pi*(mv)^2)))$
   5 defrule(fa,f(mv),1/(%pi*mv)+3/(%pi*(mv^2))^2)$
   6 batch("randi.mac")$
   7 randi:apply1(randi,fg1)$
   8 randi:apply1(randi,fg2)$
   9 randi:expand(randi)$
  10 randi:trigreduce(randi)$
  11 randi:rncombine(randi)$
  12 randi:map(factor,randi)$
  13 randi:rncombine(randi);
  14 randi:apply1(randi,fa,ga);
  15 limit(randi,w,16000000 * %pi);
  16 plot2d (d15, [w,float(%pi)*15700000,float(%pi)*16100000], [logy, true]);
  17 plot2d (d15, [w,float(%pi)*1000000,float(%pi)*15700000], [logy, true]);
  18 plot2d (d15, [w,float(%pi)*100000,float(%pi)*1000000], [logy, true]);
  19 plot2d (d15, [w,float(%pi)*10000,float(%pi)*100000], [logy, true]);
  20 plot2d (d15, [w,float(%pi)*1000,float(%pi)*10000], [logy, true]);
  21 plot2d (d15, [w,float(%pi)*100,float(%pi)*1000], [logy, true]);