|
DataMuseum.dkPresents historical artifacts from the history of: RC4000/8000/9000 |
This is an automatic "excavation" of a thematic subset of
See our Wiki for more about RC4000/8000/9000 Excavated with: AutoArchaeologist - Free & Open Source Software. |
top - metrics - download
Length: 36864 (0x9000)
Types: TextFile
Names: »elisout«
└─⟦621cfb9a2⟧ Bits:30002817 RC8000 Dump tape fra HCØ. Detaljer om "HC8000" projekt.
└─⟦0364f57e3⟧
└─⟦this⟧ »elisout«
└─⟦667bb35d6⟧ Bits:30007480 RC8000 Dump tape fra HCØ.
└─⟦4334b4c0b⟧
└─⟦this⟧ »elisout«
elis4 d.820428.1634
1
1 begin
2 real array field raf;
3 real procedure h1(z) ;
4 value z ; real z ;
5 begin
6 real b1,b2 ;
7 if z=-2.15 then b1:=0.8496786712 else
8 if z=-6.65 then b1:=-0.8496786712 else
9 b1:=0.6377458471*(1/(0.5625+(4.4+z)**2)-0.177777778)*
10 ln(abs((6.65+z)/(2.15+z)))+
11 2.124196678*(4.4+z)/(0.5625+(4.4+z)**2) ;
12 if z=-7 then b2:=1.274518007 else
13 if z=-10 then b2:=-1.274518007 else
14 b2:=0.4251639*(1/(0.25+(8.5+z)**2)-0.4)*ln(abs((10+z)/(7+z)))+
15 2.124196678*(8.5+z)/(0.25+(8.5+z)**2) ;
16 h1:=b1+b2 ;
17 end ;
18 real procedure rez(x,d1,d2,z1,z2) ;
19 value x,d1,d2,z1,z2 ; real x,d1,d2,z1,z2 ;
20 begin
21 real h1,h2,theta,h3 ;
22 h1:=((z1*z1-z2*z2)*x*x-d1**2+d2**2) ;
23 h2:=2*(z1*z2*x*x-d1*d2) ;
24 h3:=h1*h1+h2*h2 ;
25 theta:=arg(h1,h2) ;
26 if theta<0 then theta:=(theta+2*pi)/2 else theta:=theta/2 ;
27 rez:=x*(z1*cos(theta)+z2*sin(theta))/sqrt(sqrt(h3)) ;
28 slut1:
29 end ;
30 real procedure red(x,d1,d2,z1,z2) ;
31 value x,d1,d2,z1,z2 ; real x,d1,d2,z1,z2 ;
32 begin
33 real h1,h2,theta,h3 ;
34 h1:=((z1*z1-z2*z2)*x*x-d1**2+d2**2) ;
35 h2:=2*(z1*z2*x*x-d1*d2) ;
36 h3:=h1*h1+h2*h2 ;
37 theta:=arg(h1,h2) ;
38 if theta<0 then theta:=(theta+2*pi)/2 else theta:=theta/2 ;
39 red:=(d1*cos(theta)+d2*sin(theta))/(sqrt(sqrt(h3))) ;
40 slut2:
41 end ;
42 real procedure alfa(x) ;
43 value x ; real x ;
44 begin
45 real z ; z:=abs(x) ;
46 if z<=2.15 or z>=10 then alfa:=0
47 else
48 if z<6.65 then alfa:=sign(x)*0.63774595*
49 (1/((z-4.4)*(z-4.4)+0.5625)-0.177777778)
50 else
51 if z<=7 then alfa:=0
52 else
53 alfa:=sign(x)*0.42516392*
54 (1/((z-8.5)*(z-8.5)+0.25)-0.4) ;
55 end ;
56 integer i,j,k,l1,l2,l3,co,m ; real t,a,hj,d1hj,d2hj,a1,a2,a3,a4,res,hres ;
57 array int(-2304:4608),d1(0:2304),d2(0:2304),z1(0:2304),z2(0:2304)
58 ,hj1(0:4608),ada(0:2304),rezm(0:2304),redm(0:2304),alf(49:235),kff(0:6)
59 ,help(0:2304) ;
60 zone zo(1152,1,stderror) ;
61 t:=0.37332 ; a:=100/2304 ; co:=2304 ;
62 k:=1 ;
63 kff(0):=(41/140)*a ; kff(1):=(216/140)*a ; kff(2):=(27/140)*a ;
64 kff(3):=(272/140)*a ; kff(4):=kff(2) ; kff(5):=kff(1) ;
65 kff(6):=2*kff(0) ;
66 ada(0):=0 ;
67 for j:=1 step 1 until 1152 do ada(j):=adapint(50,70,h1(a*j+hj)-h1(hj-a*j)
68 ,hj,1.0'-3)+adapint(70,90,h1(a*j+hj)-h1(hj-a*j),hj,1.0'-3)
69 +adapint(90,110,h1(a*j+hj)-h1(hj-a*j),hj,1.0'-3) ;
70 for j:=49 step 6 until 229 do
71 for m:=1 step 1 until 6 do
72 begin
73 alf(j+m):=kff(m)*alfa((j+m)*a) ;
74 end ;
75
75 raf:= -4;
76 for i:=-co step 1 until 4608 do
77 int(i):=h1(i*a) ;
78 for j:=0 step 1 until 4608 do
79 begin
80 if j<2304 then hj1(j):=exp(j*a/t)+1
81 else
82 hj1(j):=1.0'100 ;
83 end ;
84 if k=1 then
85 begin
86 open(zo,4,<:gs3733:>,0) ;
87 inrec(zo,1152) ;
88 for j:=1 step 1 until 1152 do d1(j-1):=zo(j) ;
89 inrec(zo,1152) ;
90 for j:=1 step 1 until 1152 do d1(j+1151):=zo(j) ;
91 d1(2304):=d1(2303) ;
92 inrec(zo,1152) ;
93 for j:=1 step 1 until 1152 do d2(j-1):=zo(j) ;
94 inrec(zo,1152) ;
95 for j:=1 step 1 until 1152 do d2(j+1151):=zo(j) ;
96 d2(2304):=d2(2303) ;
97 inrec(zo,1152) ;
98 for j:=1 step 1 until 1152 do z1(j-1):=zo(j) ;
99 inrec(zo,1152) ;
100 for j:=1 step 1 until 1152 do z1(j+1151):=zo(j) ;
101 z1(2304):=z1(2303) ;
102 inrec(zo,1152) ;
103 for j:=1 step 1 until 1152 do z2(j-1):=zo(j) ;
104 inrec(zo,1152) ;
105 for j:=1 step 1 until 1152 do z2(j+1151):=zo(j) ;
106 z2(2304):=z2(2303) ;
107 close(zo,true) ;
108 end
109 else
110 begin
111 d1(0):= d2(0):= z2(0):= 0;
112 tofrom(d2,d2.raf,4*co);
113 tofrom(z2,z2.raf,4*co);
114 for i:=0 step 1 until co do
115 begin
116 if i<300 then
117 begin
118 d1(i):=2.68 ; z1(i):=2 ;
119 end
120 else
121 z1(i):=1 ; d1(i):=0 ;
122 end ;
123 end ;
124 a1:=(d1(30)*z1(30)+d2(30)*z2(30))/(z1(30)**2+z2(30)**2) ;
125 a2:=1.25/a1 ;
126 for i:=0 step 1 until co do
127 begin
128 d1(i):=a2*d1(i) ; d2(i):=a2*d2(i) ;
129 end ;
130 for i:=1 step 1 until 6 do
131 begin
132 for j:=0 step 1 until co do
133 begin
134 if j=0 then a1:=1.0'-2 else a1:=a*j ;
135 rezm(j):=rez(a1,d1(j),d2(j),z1(j),z2(j)) ;
136
136 redm(j):=red(a1,d1(j),d2(j),z1(j),z2(j)) ;
137 end ;
138 for k:=0 step 6 until co-6 do
139 for m:=1 step 1 until 6 do
140 begin
141 help(k+m):=kff(m)*redm(k+m) ;
142 end ;
143 for j:=0 step 1 until co do
144 begin
145 res:=0 ;
146 for k:=0 step 6 until co-6 do
147 for m:=1 step 1 until 6 do
148 begin
149 l3:=k+m ;
150 l1:=l3+j ; l2:=l3-j ;a2:=int(l1)+int(l2) ; a4:=hj1(l3) ;
151 if l1>co then a3:=int(-co) else a3:=int(-l1) ;
152 res:=res+(a2-
153 (0.11*(a4-2)+int(-l2)+a3+a2)/(a4))*help(l3) ;
154 end ;
155 d1(j):=res-kff(0)*(int(co+j)+int(co-j)-0.11)*redm(co) ;
156 res:=0 ;
157 for k:=49 step 6 until 229 do
158 for m:=1 step 1 until 6 do
159 begin
160 l3:=k+m ;
161 l1:=j+l3 ; l2:=l3-j ;
162 if l1>co then a1:=redm(co) else a1:=redm(l1) ;
163 if l2>=0 then
164 a2:=hj1(l2) else a2:=(1/(hj1(-l2)-1))+1 ;
165 if -l2<0 then a3:=-redm(l2) else a3:=redm(-l2) ;
166 res:=res+((a1+a3)/(hj1(l3)-2)+a1/(hj1(l1))+a3/(a2))*alf(l3) ;
167 end ;
168 d2(j):=pi*res ;
169 end ;
170 for k:=0 step 6 until co-6 do
171 for m:=1 step 1 until 6 do
172 begin
173 help(k+m):=kff(m)*rezm(k+m) ;
174 end ;
175 res:=0 ;
176 for k:=0 step 6 until co-6 do
177 for m:=1 step 1 until 6 do
178 begin
179 l3:=k+m ;
180 a1:=l3*a+1.0'-2 ; a2:=l3*a-1.0'-2 ;
181 res:=res+(h1(a1)-h1(a2)-
182 (h1(-a1)-h1(-a2)+h1(a1)-
183 h1(a2))/hj1(l3))*help(l3) ;
184 end ;
185 a4:=1.0'-2 ;
186 hres:=kff(0)*(h1(a4)-h1(-a4))*rezm(0) ;
187 z1(0):=1-1.0'2*(res+hres) ;
188 for j:=0 step 1 until co do
189 begin
190 if j=0 then goto nexte ;
191 if j>=1153 then
192 begin
193 z1(j):=1 ; goto nexte ;
194 end ;
195 res:=0 ;
196 for k:=0 step 6 until 1146 do
197 for m:=1 step 1 until 6 do
198 begin
199 l3:=k+m ;
200 l1:=l3+j ; l2:=l3-j ; a3:=int(l1)-int(l2) ;
201
201 if l1>co then a4:=int(-co) else a4:=int(-l1) ;
202 res:=res+(a3-
203 (a4-int(-l2)+a3)/(hj1(l3)))*help(l3) ;
204 end ;
205 hres:=kff(0)*((int(j)-int(-j))*rezm(0)-
206 (int(1152+j)-int(1152-j))*rezm(1152)) ;
207 z1(j):=1-((res+hres+ada(j))/(j*a)) ;
208 nexte:res:=0 ;
209 for k:=49 step 6 until 229 do
210 for m:=1 step 1 until 6 do
211 begin
212 l3:=k+m ;
213 l1:=l3+j ; l2:=l3-j ;
214 if l2>=0 then a1:=hj1(l2) else a1:=(1/(hj1(-l2)-1))+1 ;
215 if l1>co then a2:=rezm(co) else a2:=rezm(l1) ;
216 res:=res+((
217 rezm(abs(l2))+a2)/(hj1(l3)-2)
218 +rezm(abs(l2))/(a1)+a2/(hj1(l1)))*alf(l3) ;
219 end ;
220 if j=0 then z2(0):=1.0'2*pi*res else
221 z2(j):=(pi/(j*a))*res ;
222 end ;
223 write(out,<:<10>:>,z1(0),z2(0),d1(0),d2(0),<:<10>:>) ;
224 for j:=0 step 1 until co do
225 begin
226 if j<250 or j>2295 then
227 begin
228 hj:=z1(j)**2+z2(j)**2 ;
229 write(out,<:<10>:>,a*j,z1(j),z2(j),
230 (d1(j)*z1(j)+d2(j)*z2(j))/hj,
231 (d2(j)*z1(j)-d1(j)*z2(j))/hj) ;
232 end ;
233 end ;
234 a1:=(d1(30)*z1(30)+d2(30)*z2(30))/(z1(30)**2+z2(30)**2) ;
235 if a1<=1.37 then goto slut ;
236 end ;
237 slut:
238 begin
239 open(zo,4,<:gs3733:>,0) ;
240 outrec(zo,1152) ;
241 for i:=1 step 1 until 1152 do zo(i):=d1(i-1) ;
242 outrec(zo,1152) ;
243 for i:=1 step 1 until 1152 do zo(i):=d1(i+1151) ;
244 outrec(zo,1152) ;
245 for i:=1 step 1 until 1152 do zo(i):=d2(i-1) ;
246 outrec(zo,1152) ;
247 for i:=1 step 1 until 1152 do zo(i):=d2(i+1151) ;
248 outrec(zo,1152) ;
249 for i:=1 step 1 until 1152 do zo(i):=z1(i-1) ;
250 outrec(zo,1152) ;
251 for i:=1 step 1 until 1152 do zo(i):=z1(i+1151) ;
252 outrec(zo,1152) ;
253 for i:=1 step 1 until 1152 do zo(i):=z2(i-1) ;
254 outrec(zo,1152) ;
255 for i:=1 step 1 until 1152 do zo(i):=z2(i+1151) ;
256 close(zo,true) ;
257 end ;
258 end
algol end 47
2.0362 0.1053 2.4238 -0.0000
0.0000 2.0362 0.1053 1.1872 -0.0614
0.0434 2.0414 0.0244 1.1872 -0.0142
0.0868 2.0415 0.0123 1.1874 -0.0072
0.1302 2.0417 0.0084 1.1875 -0.0049
0.1736 2.0419 0.0064 1.1877 -0.0038
0.2170 2.0421 0.0053 1.1879 -0.0032
0.2604 2.0424 0.0046 1.1882 -0.0028
0.3038 2.0428 0.0041 1.1884 -0.0025
0.3472 2.0433 0.0038 1.1888 -0.0024
0.3906 2.0438 0.0036 1.1892 -0.0023
0.4340 2.0443 0.0034 1.1896 -0.0022
0.4774 2.0450 0.0034 1.1900 -0.0022
0.5208 2.0457 0.0033 1.1906 -0.0022
0.5642 2.0464 0.0033 1.1911 -0.0023
0.6076 2.0472 0.0033 1.1917 -0.0024
0.6510 2.0481 0.0034 1.1923 -0.0025
0.6944 2.0491 0.0035 1.1930 -0.0026
0.7378 2.0501 0.0037 1.1937 -0.0028
0.7813 2.0512 0.0038 1.1944 -0.0030
0.8247 2.0524 0.0041 1.1951 -0.0033
0.8681 2.0536 0.0044 1.1959 -0.0036
0.9115 2.0550 0.0047 1.1967 -0.0040
0.9549 2.0562 0.0051 1.1977 -0.0044
0.9983 2.0575 0.0058 1.1988 -0.0053
1.0417 2.0587 0.0058 1.2000 -0.0053
1.0851 2.0600 0.0067 1.2012 -0.0064
1.1285 2.0614 0.0068 1.2026 -0.0066
1.1719 2.0628 0.0077 1.2041 -0.0076
1.2153 2.0642 0.0074 1.2056 -0.0075
1.2587 2.0656 0.0089 1.2072 -0.0093
1.3021 2.0672 0.0082 1.2089 -0.0085
1.3455 2.0687 0.0097 1.2107 -0.0105
1.3889 2.0703 0.0095 1.2126 -0.0103
1.4323 2.0720 0.0107 1.2145 -0.0120
1.4757 2.0737 0.0101 1.2166 -0.0113
1.5191 2.0755 0.0124 1.2187 -0.0145
1.5625 2.0773 0.0109 1.2209 -0.0126
1.6059 2.0792 0.0132 1.2232 -0.0159
1.6493 2.0811 0.0126 1.2257 -0.0152
1.6927 2.0831 0.0144 1.2282 -0.0178
1.7361 2.0851 0.0132 1.2309 -0.0164
1.7795 2.0872 0.0166 1.2336 -0.0213
1.8229 2.0893 0.0143 1.2365 -0.0181
1.8663 2.0915 0.0175 1.2395 -0.0231
1.9097 2.0937 0.0164 1.2427 -0.0217
1.9531 2.0960 0.0189 1.2459 -0.0257
1.9965 2.0983 0.0171 1.2494 -0.0232
2.0399 2.1006 0.0217 1.2530 -0.0306
2.0833 2.1030 0.0183 1.2569 -0.0255
2.1267 2.1055 0.0227 1.2608 -0.0328
2.1701 2.1079 0.0211 1.2651 -0.0305
2.2135 2.1104 0.0244 1.2694 -0.0362
2.2569 2.1130 0.0217 1.2742 -0.0320
2.3003 2.1155 0.0279 1.2789 -0.0429
2.3438 2.1181 0.0229 1.2844 -0.0348
2.3872 2.1207 0.0288 1.2897 -0.0453
2.4306 2.1233 0.0262 1.2956 -0.0413
2.4740 2.1259 0.0304 1.3016 -0.0492
2.5174 2.1286 0.0263 1.3084 -0.0424
2.5608 2.1313 0.0345 1.3149 -0.0577
2.6042 2.1339 0.0271 1.3227 -0.0449
2.6476 2.1366 0.0346 1.3300 -0.0593
2.6910 2.1394 0.0305 1.3385 -0.0523
2.7344 2.1421 0.0356 1.3468 -0.0626
2.7778 2.1450 0.0294 1.3563 -0.0514
2.8212 2.1479 0.0392 1.3652 -0.0713
2.8646 2.1509 0.0290 1.3758 -0.0520
2.9080 2.1540 0.0375 1.3857 -0.0696
2.9514 2.1573 0.0315 1.3969 -0.0585
2.9948 2.1608 0.0365 1.4077 -0.0695
3.0382 2.1645 0.0283 1.4196 -0.0535
3.0816 2.1686 0.0383 1.4307 -0.0751
3.1250 2.1729 0.0262 1.4433 -0.0505
3.1684 2.1777 0.0341 1.4550 -0.0679
3.2118 2.1830 0.0271 1.4676 -0.0538
3.2552 2.1888 0.0310 1.4797 -0.0629
3.2986 2.1952 0.0226 1.4926 -0.0453
3.3420 2.2023 0.0307 1.5049 -0.0637
3.3854 2.2103 0.0198 1.5185 -0.0398
3.4288 2.2198 0.0264 1.5320 -0.0527
3.4722 2.2294 0.0207 1.5458 -0.0395
3.5156 2.2393 0.0277 1.5590 -0.0417
3.5590 2.2492 0.0200 1.5724 -0.0284
3.6024 2.2597 0.0316 1.5853 -0.0359
3.6458 2.2702 0.0239 1.5986 -0.0181
3.6892 2.2810 0.0353 1.6115 -0.0203
3.7326 2.2920 0.0278 1.6247 -0.0143
3.7760 2.3033 0.0441 1.6378 -0.0047
3.8194 2.3148 0.0332 1.6509 -0.0015
3.8628 2.3265 0.0507 1.6642 0.0032
3.9063 2.3383 0.0440 1.6776 0.0126
3.9497 2.3505 0.0612 1.6913 0.0206
3.9931 2.3627 0.0500 1.7047 0.0175
4.0365 2.3752 0.0782 1.7193 0.0409
4.0799 2.3878 0.0617 1.7326 0.0325
4.1233 2.4007 0.0879 1.7477 0.0502
4.1667 2.4136 0.0806 1.7617 0.0518
4.2101 2.4268 0.1062 1.7774 0.0712
4.2535 2.4401 0.0890 1.7909 0.0577
4.2969 2.4536 0.1335 1.8088 0.0988
4.3403 2.4671 0.1084 1.8216 0.0767
4.3837 2.4808 0.1473 1.8397 0.1100
4.4271 2.4945 0.1377 1.8539 0.1039
4.4705 2.5083 0.1756 1.8729 0.1374
4.5139 2.5222 0.1496 1.8850 0.1110
4.5573 2.5360 0.2168 1.9086 0.1762
4.6007 2.5498 0.1791 1.9183 0.1369
4.6441 2.5635 0.2365 1.9413 0.1907
4.6875 2.5770 0.2228 1.9538 0.1761
4.7309 2.5903 0.2787 1.9777 0.2286
4.7743 2.6032 0.2397 1.9845 0.1853
4.8177 2.6157 0.3395 2.0182 0.2845
4.8611 2.6275 0.2833 2.0184 0.2222
4.9045 2.6388 0.3672 2.0494 0.3044
4.9479 2.6491 0.3472 2.0553 0.2797
4.9913 2.6585 0.4281 2.0868 0.3585
5.0347 2.6666 0.3705 2.0791 0.2931
5.0781 2.6732 0.5143 2.1304 0.4388
5.1215 2.6781 0.4322 2.1078 0.3475
5.1649 2.6811 0.5504 2.1512 0.4675
5.2083 2.6818 0.5203 2.1400 0.4316
5.2517 2.6801 0.6308 2.1825 0.5441
5.2951 2.6755 0.5485 2.1423 0.4534
5.3385 2.6680 0.7398 2.2221 0.6544
5.3819 2.6573 0.6260 2.1534 0.5327
5.4253 2.6434 0.7743 2.2140 0.6937
5.4688 2.6260 0.7298 2.1695 0.6489
5.5122 2.6052 0.8584 2.2239 0.7924
5.5556 2.5812 0.7518 2.1302 0.6826
5.5990 2.5543 0.9629 2.2432 0.9236
5.6424 2.5246 0.8255 2.1107 0.7851
5.6858 2.4928 0.9727 2.1855 0.9667
5.7292 2.4590 0.9146 2.1002 0.9187
5.7726 2.4241 1.0280 2.1579 1.0679
5.8160 2.3883 0.9145 2.0146 0.9557
5.8594 2.3522 1.0914 2.1413 1.1879
5.9028 2.3162 0.9589 1.9666 1.0543
5.9462 2.2808 1.0684 2.0383 1.2148
5.9896 2.2463 1.0087 1.9324 1.1672
6.0330 2.2130 1.0842 1.9789 1.2880
6.0764 2.1809 0.9875 1.8256 1.1861
6.1198 2.1504 1.1050 1.9337 1.3679
6.1632 2.1215 1.0008 1.7651 1.2509
6.2066 2.0942 1.0674 1.8178 1.3658
6.2500 2.0684 1.0176 1.7199 1.3193
6.2934 2.0443 1.0601 1.7484 1.3992
6.3368 2.0216 0.9891 1.6192 1.3112
6.3802 2.0004 1.0584 1.6949 1.4351
6.4236 1.9806 0.9849 1.5615 1.3375
6.4670 1.9622 1.0222 1.5933 1.4068
6.5104 1.9449 0.9844 1.5180 1.3647
6.5538 1.9287 1.0070 1.5321 1.4075
6.5972 1.9136 0.9575 1.4380 1.3348
6.6406 1.8995 0.9967 1.4855 1.4096
6.6840 1.8863 0.9463 1.3930 1.3326
6.7274 1.8739 0.9673 1.4096 1.3633
6.7708 1.8622 0.9382 1.3605 1.3309
6.8142 1.8512 0.9504 1.3662 1.3425
6.8576 1.8408 0.9157 1.3066 1.2869
6.9010 1.8310 0.9373 1.3358 1.3226
6.9444 1.8216 0.9015 1.2814 1.2685
6.9878 1.8127 0.9141 1.2897 1.2665
7.0313 1.8043 0.8896 1.2674 1.2505
7.0747 1.7963 0.8965 1.2705 1.2373
7.1181 1.7888 0.8703 1.2421 1.2018
7.1615 1.7817 0.8810 1.2618 1.2099
7.2049 1.7751 0.8537 1.2391 1.1807
7.2483 1.7691 0.8606 1.2439 1.1585
7.2917 1.7636 0.8388 1.2418 1.1597
7.3351 1.7588 0.8410 1.2442 1.1360
7.3785 1.7546 0.8210 1.2352 1.1182
7.4219 1.7510 0.8235 1.2471 1.1134
7.4653 1.7481 0.8037 1.2408 1.0995
7.5087 1.7458 0.8046 1.2429 1.0750
7.5521 1.7440 0.7883 1.2467 1.0772
7.5955 1.7428 0.7862 1.2484 1.0543
7.6389 1.7422 0.7727 1.2485 1.0419
7.6823 1.7420 0.7697 1.2554 1.0290
7.7257 1.7424 0.7572 1.2584 1.0181
7.7691 1.7433 0.7539 1.2620 0.9939
7.8125 1.7447 0.7431 1.2710 0.9892
7.8559 1.7466 0.7382 1.2768 0.9657
7.8993 1.7492 0.7299 1.2875 0.9539
7.9427 1.7524 0.7243 1.2983 0.9333
7.9861 1.7559 0.7170 1.3127 0.9232
8.0295 1.7598 0.7126 1.3249 0.9012
8.0729 1.7639 0.7056 1.3420 0.8925
8.1163 1.7686 0.7012 1.3575 0.8727
8.1597 1.7737 0.6954 1.3765 0.8601
8.2031 1.7795 0.6909 1.3961 0.8423
8.2465 1.7861 0.6856 1.4200 0.8287
8.2899 1.7941 0.6826 1.4475 0.8109
8.3333 1.8019 0.6782 1.4760 0.7990
8.3767 1.8102 0.6777 1.5084 0.7906
8.4201 1.8186 0.6723 1.5350 0.7725
8.4635 1.8275 0.6768 1.5769 0.7811
8.5069 1.8366 0.6695 1.6011 0.7557
8.5503 1.8461 0.6754 1.6456 0.7680
8.5938 1.8559 0.6711 1.6754 0.7510
8.6372 1.8661 0.6791 1.7239 0.7693
8.6806 1.8766 0.6698 1.7464 0.7347
8.7240 1.8876 0.6893 1.8130 0.7880
8.7674 1.8988 0.6753 1.8288 0.7376
8.8108 1.9104 0.6945 1.8949 0.7879
8.8542 1.9222 0.6886 1.9235 0.7615
8.8976 1.9344 0.7102 1.9932 0.8165
8.9410 1.9468 0.6932 2.0052 0.7558
8.9844 1.9595 0.7376 2.1086 0.8764
9.0278 1.9722 0.7115 2.1051 0.7882
9.0712 1.9850 0.7532 2.2026 0.8974
9.1146 1.9976 0.7439 2.2219 0.8589
9.1580 2.0100 0.7891 2.3222 0.9728
9.2014 2.0218 0.7581 2.3035 0.8741
9.2448 2.0330 0.8463 2.4669 1.1026
9.2882 2.0428 0.7980 2.4104 0.9610
9.3316 2.0512 0.8787 2.5554 1.1653
9.3750 2.0573 0.8623 2.5372 1.1159
9.4184 2.0608 0.9462 2.6802 1.3226
9.4618 2.0608 0.8898 2.5744 1.1816
9.5052 2.0570 1.0440 2.8339 1.5620
9.5486 2.0483 0.9582 2.6425 1.3719
9.5920 2.0345 1.0866 2.8476 1.6976
9.6354 2.0151 1.0538 2.7266 1.6581
9.6788 1.9904 1.1688 2.9050 1.9579
9.7222 1.9604 1.0787 2.6274 1.8163
9.7656 1.9261 1.2625 2.9740 2.2810
9.8090 1.8884 1.1437 2.5904 2.1139
9.8524 1.8486 1.2642 2.8079 2.4600
9.8958 1.8078 1.2117 2.5719 2.4505
9.9392 1.7674 1.2942 2.7164 2.7172
9.9826 1.7279 1.2004 2.3449 2.6248
10.0260 1.6904 1.3188 2.6337 2.9674
10.0694 1.6550 1.2177 2.2241 2.8586
10.1128 1.6223 1.2790 2.3594 3.0797
10.1563 1.5919 1.2314 2.1248 3.0683
10.1997 1.5642 1.2653 2.1867 3.2191
10.2431 1.5387 1.2004 1.8839 3.1518
10.2865 1.5156 1.2531 2.0390 3.3373
10.3299 1.4943 1.1896 1.7375 3.2565
10.3733 1.4751 1.2121 1.7869 3.3648
10.4167 1.4574 1.1795 1.6126 3.3405
10.4601 1.4413 1.1889 1.6168 3.4073
10.5035 1.4265 1.1495 1.4129 3.3509
10.5469 1.4130 1.1688 1.4713 3.4374
10.5903 1.4005 1.1317 1.2772 3.3736
10.6337 1.3891 1.1365 1.2763 3.4151
10.6771 1.3785 1.1153 1.1589 3.3852
10.7205 1.3687 1.1141 1.1343 3.4029
10.7639 1.3596 1.0907 1.0073 3.3566
10.8073 1.3513 1.0940 1.0104 3.3838
99.6528 1.0000 0.1018 -0.2971 -0.0028
99.6962 1.0000 0.1017 -0.2970 -0.0028
99.7396 1.0000 0.1017 -0.2970 -0.0028
99.7830 1.0000 0.1016 -0.2969 -0.0028
99.8264 1.0000 0.1016 -0.2968 -0.0028
99.8698 1.0000 0.1016 -0.2967 -0.0028
99.9132 1.0000 0.1015 -0.2967 -0.0028
99.9566 1.0000 0.1015 -0.2966 -0.0028
10.0000'1 1.0000 0.1014 -0.2965 -0.0028
end 34
heatcvscd d.820427.2201
1
1 comment:
1 DETTE PROGRAM UDREGNER VARMELEDNINGEN FOR BLY ;
1
1 begin
2 real procedure simp(l,k,j,f) ;
3 value l,k,j ; integer k,j ; real l ; array f ;
4 begin
5 integer i ; real x,y,res ; res:=0 ;
6 x:=l/(3*k) ; y:=4*x ;
7 for i:=(j+1) step 2 until (j+k-1) do
8 begin
9 res:=res+y*f(i) ;
10 end ;
11 res:=res+(f(j)+f(j+k))*x ; y:=2*x ;
12 for i:=(j+2) step 2 until (j+k-2) do
13 begin
14 res:=res+y*f(i) ;
15 end ;
16 simp:=res ;
17 end ;
18 real procedure rez(x,d1,d2,z1,z2) ;
19 value x,d1,d2,z1,z2 ; real x,d1,d2,z1,z2 ;
20 begin
21 real h1,h2,h3,theta ;
22 h1:=((z1**2-z2**2)*(x**2)-d1**2+d2**2) ;
23 h2:=2*(z1*z2*x*x-d1*d2) ;
24 h3:=h1**2+h2**2 ;
25 theta:=arg(h1,h2) ;
26 if theta<0 then theta:=(theta+2*pi)/2 else theta:=theta/2 ;
27 rez:=x*(z1*cos(theta)+z2*sin(theta))/sqrt(sqrt(h3)) ;
28 end ;
29 real procedure red(x,d1,d2,z1,z2) ;
30 value x,d1,d2,z1,z2 ; real x,d1,d2,z1,z2 ;
31 begin
32 real h1,h2,h3,theta ;
33 h1:=((z1**2-z2**2)*(x**2)-d1**2+d2**2) ;
34 h2:=2*(z1*z2*(x**2)-d1*d2) ;
35 h3:=h1**2+h2**2 ;
36 theta:=arg(h1,h2) ;
37 if theta<0 then theta:=(theta+2*pi)/2 else theta:=theta/2 ;
38 red:=(d1*cos(theta)+d2*sin(theta))/sqrt(sqrt(h3)) ;
39 end ;
40 real procedure rn2(x,d1,d2,z1,z2) ;
41 value x,d1,d2,z1,z2 ; real x,d1,d2,z1,z2 ;
42 begin
43 real h1,h2,h3,theta ;
44 h1:=((z1**2-z2**2)*(x**2)-d1**2+d2**2) ;
45 h2:=2*(z1*z2*(x**2)-d1*d2) ;
46 h3:=h1**2+h2**2 ;
47 theta:=arg(h1,h2) ;
48 if theta<0 then theta:=(theta+2*pi)/2 else theta:=theta/2 ;
49 rn2:=-(d2*cos(theta)-d1*sin(theta))/sqrt(sqrt(h3)) ;
50 end ;
51 real procedure alf(x,y) ;
52 value x,y ; real x,y ;
53 begin
54 real z1,z2 ; z1:=abs(x) ; z2:=abs(y) ;
55 if z1<=0 or z1>=10 then alf:=0 else
56 if z1<6.65 then
57 begin
58 real hj1,hj2 ; hj1:=2.74905'-3*z1*z1/sinh(z2/2) ;
59 hj2:=0.63774595*
60 (1/((z1-4.4)*(z1-4.4)+0.5625)-0.1777777778)/sinh(z2/2) ;
61 if hj1>hj2 then alf:=hj1 else alf:=hj2 ;
62 end
63 else
64 if z1<=7 then alf:=0 else
65 alf:=0.42516392*
66 (1/((z1-8.5)*(z1-8.5)+0.25)-0.4)/sinh(z2/2) ;
67 end ;
68 real procedure cosh(x) ;
69 value x ; real x ;
70 begin
71 real a ; a:=exp(x) ; cosh:=(a+1/a)/2 ;
72 end ;
73 integer i,j,k,n,l ; real t,tauf,q,q1 ;
74 array kohy(0:510),f(0:510),e1(0:510),e(0:510),
75 x1(0:511),iz(1:510),ifi(1:510),n1(1:510),n2(1:510),
76 r2(1:510),nsqr(0:510) ;
77 zone zo(512,1,stderror),zo1(1152,1,stderror);
78 n:=510 ; t:=0.37332 ;
79 k:=entier(40*t*1152/50)+2 ;
80 nsqr(0):=0 ;
81 for i:=0 step 1 until n do
82 begin
83 e(i):=i*40*t/n ; e1(i):=e(i)/t ;
84 kohy(i):=cosh(e1(i)/2) ;
85 end ;
86 begin
87 array rz(1:n),rfi(1:n),e2(1:k),m(1:k),y(1:k),a(0:4) ;
88 a(0):=3 ; e2(1):=0.01 ;
89 for i:=2 step 1 until k do e2(i):=(i-1)*50/1152 ;
90 open(zo1,4,<:gs3733:>,0) ;
91 inrec(zo1,1152) ;
92 for i:=1 step 1 until k do y(i):=zo1(i) ;
93 spln3(a,e2,y,m,k) ;
94 for i:=1 step 1 until n do rfi(i):=spln3val(e(i),e2,y,m,k) ;
95 inrec(zo1,1152) ;
96 inrec(zo1,1152) ;
97 for i:=1 step 1 until k do y(i):=zo1(i) ;
98 spln3(a,e2,y,m,k) ;
99 for i:=1 step 1 until n do ifi(i):=spln3val(e(i),e2,y,m,k) ;
100 inrec(zo1,1152) ;
101 inrec(zo1,1152) ;
102 for i:=1 step 1 until k do y(i):=zo1(i) ;
103 spln3(a,e2,y,m,k) ;
104 for i:=1 step 1 until n do rz(i):=spln3val(e(i),e2,y,m,k) ;
105 inrec(zo1,1152) ;
106 inrec(zo1,1152) ;
107 for i:=1 step 1 until k do y(i):=zo1(i) ;
108 spln3(a,e2,y,m,k) ;
109 for i:=1 step 1 until n do iz(i):=spln3val(e(i),e2,y,m,k) ;
110 close(zo1,true) ;
111 for i:=1 step 1 until n do
112 begin
113 n1(i):=rez(e(i),rfi(i),ifi(i),rz(i),iz(i)) ;
114 n2(i):=rn2(e(i),rfi(i),ifi(i),rz(i),iz(i)) ;
115 r2(i):=red(e(i),rfi(i),ifi(i),rz(i),iz(i)) ;
116 nsqr(i):=n1(i)*n1(i)-r2(i)*r2(i) ;
117 end ;
118 end ;
119 for i:=1 step 1 until n do
120 begin
121 array hjm(0:n) ;
122 for j:=1 step 1 until n do
123 begin
124 real hj1,hj2 ; hj1:=n1(i)*n1(j) ; hj2:=r2(i)*r2(j) ;
125 hjm(j):=((hj1-hj2)*alf(e(i)-e(j),e1(i)-e1(j))+(hj1+hj2)*
126 alf(e(i)+e(j),e1(i)+e1(j)))/kohy(j) ;
127 end ;
128 hjm(0):=0 ;
129 f(i):=0.5*pi*kohy(i)*simp(40*t,n,0,hjm)/t ;
130 end ;
131 for tauf:=0 do
132 begin
133 array dr(0:510) ;
134 q:=0 ;
135 f(0):=dr(0):=0 ;
136 for i:=1 step 1 until n do
137 begin
138 f(i):=f(i)+nsqr(i)*tauf*0.0331726/t ;
139 dr(i):=-(e1(i)*nsqr(i))/(f(i)*kohy(i)) ;
140 end ;
141 open(zo,4,<:min99:>,0) ;
142 for i:=0 step 1 until n do
143 begin
144 outrec(zo,512) ;
145 if i=0 then
146 begin
147 cleararray(zo) ; goto nexte ;
148 end ;
149 for j:=1 step 2 until (n-1) do
150 begin
151 zo(j+1):=-2*e1(1)*pi*nsqr(i)*nsqr(j)*(1+n2(i)*n2(j)/(n1(i)*n1(j)))*
152 (alf(e(j)-e(i),e1(j)-e1(i))-alf(e(j)+e(i),e1(j)+e1(i)))/(3*f(i)) ;
153 end ;
154 for j:=2 step 2 until (n-2) do
155 begin
156 zo(j+1):=-e1(1)*pi*nsqr(i)*nsqr(j)*(1+n2(i)*n2(j)/(n1(i)*n1(j)))*
157 (alf(e(j)-e(i),e1(j)-e1(i))-alf(e(j)+e(i),e1(j)+e1(i)))/(3*f(i)) ;
158 end ;
159 zo(1):=0 ;
160 zo(n+1):=-e1(1)*pi*nsqr(i)*nsqr(n)*(1+n2(i)*n2(n)/(n1(i)*n1(n)))*
161 (alf(e(n)-e(i),e1(n)-e1(i))-alf(e(n)+e(i),e1(n)+e1(i)))/(6*f(i)) ;
162 nexte:zo(i+1):=1+zo(i+1) ; zo(512):=0 ;
163 dr(i):=dr(i)/zo(i+1) ;
164 for j:=1 step 1 until (n+1) do
165 begin
166 zo(j):=zo(j)/zo(i+1) ;
167 end ;
168 zo(i+1):=0 ;
169 end ;
170 cleararray(x1) ;
171 for i:=1 step 1 until 100 do
172 begin
173 array x2(0:511),x3(0:511) ;
174 l:=i ;
175 setposition(zo,0,0) ;
176 for j:=0 step 1 until n do
177 begin
178 inrec(zo,512) ;
179 x2(j):=dr(j)-prik(512,zo,x1) ;
180 x3(j):=e1(j)*nsqr(j)*x2(j)/kohy(j) ;
181 x1(j):=x2(j) ;
182 end ;
183 q1:=0.0759908877*simp(40,n,0,x3) ;
184 write(out,<<dddd>,l,<<-d.dddd'-dd>,t/0.6222,<: :>,q1,<:<10>:>) ;
185 if abs(abs(q/q1)-1)<0.0001 then goto slut ;
186 q:=q1 ;
187 end ;
188 slut:close(zo,true) ;
189 write(out,<<dddd>,<:<10><10>totale antal punkter=:>,n) ;
190 write(out,<<-d.dddd'-dd>,
191 <: temperatur=:>,t,<: tauf=:>,tauf,
192 <:<10><10> T*TAUk=:>,-q1,<: antal itterationer=:>,l,<:<10><10>:>) ;
193 for j:=0 step 1 until n do
194 write(out,<<-d.dddd'-dd>,-x1(j)*kohy(j),
195 if j mod 5=0 then <:<10>:> else <: :>) ;
196 end ;
197 end
algol end 62
1 6.0000' -1 -1.1789
2 6.0000' -1 -1.1084
3 6.0000' -1 -1.1105
4 6.0000' -1 -1.1105
totale antal punkter= 510 temperatur= 3.7332' -1 tauf= 0.0000
T*TAUk= 1.1105 antal itterationer= 4.0000
0.0000
1.6022' -3 -1.6561' -3 1.0899' -3 2.7963' -3 9.3556' -4
2.7654' -3 3.8860' -3 4.6040' -3 6.2888' -3 8.1125' -3
1.0259' -2 1.2930' -2 1.6101' -2 1.9853' -2 2.4302' -2
2.9509' -2 3.5607' -2 4.2662' -2 5.0845' -2 6.0296' -2
7.1221' -2 8.3769' -2 9.8205' -2 1.1476' -1 1.3383' -1
1.5570' -1 1.8120' -1 2.1065' -1 2.4424' -1 2.8764' -1
3.3914' -1 3.8725' -1 4.6798' -1 5.7015' -1 6.1517' -1
7.1643' -1 9.0155' -1 1.0247 1.2950 1.9025
4.0229 1.3233' 1 1.8349' 1 2.1538' 1 2.3729' 1
2.5289' 1 2.6381' 1 2.7156' 1 2.7662' 1 2.7978' 1
2.8121' 1 2.8149' 1 2.8053' 1 2.7885' 1 2.7631' 1
2.7329' 1 2.6964' 1 2.6572' 1 2.6132' 1 2.5679' 1
2.5192' 1 2.4697' 1 2.4180' 1 2.3663' 1 2.3131' 1
2.2604' 1 2.2069' 1 2.1540' 1 2.1011' 1 2.0491' 1
1.9974' 1 1.9468' 1 1.8970' 1 1.8484' 1 1.8009' 1
1.7548' 1 1.7100' 1 1.6668' 1 1.6252' 1 1.5851' 1
1.5469' 1 1.5105' 1 1.4759' 1 1.4434' 1 1.4128' 1
1.3844' 1 1.3582' 1 1.3341' 1 1.3124' 1 1.2930' 1
1.2761' 1 1.2617' 1 1.2500' 1 1.2408' 1 1.2345' 1
1.2310' 1 1.2304' 1 1.2328' 1 1.2384' 1 1.2472' 1
1.2593' 1 1.2749' 1 1.2939' 1 1.3166' 1 1.3431' 1
1.3733' 1 1.4074' 1 1.4455' 1 1.4877' 1 1.5339' 1
1.5842' 1 1.6385' 1 1.6969' 1 1.7591' 1 1.8251' 1
1.8946' 1 1.9673' 1 2.0431' 1 2.1213' 1 2.1958' 1
2.2634' 1 2.3242' 1 2.3790' 1 2.4278' 1 2.4706' 1
2.5072' 1 2.5375' 1 2.5613' 1 2.5789' 1 2.5902' 1
2.5957' 1 2.5956' 1 2.5905' 1 2.5807' 1 2.5670' 1
2.5497' 1 2.5294' 1 2.5065' 1 2.4818' 1 2.4553' 1
2.4277' 1 2.3992' 1 2.3701' 1 2.3407' 1 2.3115' 1
2.2819' 1 2.2528' 1 2.2241' 1 2.1959' 1 2.1680' 1
2.1411' 1 2.1144' 1 2.0888' 1 2.0637' 1 2.0393' 1
2.0155' 1 1.9928' 1 1.9704' 1 1.9491' 1 1.9281' 1
1.9081' 1 1.8884' 1 1.8696' 1 1.8512' 1 1.8341' 1
1.8172' 1 1.8019' 1 1.7870' 1 1.7730' 1 1.7599' 1
1.7475' 1 1.7352' 1 1.7252' 1 1.7147' 1 1.7046' 1
1.6963' 1 1.6885' 1 1.6804' 1 1.6750' 1 1.6681' 1
1.6623' 1 1.6598' 1 1.6558' 1 1.6508' 1 1.6511' 1
1.6492' 1 1.6480' 1 1.6487' 1 1.6458' 1 1.6442' 1
1.6488' 1 1.6473' 1 1.6456' 1 1.6505' 1 1.6509' 1
1.6517' 1 1.6567' 1 1.6552' 1 1.6564' 1 1.6639' 1
1.6623' 1 1.6616' 1 1.6671' 1 1.6665' 1 1.6669' 1
1.6702' 1 1.6666' 1 1.6665' 1 1.6704' 1 1.6656' 1
1.6629' 1 1.6640' 1 1.6597' 1 1.6569' 1 1.6554' 1
1.6488' 1 1.6456' 1 1.6439' 1 1.6365' 1 1.6315' 1
1.6282' 1 1.6215' 1 1.6165' 1 1.6119' 1 1.6046' 1
1.6001' 1 1.5959' 1 1.5890' 1 1.5843' 1 1.5805' 1
1.5750' 1 1.5714' 1 1.5682' 1 1.5640' 1 1.5621' 1
1.5606' 1 1.5582' 1 1.5582' 1 1.5587' 1 1.5590' 1
1.5614' 1 1.5644' 1 1.5676' 1 1.5729' 1 1.5787' 1
1.5840' 1 1.5912' 1 1.5986' 1 1.6062' 1 1.6157' 1
1.6258' 1 1.6362' 1 1.6482' 1 1.6606' 1 1.6734' 1
1.6876' 1 1.7020' 1 1.7165' 1 1.7321' 1 1.7478' 1
1.7636' 1 1.7800' 1 1.7964' 1 1.8127' 1 1.8296' 1
1.8461' 1 1.8629' 1 1.8799' 1 1.9120' 1 1.9364' 1
1.9623' 1 1.9810' 1 2.0009' 1 2.0170' 1 2.0335' 1
2.0475' 1 2.0615' 1 2.0736' 1 2.0853' 1 2.0955' 1
2.1043' 1 2.1114' 1 2.1171' 1 2.1213' 1 2.1243' 1
2.1261' 1 2.1265' 1 2.1259' 1 2.1244' 1 2.1214' 1
2.1176' 1 2.1129' 1 2.1072' 1 2.1004' 1 2.0933' 1
2.0843' 1 2.0756' 1 2.0660' 1 2.0550' 1 2.0436' 1
2.0324' 1 2.0192' 1 2.0067' 1 1.9933' 1 1.9785' 1
1.9646' 1 1.9510' 1 1.9340' 1 1.9193' 1 1.9042' 1
1.8875' 1 1.8719' 1 1.8562' 1 1.8373' 1 1.8232' 1
1.8078' 1 1.7886' 1 1.7735' 1 1.7593' 1 1.7413' 1
1.7287' 1 1.7139' 1 1.6962' 1 1.6877' 1 1.6775' 1
1.6603' 1 1.6540' 1 1.6472' 1 1.6365' 1 1.6348' 1
1.6296' 1 1.6200' 1 1.6269' 1 1.6275' 1 1.6193' 1
1.6263' 1 1.6306' 1 1.6288' 1 1.6377' 1 1.6394' 1
1.6367' 1 1.6498' 1 1.6531' 1 1.6484' 1 1.6568' 1
1.6600' 1 1.6582' 1 1.6639' 1 1.6620' 1 1.6585' 1
1.6642' 1 1.6617' 1 1.6552' 1 1.6567' 1 1.6531' 1
1.6482' 1 1.6464' 1 1.6403' 1 1.6340' 1 1.6326' 1
1.6257' 1 1.6185' 1 1.6146' 1 1.6086' 1 1.6020' 1
1.5976' 1 1.5906' 1 1.5849' 1 1.5809' 1 1.5751' 1
1.5695' 1 1.5660' 1 1.5615' 1 1.5578' 1 1.5549' 1
1.5519' 1 1.5496' 1 1.5487' 1 1.5474' 1 1.5472' 1
1.5478' 1 1.5491' 1 1.5508' 1 1.5534' 1 1.5564' 1
1.5603' 1 1.5648' 1 1.5703' 1 1.5764' 1 1.5835' 1
1.5913' 1 1.6001' 1 1.6096' 1 1.6200' 1 1.6312' 1
1.6431' 1 1.6558' 1 1.6692' 1 1.6833' 1 1.6979' 1
1.7131' 1 1.7287' 1 1.7448' 1 1.7611' 1 1.7777' 1
1.7944' 1 1.8113' 1 1.8279' 1 1.8447' 1 1.8610' 1
1.8773' 1 1.8929' 1 1.9082' 1 1.9228' 1 1.9370' 1
1.9502' 1 1.9630' 1 1.9747' 1 1.9858' 1 1.9957' 1
2.0050' 1 2.0131' 1 2.0206' 1 2.0267' 1 2.0323' 1
2.0366' 1 2.0403' 1 2.0428' 1 2.0447' 1 2.0454' 1
2.0457' 1 2.0448' 1 2.0435' 1 2.0412' 1 2.0385' 1
2.0349' 1 2.0310' 1 2.0263' 1 2.0213' 1 2.0155' 1
2.0097' 1 2.0031' 1 1.9964' 1 1.9892' 1 1.9820' 1
1.9742' 1 1.9665' 1 1.9583' 1 1.9502' 1 1.9418' 1
1.9335' 1 1.9249' 1 1.9166' 1 1.9079' 1 1.8995' 1
1.8910' 1 1.8827' 1 1.8742' 1 1.8662' 1 1.8580' 1
1.8501' 1 1.8423' 1 1.8349' 1 1.8274' 1 1.8204' 1
1.8134' 1 1.8068' 1 1.8003' 1 1.7942' 1 1.7882' 1
1.7828' 1 1.7773' 1 1.7722' 1 1.7675' 1 1.7630' 1
1.7587' 1 1.7549' 1 1.7511' 1 1.7477' 1 1.7446' 1
1.7417' 1 1.7389' 1 1.7367' 1 1.7343' 1 1.7323' 1
1.7306' 1 1.7289' 1 1.7273' 1 1.7262' 1 1.7249' 1
1.7237' 1 1.7230' 1 1.7220' 1 1.7213' 1 1.7208' 1
1.7200' 1 1.7194' 1 1.7192' 1 1.7184' 1 1.7180' 1
1.7177' 1 1.7172' 1 1.7167' 1 1.7166' 1 1.7160' 1
1.7158' 1 1.7157' 1 1.7153' 1 1.7150' 1 1.7152' 1
end 47
▶EOF◀