$SPAD/src/input fresnel.input

Timothy Daly December 3, 2016

Abstract

Contents

2

— * —

)set break resume
)spool fresnel.output
)set message test on
)set message auto off
)clear all
--S 1 of 4
fC(z:FLOAT):FLOAT ==
  z < 0 => error "fC not defined for negative argument"
  z = 0 => 0
  n := 10
  sqrt((2/%pi)*z)*_
   reduce(+,[(-1)^k*z^(2*k)/(factorial(2*k)*(4*k+1)) for k in 0..n])
--E 1
--S 2 of 4
--[[z/100.0,fC(z/100.0)] for z in 1..5000]
--E 2
--S 3 of 4
fS(z:FLOAT):FLOAT ==
  z < 0 => error "fS not defined for negative argument"
  z = 0 => 0
  n := 10
  sqrt((2/%pi)*z)*_
   reduce(+,[(-1)^k*(z^(2*k+1))/(factorial(2*k+1)*(4*k+3)) for k in 0..n])
--E 3
--S 4 of 4
--[[z/100.0,fS(z/100.0)] for z in 1..5000]
--E 4
--S
pearceyC:=_
[ [0.00, 0.0000000], [0.25, 0.3964561], [0.50, 0.5502472], [0.75, 0.6531193],_
  [1.00, 0.7217059], [1.25, 0.762404],  [1.50, 0.779084],  [1.75, 0.774978],_
  [2.00, 0.753302],  [2.25, 0.717446],  [2.50, 0.670986],  [2.75, 0.617615],_
  [3.00, 0.561020],  [3.25, 0.504745],  [3.50, 0.452047],  [3.75, 0.405762],_
  [4.00, 0.368193],  [4.25, 0.341021],  [4.50, 0.325249],  [4.75, 0.321186],_
  [5.00, 0.328457],  [5.25, 0.346058],  [5.50, 0.372439],  [5.75, 0.405610],_
  [6.00, 0.443274],  [6.25, 0.482966],  [6.50, 0.522202],  [6.75, 0.558620],_
  [7.00, 0.590116],  [7.25, 0.614951],  [7.50, 0.631845],  [7.75, 0.640034],_
  [8.00, 0.639301],  [8.25, 0.629969],  [8.50, 0.612868],  [8.75, 0.589271],_
  [9.00, 0.560804],  [9.25, 0.529344],  [9.50, 0.496895],  [9.75, 0.465469],_
 [10.00, 0.436964], [10.25, 0.413053], [10.50, 0.395087], [10.75, 0.384027],_
 [11.00, 0.380390], [11.25, 0.384231], [11.50, 0.395149], [11.75, 0.412319],_
 [12.00, 0.434555], [12.25, 0.460384], [12.50, 0.488146], [12.75, 0.516096],_

[13.00, 0.542511], [13.25, 0.565798], [13.50, 0.584583], [13.75, 0.597795],_
[14.00, 0.604721], [14.25, 0.605048], [14.50, 0.598871], [14.75, 0.586682],_
[15.00, 0.569335], [15.25, 0.547984], [15.50, 0.524009], [15.75, 0.498930],_
[16.00, 0.474310], [16.25, 0.451659], [16.50, 0.432343], [16.75, 0.417502],_
[17.00, 0.407985], [17.25, 0.404300], [17.50, 0.406589], [17.75, 0.414627],_
[18.00, 0.427837], [18.25, 0.445331], [18.50, 0.465972], [18.75, 0.488443],_
[19.00, 0.511332], [19.25, 0.533222], [19.50, 0.552774], [19.75, 0.568812],_
[20.00, 0.580389], [20.25, 0.586847], [20.50, 0.587849], [20.75, 0.583401],_
[21.00, 0.573842], [21.25, 0.559824], [21.50, 0.542266], [21.75, 0.522293],_
[22.00, 0.501167], [22.25, 0.480207], [22.50, 0.460707], [22.75, 0.443854],_
[23.00, 0.430662], [23.25, 0.421906], [23.50, 0.418080], [23.75, 0.419367],_
[24.00, 0.425635], [24.25, 0.436444], [24.50, 0.451078], [24.75, 0.468594],_
[25.00, 0.487880], [25.25, 0.507725], [25.50, 0.526896], [25.75, 0.544215],_
[26.00, 0.558626], [26.25, 0.569272], [26.50, 0.575524], [26.75, 0.577038],_
[27.00, 0.573766], [27.25, 0.565954], [27.50, 0.554127], [27.75, 0.539054],_
[28.00, 0.521695], [28.25, 0.503146], [28.50, 0.484566], [28.75, 0.467104],_
[29.00, 0.451832], [29.25, 0.439675], [29.50, 0.431359], [29.75, 0.427366],_
[30.00, 0.427908], [30.25, 0.432913], [30.50, 0.442034], [30.75, 0.454673],_
[31.00, 0.470019], [31.25, 0.487100], [31.50, 0.504844], [31.75, 0.522148],_
[32.00, 0.537944], [32.25, 0.551266], [32.50, 0.561307], [32.75, 0.567471],_
[33.00, 0.569407], [33.25, 0.567026], [33.50, 0.560508], [33.75, 0.550288],_
[34.00, 0.537026], [34.25, 0.521566], [34.50, 0.504881], [34.75, 0.488015],_
[35.00, 0.472012], [35.25, 0.457857], [35.50, 0.446415], [35.75, 0.438375],_
[36.00, 0.434212], [36.25, 0.434156], [36.50, 0.438182], [36.75, 0.446014],_
[37.00, 0.457140], [37.25, 0.470848], [37.50, 0.486272], [37.75, 0.502444],_
[38.00, 0.518359], [38.25, 0.533031], [38.50, 0.545560], [38.75, 0.555182],_
[39.00, 0.561321], [39.25, 0.563619], [39.50, 0.561957], [39.75, 0.556463],_
[40.00, 0.547503], [40.25, 0.535653], [40.50, 0.521665], [40.75, 0.506420],_
[41.00, 0.490870], [41.25, 0.475980], [41.50, 0.462670], [41.75, 0.451755],_
[42.00, 0.443897], [42.25, 0.439565], [42.50, 0.439006], [42.75, 0.442234],_
[43.00, 0.449025], [43.25, 0.458938], [43.50, 0.471341], [43.75, 0.485450],_
[44.00, 0.500382], [44.25, 0.515205], [44.50, 0.529002], [44.75, 0.540923],_
[45.00, 0.550239], [45.25, 0.556387], [45.50, 0.559004], [45.75, 0.557947],_
[46.00, 0.553301], [46.25, 0.545374], [46.50, 0.534676], [46.75, 0.521883],_
[47.00, 0.507802], [47.25, 0.493312], [47.50, 0.479313], [47.75, 0.466670],_
[48.00, 0.456160], [48.25, 0.448425], [48.50, 0.443930], [48.75, 0.442936],_
[49.00, 0.445486], [49.25, 0.451406], [49.50, 0.460311], [49.75, 0.471633],_
[50.00, 0.484658]]
--E
--S
[[x.1,x.2,fC(x.1),fC(x.1)-x.2] for x in perceyC]
--E
--S
pearceyS:=_
[ [0.00, 0.0000000], [0.25, 0.0330970], [0.50, 0.0923658], [0.75, 0.1659294],_
 [1.00, 0.2475583], [1.25, 0.332216],  [1.50, 0.415483],  [1.75, 0.493469],_
 [2.00, 0.562849],  [2.25, 0.620944],  [2.50, 0.665787],  [2.75, 0.696174],_
 [3.00, 0.711685],  [3.25, 0.712666],  [3.50, 0.700180],  [3.75, 0.675925],_

 [4.00, 0.642119],  [4.25, 0.601362],  [4.50, 0.556489],  [4.75, 0.510408],_
 [5.00, 0.465942],  [5.25, 0.425677],  [5.50, 0.391834],  [5.75, 0.366161],_
 [6.00, 0.349852],  [6.25, 0.343503],  [6.50, 0.347099],  [6.75, 0.360040],_
 [7.00, 0.381195],  [7.25, 0.408982],  [7.50, 0.441485],  [7.75, 0.476568],_
 [8.00, 0.512010],  [8.25, 0.545638],  [8.50, 0.575457],  [8.75, 0.599758],_
 [9.00, 0.617214],  [9.25, 0.626948],  [9.50, 0.628573],  [9.75, 0.622204],_
[10.00, 0.608436], [10.25, 0.588297], [10.50, 0.563176], [10.75, 0.534731],_
[11.00, 0.504784], [11.25, 0.475208], [11.50, 0.447809], [11.75, 0.424220],_
[12.00, 0.405810], [12.25, 0.393601], [12.50, 0.388217], [12.75, 0.389852],_
[13.00, 0.398268], [13.25, 0.412817], [13.50, 0.432489], [13.75, 0.455978],_
[14.00, 0.481770], [14.25, 0.508236], [14.50, 0.533736], [14.75, 0.556716],_
[15.00, 0.575803], [15.25, 0.589887], [15.50, 0.598183], [15.75, 0.600273],_
[16.00, 0.596126], [16.25, 0.586095], [16.50, 0.570890], [16.75, 0.551526],_
[17.00, 0.529259], [17.25, 0.505505], [17.50, 0.481750], [17.75, 0.459460],_
[18.00, 0.439989], [18.25, 0.424500], [18.50, 0.413893], [18.75, 0.408757],_
[19.00, 0.409336], [19.25, 0.415520], [19.50, 0.426853], [19.75, 0.442571],_
[20.00, 0.461646], [20.25, 0.482860], [20.50, 0.504875], [20.75, 0.526323],_
[21.00, 0.545885], [21.25, 0.562375], [21.50, 0.574811], [21.75, 0.582472],_
[22.00, 0.584939], [22.25, 0.582119], [22.50, 0.574246], [22.75, 0.561862],_
[23.00, 0.545782], [23.25, 0.527040], [23.50, 0.506824], [23.75, 0.486399],_
[24.00, 0.467029], [24.25, 0.449901], [24.50, 0.436051], [24.75, 0.426303],_
[25.00, 0.421217], [25.25, 0.421062], [25.50, 0.425797], [25.75, 0.435083],_
[26.00, 0.448300], [26.25, 0.464594], [26.50, 0.482927], [26.75, 0.502146],_
[27.00, 0.521054], [27.25, 0.538483], [27.50, 0.553369], [27.75, 0.564814],_
[28.00, 0.572142], [28.25, 0.574935], [28.50, 0.573060], [28.75, 0.566674],_
[29.00, 0.556212], [29.25, 0.542357], [29.50, 0.525995], [29.75, 0.508160],_
[30.00, 0.489969], [30.25, 0.472549], [30.50, 0.456974], [30.75, 0.444193],_
[31.00, 0.434973], [31.25, 0.429857], [31.50, 0.429129], [31.75, 0.432799],_
[32.00, 0.440605], [32.25, 0.452031], [32.50, 0.466343], [32.75, 0.482632],_
[33.00, 0.499873], [33.25, 0.516992], [33.50, 0.532930], [33.75, 0.546708],_
[34.00, 0.557490], [34.25, 0.564629], [34.50, 0.567709], [34.75, 0.566570],_
[35.00, 0.561313], [35.25, 0.552293], [35.50, 0.540094], [35.75, 0.525495],_
[36.00, 0.509417], [36.25, 0.492866], [36.50, 0.476871], [36.75, 0.462420],_
[37.00, 0.450396], [37.25, 0.441528], [37.50, 0.436345], [37.75, 0.435144],_
[38.00, 0.437971], [38.25, 0.444626], [38.50, 0.454670], [38.75, 0.467461],_
[39.00, 0.482187], [39.25, 0.497924], [39.50, 0.513690], [39.75, 0.528507],_
[40.00, 0.541464], [40.25, 0.551768], [40.50, 0.558799], [40.75, 0.562140],_
[41.00, 0.561608], [41.25, 0.557258], [41.50, 0.549384], [41.75, 0.538494],_
[42.00, 0.525282], [42.25, 0.510580], [42.50, 0.495309], [42.75, 0.480418],_
[43.00, 0.466829], [43.25, 0.455375], [43.50, 0.446755], [43.75, 0.441487],_
[44.00, 0.439878], [44.25, 0.442007], [44.50, 0.447720], [44.75, 0.456645],_
[45.00, 0.468209], [45.25, 0.481681], [45.50, 0.496215], [45.75, 0.510904],_
[46.00, 0.524837], [46.25, 0.537153], [46.50, 0.547099], [46.75, 0.554070],_
[47.00, 0.557650], [47.25, 0.557635], [47.50, 0.554044], [47.75, 0.547120],_
[48.00, 0.537309], [48.25, 0.525234], [48.50, 0.511657], [48.75, 0.497426],_
[49.00, 0.483428], [49.25, 0.470529], [49.50, 0.459523], [49.75, 0.451084],_
[50.00, 0.445722]]
--E
--S

[[x.1,x.2,fS(x.1),fS(x.1)-x.2] for x in perceyS]
--E
)spool
)lisp (bye)

———-

References^{[1] nothing}

7