1//<-- CLI SHELL MODE -->
2// =============================================================================
3// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
4// Copyright (C) ????-2008 - INRIA
5//
6//  This file is distributed under the same license as the Scilab package.
7// =============================================================================
8warning("off");
9
10function r=testresult(eps)
11  if argn(2)<=0 then eps=1d-10;end
12  if or(isnan(real(y))<>isnan(real(yref)))  then r=%t,return,end;
13  if or(isnan(imag(y))<>isnan(imag(yref)))  then r=%t,return,end;
14  if or(size(y)<>size(yref)) then r=%t,return,end;
15  if yref==[] then r=%f,return,end;
16  k=find(abs(yref)==0);
17  r=or(abs(y(k)-yref(k))>eps);
18  k=find(abs(yref)<>0);
19  r=r|or(max(abs(y(k)-yref(k))./abs(yref(k)))>eps);
20endfunction
21ieee(1)
22
23// ------------------------------ bessely(alpha,x)------------------------------
24//test of numerical precision reference computed by maple
25exec(SCI+'/modules/special_functions/tests/unit_tests/bessely_ref.sce',-1)
26besselY=bessely(alpha',x);
27if norm((besselY-besselY_ref)./abs(besselY_ref),1)>1d6*%eps then pause,end
28
29alpha=[];
30x=[];
31yref=[];
32y=bessely(alpha,x);
33if testresult() then pause,end
34yref=[];
35y=bessely(alpha,x,1);
36if testresult() then pause,end
37alpha=[];
38x=[0];
39yref=[];
40y=bessely(alpha,x);
41if testresult() then pause,end
42yref=[];
43y=bessely(alpha,x,1);
44if testresult() then pause,end
45alpha=[];
46x=[9.99999999999999955e-08];
47yref=[];
48y=bessely(alpha,x);
49if testresult() then pause,end
50yref=[];
51y=bessely(alpha,x,1);
52if testresult() then pause,end
53alpha=[];
54x=[1];
55yref=[];
56y=bessely(alpha,x);
57if testresult() then pause,end
58yref=[];
59y=bessely(alpha,x,1);
60if testresult() then pause,end
61alpha=[];
62x=[3 4 5];
63yref=[];
64y=bessely(alpha,x);
65if testresult() then pause,end
66yref=[];
67y=bessely(alpha,x,1);
68if testresult() then pause,end
69alpha=[];
70x=[-7.29999999999999982 8.40000000000000036];
71yref=[];
72y=bessely(alpha,x);
73if testresult() then pause,end
74yref=[];
75y=bessely(alpha,x,1);
76if testresult() then pause,end
77alpha=[];
78x=[3
794
805];
81yref=[];
82y=bessely(alpha,x);
83if testresult() then pause,end
84yref=[];
85y=bessely(alpha,x,1);
86if testresult() then pause,end
87alpha=[];
88x=[-7.29999999999999982
89 8.40000000000000036];
90yref=[];
91y=bessely(alpha,x);
92if testresult() then pause,end
93yref=[];
94y=bessely(alpha,x,1);
95if testresult() then pause,end
96alpha=[];
97x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
98yref=[];
99y=bessely(alpha,x);
100if testresult() then pause,end
101yref=[];
102y=bessely(alpha,x,1);
103if testresult() then pause,end
104alpha=[0];
105x=[];
106yref=[];
107y=bessely(alpha,x);
108if testresult() then pause,end
109yref=[];
110y=bessely(alpha,x,1);
111if testresult() then pause,end
112alpha=[0];
113x=[0];
114yref=[-%inf];
115y=bessely(alpha,x);
116if testresult() then pause,end
117yref=[-%inf];
118y=bessely(alpha,x,1);
119if testresult() then pause,end
120alpha=[0];
121x=[9.99999999999999955e-08];
122yref=[-10.3349026794206473];
123y=bessely(alpha,x);
124if testresult() then pause,end
125yref=[-10.3349026794206473];
126y=bessely(alpha,x,1);
127if testresult() then pause,end
128alpha=[0];
129x=[1];
130yref=[0.0882569642156769973];
131y=bessely(alpha,x);
132if testresult() then pause,end
133yref=[0.0882569642156769696];
134y=bessely(alpha,x,1);
135if testresult() then pause,end
136alpha=[0];
137x=[3 4 5];
138yref=[0.376850010012790504 -0.0169407393250650198 -0.308517625249033867];
139y=bessely(alpha,x);
140if testresult() then pause,end
141yref=[0.376850010012790504 -0.0169407393250650129 -0.308517625249033811];
142y=bessely(alpha,x,1);
143if testresult() then pause,end
144alpha=[0];
145x=[-7.29999999999999982 8.40000000000000036];
146yref=[0.0627738863740378017+0.576433895270028751*%i 0.266221867363905018+0*%i ];
147y=bessely(alpha,x);
148if testresult() then pause,end
149yref=[0.0627738863740374964+0.576433895270028751*%i 0.266221867363905074+0*%i ];
150y=bessely(alpha,x,1);
151if testresult() then pause,end
152alpha=[0];
153x=[3
1544
1555];
156yref=[ 0.376850010012790504
157-0.0169407393250650198
158 -0.308517625249033867];
159y=bessely(alpha,x);
160if testresult() then pause,end
161yref=[ 0.376850010012790504
162-0.0169407393250650129
163 -0.308517625249033811];
164y=bessely(alpha,x,1);
165if testresult() then pause,end
166alpha=[0];
167x=[-7.29999999999999982
168 8.40000000000000036];
169yref=[0.0627738863740378017+0.576433895270028751*%i
170 0.266221867363905018+0*%i ];
171y=bessely(alpha,x);
172if testresult() then pause,end
173yref=[0.0627738863740374964+0.576433895270028751*%i
174 0.266221867363905074+0*%i ];
175y=bessely(alpha,x,1);
176if testresult() then pause,end
177alpha=[0];
178x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
179yref=[0.0594911047067641363+0.570852419787944876*%i 0.266276386958715638+0.00106971720159772243*%i ];
180y=bessely(alpha,x);
181if testresult() then pause,end
182yref=[0.0583131019071235768+0.55954878453009449*%i 0.261003761229849252+0.00104853538181959999*%i ];
183y=bessely(alpha,x,1);
184if testresult() then pause,end
185alpha=[1];
186x=[];
187yref=[];
188y=bessely(alpha,x);
189if testresult() then pause,end
190yref=[];
191y=bessely(alpha,x,1);
192if testresult() then pause,end
193alpha=[1];
194x=[0];
195yref=[-%inf];
196y=bessely(alpha,x);
197if testresult() then pause,end
198yref=[-%inf];
199y=bessely(alpha,x,1);
200if testresult() then pause,end
201alpha=[1];
202x=[9.99999999999999955e-08];
203yref=[-6366197.72367634624];
204y=bessely(alpha,x);
205if testresult() then pause,end
206yref=[-6366197.72367634624];
207y=bessely(alpha,x,1);
208if testresult() then pause,end
209alpha=[1];
210x=[1];
211yref=[-0.781212821300288907];
212y=bessely(alpha,x);
213if testresult() then pause,end
214yref=[-0.781212821300288907];
215y=bessely(alpha,x,1);
216if testresult() then pause,end
217alpha=[1];
218x=[3 4 5];
219yref=[0.324674424791800087 0.397925710557099965 0.147863143391226859];
220y=bessely(alpha,x);
221if testresult() then pause,end
222yref=[0.324674424791800031 0.397925710557100021 0.147863143391226803];
223y=bessely(alpha,x,1);
224if testresult() then pause,end
225alpha=[1];
226x=[-7.29999999999999982 8.40000000000000036];
227yref=[0.28459437186807196-0.165140860986515869*%i -0.0534845083974161278+0*%i ];
228y=bessely(alpha,x);
229if testresult() then pause,end
230yref=[0.28459437186807196-0.165140860986515869*%i -0.0534845083974161833+0*%i ];
231y=bessely(alpha,x,1);
232if testresult() then pause,end
233alpha=[1];
234x=[3
2354
2365];
237yref=[0.324674424791800087
2380.397925710557099965
2390.147863143391226859];
240y=bessely(alpha,x);
241if testresult() then pause,end
242yref=[0.324674424791800031
2430.397925710557100021
2440.147863143391226803];
245y=bessely(alpha,x,1);
246if testresult() then pause,end
247alpha=[1];
248x=[-7.29999999999999982
249 8.40000000000000036];
250yref=[ 0.28459437186807196-0.165140860986515869*%i
251-0.0534845083974161278+0*%i ];
252y=bessely(alpha,x);
253if testresult() then pause,end
254yref=[ 0.28459437186807196-0.165140860986515869*%i
255-0.0534845083974161833+0*%i ];
256y=bessely(alpha,x,1);
257if testresult() then pause,end
258alpha=[1];
259x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
260yref=[0.273570530695402159-0.163153078769770365*%i -0.053488563419309812+0.00545213754945528262*%i ];
261y=bessely(alpha,x);
262if testresult() then pause,end
263yref=[0.268153471243458119-0.159922431356041589*%i -0.0524294189006919362+0.00534417799266206917*%i ];
264y=bessely(alpha,x,1);
265if testresult() then pause,end
266alpha=[1 2 3];
267x=[];
268yref=[];
269y=bessely(alpha,x);
270if testresult() then pause,end
271yref=[];
272y=bessely(alpha,x,1);
273if testresult() then pause,end
274alpha=[1 2 3];
275x=[0];
276yref=[-%inf -%inf -%inf];
277y=bessely(alpha,x);
278if testresult() then pause,end
279yref=[-%inf -%inf -%inf];
280y=bessely(alpha,x,1);
281if testresult() then pause,end
282alpha=[1 2 3];
283x=[9.99999999999999955e-08];
284yref=[-6366197.72367634624 -127323954473516.594 -5.09295817894065707e+21];
285y=bessely(alpha,x);
286if testresult() then pause,end
287yref=[-6366197.72367634624 -127323954473516.578 -5.09295817894065812e+21];
288y=bessely(alpha,x,1);
289if testresult() then pause,end
290alpha=[1 2 3];
291x=[1];
292yref=[-0.781212821300288907 -1.65068260681625478 -5.82151760596473089];
293y=bessely(alpha,x);
294if testresult() then pause,end
295yref=[-0.781212821300288907 -1.65068260681625478 -5.82151760596473089];
296y=bessely(alpha,x,1);
297if testresult() then pause,end
298alpha=[1 2 3];
299x=[3 4 5];
300yref=[0.324674424791800087 0.215903594603614996 0.146267162693192898];
301y=bessely(alpha,x);
302if testresult() then pause,end
303yref=[0.324674424791800031 0.215903594603615023 0.14626716269319287];
304y=bessely(alpha,x,1);
305if testresult() then pause,end
306alpha=[1 2 3];
307x=[3
3084
3095];
310yref=[0.324674424791800087 -0.16040039348492377 -0.538541616105031817
3110.397925710557099965 0.215903594603615023 -0.182022115953484942
3120.147863143391226859 0.367662882605524655 0.146267162693192898];
313y=bessely(alpha,x);
314if testresult() then pause,end
315yref=[0.324674424791800031 -0.160400393484923798 -0.538541616105031817
3160.397925710557100021 0.215903594603615051 -0.182022115953484942
3170.147863143391226803 0.367662882605524544 0.14626716269319287];
318y=bessely(alpha,x,1);
319if testresult() then pause,end
320alpha=[1 2 3];
321x=[-7.29999999999999982
322 8.40000000000000036];
323yref=[ 0.28459437186807196-0.165140860986515869*%i -0.140744947159810951-0.531189823766873737*%i -0.207473852876394749+0.456203778119049386*%i
324-0.0534845083974161278+0*%i -0.278956274125194581+0*%i -0.0793518126145812824+0*%i ];
325y=bessely(alpha,x);
326if testresult() then pause,end
327yref=[ 0.28459437186807196-0.165140860986515869*%i -0.140744947159811229-0.531189823766873737*%i -0.207473852876394638+0.456203778119049441*%i
328-0.0534845083974161833+0*%i -0.278956274125194637+0*%i -0.0793518126145812824+0*%i ];
329y=bessely(alpha,x,1);
330if testresult() then pause,end
331alpha=[1.19999999999999996 2.20000000000000018 3.20000000000000018];
332x=[];
333yref=[];
334y=bessely(alpha,x);
335if testresult() then pause,end
336yref=[];
337y=bessely(alpha,x,1);
338if testresult() then pause,end
339alpha=[1.19999999999999996 2.20000000000000018 3.20000000000000018];
340x=[0];
341yref=[-%inf -%inf -%inf];
342y=bessely(alpha,x);
343if testresult() then pause,end
344yref=[-%inf -%inf -%inf];
345y=bessely(alpha,x,1);
346if testresult() then pause,end
347alpha=[1.19999999999999996 2.20000000000000018 3.20000000000000018];
348x=[9.99999999999999955e-08];
349yref=[-168658652.388109088 -4047807657314589 -1.78103536921841748e+23];
350y=bessely(alpha,x);
351if testresult() then pause,end
352yref=[-168658652.388109088 -4047807657314589.5 -1.78103536921841714e+23];
353y=bessely(alpha,x,1);
354if testresult() then pause,end
355alpha=[1.19999999999999996 2.20000000000000018 3.20000000000000018];
356x=[1];
357yref=[-0.901214954842308846 -2.02168904229560242 -7.99421683125834193];
358y=bessely(alpha,x);
359if testresult() then pause,end
360yref=[-0.901214954842308735 -2.02168904229560242 -7.99421683125834193];
361y=bessely(alpha,x,1);
362if testresult() then pause,end
363alpha=[1.19999999999999996 2.20000000000000018 3.20000000000000018];
364x=[3 4 5];
365yref=[0.236105148830688782 0.137707238974328822 0.0774312728572489062];
366y=bessely(alpha,x);
367if testresult() then pause,end
368yref=[0.236105148830688782 0.137707238974328877 0.0774312728572488507];
369y=bessely(alpha,x,1);
370if testresult() then pause,end
371alpha=[1.19999999999999996 2.20000000000000018 3.20000000000000018];
372x=[3
3734
3745];
375yref=[0.236105148830688782 -0.249464037123650451 -0.601985736612042799
3760.404200660729673644 0.137707238974328822 -0.25272269785791196
3770.231079618553641175 0.350580558421466015 0.0774312728572489062];
378y=bessely(alpha,x);
379if testresult() then pause,end
380yref=[0.236105148830688782 -0.249464037123650617 -0.60198573661204291
381 0.4042006607296737 0.137707238974328877 -0.252722697857912071
3820.231079618553641009 0.350580558421465849 0.0774312728572488507];
383y=bessely(alpha,x,1);
384if testresult() then pause,end
385alpha=[1.19999999999999996 2.20000000000000018 3.20000000000000018];
386x=[-7.29999999999999982
387 8.40000000000000036];
388yref=[0.240226585254262215-0.173655347820190264*%i -0.0558641372965151817-0.434822373349611213*%i -0.206555050445403759+0.435740066003517668*%i
389-0.12841584805211595+0*%i -0.268892924189574611+0*%i -0.0124328265233755314+0*%i ];
390y=bessely(alpha,x);
391if testresult() then pause,end
392yref=[ 0.24022658525426227-0.173655347820190098*%i -0.0558641372965149458-0.434822373349611269*%i -0.206555050445404065+0.435740066003517612*%i
393-0.128415848052115866+0*%i -0.268892924189574556+0*%i -0.0124328265233754898+0*%i ];
394y=bessely(alpha,x,1);
395if testresult() then pause,end
396alpha=[0
3971
3982];
399x=[];
400yref=[];
401y=bessely(alpha,x);
402if testresult() then pause,end
403yref=[];
404y=bessely(alpha,x,1);
405if testresult() then pause,end
406alpha=[0
4071
4082];
409x=[0];
410yref=[-%inf
411-%inf
412-%inf];
413y=bessely(alpha,x);
414if testresult() then pause,end
415yref=[-%inf
416-%inf
417-%inf];
418y=bessely(alpha,x,1);
419if testresult() then pause,end
420alpha=[0
4211
4222];
423x=[9.99999999999999955e-08];
424yref=[-10.3349026794206473
425-6366197.72367634624
426-127323954473516.594];
427y=bessely(alpha,x);
428if testresult() then pause,end
429yref=[-10.3349026794206473
430-6366197.72367634624
431-127323954473516.578];
432y=bessely(alpha,x,1);
433if testresult() then pause,end
434alpha=[0
4351
4362];
437x=[1];
438yref=[0.0882569642156769973
439-0.781212821300288907
440 -1.65068260681625478];
441y=bessely(alpha,x);
442if testresult() then pause,end
443yref=[0.0882569642156769696
444-0.781212821300288907
445 -1.65068260681625478];
446y=bessely(alpha,x,1);
447if testresult() then pause,end
448alpha=[0
4491
4502];
451x=[3 4 5];
452yref=[ 0.376850010012790504 0.324674424791800087 -0.160400393484923798
453-0.0169407393250650198 0.397925710557099965 0.215903594603614996
454 -0.308517625249033867 0.147863143391226887 0.367662882605524655];
455y=bessely(alpha,x);
456if testresult() then pause,end
457yref=[ 0.376850010012790504 0.324674424791800031 -0.160400393484923853
458-0.0169407393250650129 0.397925710557100021 0.215903594603615023
459 -0.308517625249033811 0.147863143391226831 0.367662882605524544];
460y=bessely(alpha,x,1);
461if testresult() then pause,end
462alpha=[0
4631
4642];
465x=[-7.29999999999999982 8.40000000000000036];
466yref=[0.0627738863740378017+0.576433895270028751*%i 0.28459437186807196-0.165140860986515869*%i -0.140744947159810951-0.531189823766873737*%i
467 0.266221867363905018+0*%i -0.0534845083974161417+0*%i -0.278956274125194581+0*%i ];
468y=bessely(alpha,x);
469if testresult() then pause,end
470yref=[0.0627738863740374964+0.576433895270028751*%i 0.28459437186807196-0.165140860986515869*%i -0.140744947159811201-0.531189823766873737*%i
471 0.266221867363905074+0*%i -0.0534845083974161833+0*%i -0.278956274125194637+0*%i ];
472y=bessely(alpha,x,1);
473if testresult() then pause,end
474alpha=[0
4751
4762];
477x=[3
4784
4795];
480yref=[0.376850010012790504
4810.397925710557099965
4820.367662882605524655];
483y=bessely(alpha,x);
484if testresult() then pause,end
485yref=[0.376850010012790504
4860.397925710557100021
4870.367662882605524544];
488y=bessely(alpha,x,1);
489if testresult() then pause,end
490alpha=[0
4911
4922];
493x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
494yref=[0.0594911047067641363+0.570852419787944876*%i 0.273570530695402159-0.163153078769770365*%i -0.134563835842271912-0.526358624867831137*%i
495 0.266276386958715638+0.00106971720159772243*%i -0.0534885634193098258+0.00545213754945528262*%i -0.279008596242427576+0.000258725570376667352*%i ];
496y=bessely(alpha,x);
497if testresult() then pause,end
498yref=[0.0583131019071235768+0.55954878453009449*%i 0.268153471243458119-0.159922431356041589*%i -0.131899293367662818-0.515936025779016094*%i
499 0.261003761229849252+0.00104853538181959999*%i -0.0524294189006919362+0.00534417799266206917*%i -0.273483855878007609+0.000253602460833710075*%i ];
500y=bessely(alpha,x,1);
501if testresult() then pause,end
502alpha=[0.0200000000000000004
503 1.02000000000000002
504 2.02000000000000002];
505x=[];
506yref=[];
507y=bessely(alpha,x);
508if testresult() then pause,end
509yref=[];
510y=bessely(alpha,x,1);
511if testresult() then pause,end
512alpha=[0.0200000000000000004
513 1.02000000000000002
514 2.02000000000000002];
515x=[0];
516yref=[-%inf
517-%inf
518-%inf];
519y=bessely(alpha,x);
520if testresult() then pause,end
521yref=[-%inf
522-%inf
523-%inf];
524y=bessely(alpha,x,1);
525if testresult() then pause,end
526alpha=[0.0200000000000000004
527 1.02000000000000002
528 2.02000000000000002];
529x=[9.99999999999999955e-08];
530yref=[-10.5435018262277929
531-8811068.85048608482
532-179745804549905.625];
533y=bessely(alpha,x);
534if testresult() then pause,end
535yref=[-10.5435018262277929
536-8811068.85048608482
537-179745804549905.656];
538y=bessely(alpha,x,1);
539if testresult() then pause,end
540alpha=[0.0200000000000000004
541 1.02000000000000002
542 2.02000000000000002];
543x=[1];
544yref=[0.0643040307807657419
545-0.793245652497184595
546 -1.68252516187502232];
547y=bessely(alpha,x);
548if testresult() then pause,end
549yref=[0.0643040307807657419
550-0.793245652497184706
551 -1.68252516187502277];
552y=bessely(alpha,x,1);
553if testresult() then pause,end
554alpha=[0.0200000000000000004
555 1.02000000000000002
556 2.02000000000000002];
557x=[3 4 5];
558yref=[ 0.384820425921874987 0.316452158174697407 -0.169632958363080733
559-0.00447693217669133808 0.399776541248366868 0.208362968213358424
560 -0.302795423995985136 0.15685459545204769 0.366792098940420641];
561y=bessely(alpha,x);
562if testresult() then pause,end
563yref=[ 0.384820425921874987 0.316452158174697407 -0.169632958363080816
564-0.00447693217669134502 0.399776541248366868 0.208362968213358452
565 -0.302795423995985136 0.156854595452047746 0.366792098940420641];
566y=bessely(alpha,x,1);
567if testresult() then pause,end
568alpha=[0.0200000000000000004
569 1.02000000000000002
570 2.02000000000000002];
571x=[-7.29999999999999982 8.40000000000000036];
572yref=[0.0535918058658350627+0.575574764795366067*%i 0.286331738254591017-0.166506620507523545*%i -0.133607798473967337-0.529044147557647282*%i
573 0.26392041753001122+0*%i -0.0613287659781691918+0*%i -0.278814546410423747+0*%i ];
574y=bessely(alpha,x);
575if testresult() then pause,end
576yref=[0.0535918058658351945+0.575574764795366067*%i 0.286331738254590962-0.166506620507523573*%i -0.13360779847396731-0.52904414755764706*%i
577 0.26392041753001122+0*%i -0.061328765978169178+0*%i -0.278814546410423747+0*%i ];
578y=bessely(alpha,x,1);
579if testresult() then pause,end
580alpha=[0.0200000000000000004
581 1.02000000000000002
582 2.02000000000000002];
583x=[3
5844
5855];
586yref=[0.384820425921874987
5870.399776541248366868
5880.366792098940420641];
589y=bessely(alpha,x);
590if testresult() then pause,end
591yref=[0.384820425921874987
5920.399776541248366868
5930.366792098940420641];
594y=bessely(alpha,x,1);
595if testresult() then pause,end
596alpha=[0.0200000000000000004
597 1.02000000000000002
598 2.02000000000000002];
599x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
600yref=[0.0503114628539832326+0.569955434445619136*%i 0.27533846887479968-0.164682242384238758*%i -0.127380869613414494-0.524145792482789141*%i
601 0.263974678263304063+0.00123918133456109997*%i -0.061334389730632502+0.00542770518273297253*%i -0.278866950019879489+0.000114433428284946292*%i ];
602y=bessely(alpha,x);
603if testresult() then pause,end
604yref=[0.0493152291415966115+0.55866956068757101*%i 0.269886401901392015-0.161421315502212115*%i -0.12485855939972948-0.513767010410947855*%i
605 0.258747629420268188+0.00121464390012337303*%i -0.0601198874420454457+0.00532022941921519998*%i -0.273345014438587053+0.000112167494586902505*%i ];
606y=bessely(alpha,x,1);
607if testresult() then pause,end
608alpha=[-1];
609x=[];
610yref=[];
611y=bessely(alpha,x);
612if testresult() then pause,end
613yref=[];
614y=bessely(alpha,x,1);
615if testresult() then pause,end
616alpha=[-1];
617x=[0];
618yref=[-%inf];
619y=bessely(alpha,x);
620if testresult() then pause,end
621yref=[-%inf];
622y=bessely(alpha,x,1);
623if testresult() then pause,end
624alpha=[-1];
625x=[9.99999999999999955e-08];
626yref=[6366197.72367634624];
627y=bessely(alpha,x);
628if testresult() then pause,end
629yref=[6366197.72367634624];
630y=bessely(alpha,x,1);
631if testresult() then pause,end
632alpha=[-1];
633x=[1];
634yref=[0.781212821300288907];
635y=bessely(alpha,x);
636if testresult() then pause,end
637yref=[0.781212821300288907];
638y=bessely(alpha,x,1);
639if testresult() then pause,end
640alpha=[-1];
641x=[3 4 5];
642yref=[-0.324674424791800087 -0.397925710557099965 -0.147863143391226859];
643y=bessely(alpha,x);
644if testresult() then pause,end
645yref=[-0.324674424791800031 -0.397925710557100021 -0.147863143391226803];
646y=bessely(alpha,x,1);
647if testresult() then pause,end
648alpha=[-1];
649x=[-7.29999999999999982 8.40000000000000036];
650yref=[-0.28459437186807196+0.165140860986515869*%i 0.0534845083974161278+0*%i ];
651y=bessely(alpha,x);
652if testresult() then pause,end
653yref=[-0.28459437186807196+0.165140860986515869*%i 0.0534845083974161833+0*%i ];
654y=bessely(alpha,x,1);
655if testresult() then pause,end
656alpha=[-1];
657x=[3
6584
6595];
660yref=[-0.324674424791800087
661-0.397925710557099965
662-0.147863143391226859];
663y=bessely(alpha,x);
664if testresult() then pause,end
665yref=[-0.324674424791800031
666-0.397925710557100021
667-0.147863143391226803];
668y=bessely(alpha,x,1);
669if testresult() then pause,end
670alpha=[-1];
671x=[-7.29999999999999982
672 8.40000000000000036];
673yref=[ -0.28459437186807196+0.165140860986515869*%i
6740.0534845083974161278+0*%i ];
675y=bessely(alpha,x);
676if testresult() then pause,end
677yref=[ -0.28459437186807196+0.165140860986515869*%i
6780.0534845083974161833+0*%i ];
679y=bessely(alpha,x,1);
680if testresult() then pause,end
681alpha=[-1];
682x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
683yref=[-0.273570530695402159+0.163153078769770365*%i 0.053488563419309812-0.00545213754945528262*%i ];
684y=bessely(alpha,x);
685if testresult() then pause,end
686yref=[-0.268153471243458119+0.159922431356041589*%i 0.0524294189006919362-0.00534417799266206917*%i ];
687y=bessely(alpha,x,1);
688if testresult() then pause,end
689alpha=[ -1
690-1.5];
691x=[];
692yref=[];
693y=bessely(alpha,x);
694if testresult() then pause,end
695yref=[];
696y=bessely(alpha,x,1);
697if testresult() then pause,end
698alpha=[ -1
699-1.5];
700x=[0];
701yref=[-%inf
702-%inf];
703y=bessely(alpha,x);
704if testresult() then pause,end
705yref=[-%inf
706-%inf];
707y=bessely(alpha,x,1);
708if testresult() then pause,end
709alpha=[ -1
710-1.5];
711x=[9.99999999999999955e-08];
712yref=[ 6366197.72367634624
713-8.4104417400672E-12];
714y=bessely(alpha,x);
715if testresult() then pause,end
716yref=[ 6366197.72367634624
717 -8.4104417400672E-12];
718y=bessely(alpha,x,1);
719if testresult() then pause,end
720alpha=[ -1
721-1.5];
722x=[1];
723yref=[ 0.781212821300288907
724-0.240297839123427059];
725y=bessely(alpha,x);
726if testresult() then pause,end
727yref=[ 0.781212821300288907
728-0.240297839123427059];
729y=bessely(alpha,x,1);
730if testresult() then pause,end
731alpha=[ -1
732-1.5];
733x=[3 4 5];
734yref=[-0.324674424791800087 -0.477718215087092224
735-0.397925710557099965 -0.185285948354267732
736-0.147863143391226859 0.169651306144741165];
737y=bessely(alpha,x);
738if testresult() then pause,end
739yref=[-0.324674424791800031 -0.477718215087092224
740-0.397925710557100021 -0.185285948354267732
741-0.147863143391226803 0.169651306144741165];
742y=bessely(alpha,x,1);
743if testresult() then pause,end
744alpha=[ -1
745-1.5];
746x=[-7.29999999999999982 8.40000000000000036];
747yref=[ -0.28459437186807196+0.165140860986515869*%i -1.66533453693773654e-16-0.120953010973629607*%i
7480.0534845083974161278+0*%i -0.170966347821113557+0*%i ];
749y=bessely(alpha,x);
750if testresult() then pause,end
751yref=[ -0.28459437186807196+0.165140860986515869*%i -1.66533453693773654e-16-0.120953010973629607*%i
7520.0534845083974161833+0*%i -0.170966347821113557+0*%i ];
753y=bessely(alpha,x,1);
754if testresult() then pause,end
755alpha=[ -1
756-1.5];
757x=[-7.29999999999999982
758 8.40000000000000036];
759yref=[ -0.28459437186807196+0.165140860986515869*%i
760-0.170966347821113557+0*%i ];
761y=bessely(alpha,x);
762if testresult() then pause,end
763yref=[ -0.28459437186807196+0.165140860986515869*%i
764-0.170966347821113557+0*%i ];
765y=bessely(alpha,x,1);
766if testresult() then pause,end
767alpha=[ -1
768-1.5];
769x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
770yref=[-0.273570530695402159+0.163153078769770365*%i 0.00552027940993567862-0.120968618958143467*%i
771 0.053488563419309812-0.00545213754945528262*%i -0.171004326769525289-0.00409499672885600772*%i ];
772y=bessely(alpha,x);
773if testresult() then pause,end
774yref=[-0.268153471243458119+0.159922431356041589*%i 0.00541097055390155068-0.118573279814522647*%i
7750.0524294189006919362-0.00534417799266206917*%i -0.16761821422920356-0.00401391036082020345*%i ];
776y=bessely(alpha,x,1);
777if testresult() then pause,end
778alpha=[-3
779-2
780-1
781 0
782 1];
783x=[];
784yref=[];
785y=bessely(alpha,x);
786if testresult() then pause,end
787yref=[];
788y=bessely(alpha,x,1);
789if testresult() then pause,end
790alpha=[-3
791-2
792-1
793 0
794 1];
795x=[0];
796yref=[-%inf
797-%inf
798-%inf
799-%inf
800-%inf];
801y=bessely(alpha,x);
802if testresult() then pause,end
803yref=[-%inf
804-%inf
805-%inf
806-%inf
807-%inf];
808y=bessely(alpha,x,1);
809if testresult() then pause,end
810alpha=[-3
811-2
812-1
813 0
814 1];
815x=[9.99999999999999955e-08];
816yref=[5.09295817894065707e+21
817 -127323954473516.594
818 6366197.72367634624
819 -10.3349026794206473
820 -6366197.72367634624];
821y=bessely(alpha,x);
822if testresult() then pause,end
823yref=[5.09295817894065812e+21
824 -127323954473516.578
825 6366197.72367634624
826 -10.3349026794206473
827 -6366197.72367634624];
828y=bessely(alpha,x,1);
829if testresult() then pause,end
830alpha=[-3
831-2
832-1
833 0
834 1];
835x=[1];
836yref=[ 5.82151760596473089
837 -1.65068260681625478
838 0.781212821300288907
8390.0882569642156769973
840-0.781212821300288907];
841y=bessely(alpha,x);
842if testresult() then pause,end
843yref=[ 5.82151760596473089
844 -1.65068260681625478
845 0.781212821300288907
8460.0882569642156769696
847-0.781212821300288907];
848y=bessely(alpha,x,1);
849if testresult() then pause,end
850alpha=[-3
851-2
852-1
853 0
854 1];
855x=[3 4 5];
856yref=[ 0.538541616105031817 -0.160400393484923798 -0.324674424791800087 0.376850010012790504 0.324674424791800087
857 0.182022115953484942 0.215903594603614996 -0.397925710557099965 -0.0169407393250650198 0.397925710557099965
858-0.146267162693192898 0.367662882605524655 -0.147863143391226859 -0.308517625249033867 0.147863143391226859];
859y=bessely(alpha,x);
860if testresult() then pause,end
861yref=[0.538541616105031817 -0.160400393484923853 -0.324674424791800031 0.376850010012790504 0.324674424791800031
8620.182022115953484942 0.215903594603615023 -0.397925710557100021 -0.0169407393250650129 0.397925710557100021
863-0.14626716269319287 0.367662882605524544 -0.147863143391226803 -0.308517625249033811 0.147863143391226803];
864y=bessely(alpha,x,1);
865if testresult() then pause,end
866alpha=[-3
867-2
868-1
869 0
870 1];
871x=[-7.29999999999999982 8.40000000000000036];
872yref=[ 0.207473852876394749-0.456203778119049386*%i -0.140744947159810951-0.531189823766873737*%i -0.28459437186807196+0.165140860986515869*%i 0.0627738863740378017+0.576433895270028751*%i 0.28459437186807196-0.165140860986515869*%i
8730.0793518126145812824+0*%i -0.278956274125194581+0*%i 0.0534845083974161278+0*%i 0.266221867363905018+0*%i -0.0534845083974161278+0*%i ];
874y=bessely(alpha,x);
875if testresult() then pause,end
876yref=[ 0.207473852876394638-0.456203778119049441*%i -0.140744947159811201-0.531189823766873737*%i -0.28459437186807196+0.165140860986515869*%i 0.0627738863740374964+0.576433895270028751*%i 0.28459437186807196-0.165140860986515869*%i
8770.0793518126145812824+0*%i -0.278956274125194637+0*%i 0.0534845083974161833+0*%i 0.266221867363905074+0*%i -0.0534845083974161833+0*%i ];
878y=bessely(alpha,x,1);
879if testresult() then pause,end
880alpha=[-3
881-2
882-1
883 0
884 1];
885x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
886yref=[ 0.200627649489689364-0.451768607139008516*%i -0.134563835842271912-0.526358624867831137*%i -0.273570530695402159+0.163153078769770365*%i 0.0594911047067641363+0.570852419787944876*%i 0.273570530695402159-0.163153078769770365*%i
8870.0793716263712800352+0.00501260111167922684*%i -0.279008596242427576+0.000258725570376667352*%i 0.053488563419309812-0.00545213754945528262*%i 0.266276386958715638+0.00106971720159772243*%i -0.053488563419309812+0.00545213754945528262*%i ];
888y=bessely(alpha,x);
889if testresult() then pause,end
890yref=[ 0.196654955858446223-0.44282298935929687*%i -0.131899293367662818-0.515936025779016094*%i -0.268153471243458119+0.159922431356041589*%i 0.0583131019071235768+0.55954878453009449*%i 0.268153471243458119-0.159922431356041589*%i
8910.0777999628673285148+0.00491334495948420047*%i -0.273483855878007609+0.000253602460833710075*%i 0.0524294189006919362-0.00534417799266206917*%i 0.261003761229849252+0.00104853538181959999*%i -0.0524294189006919362+0.00534417799266206917*%i ];
892y=bessely(alpha,x,1);
893if testresult() then pause,end
894alpha=[-2
895-1
896 0
897 1];
898x=[];
899yref=[];
900y=bessely(alpha,x);
901if testresult() then pause,end
902yref=[];
903y=bessely(alpha,x,1);
904if testresult() then pause,end
905alpha=[-2
906-1
907 0
908 1];
909x=[0];
910yref=[-%inf
911-%inf
912-%inf
913-%inf];
914y=bessely(alpha,x);
915if testresult() then pause,end
916yref=[-%inf
917-%inf
918-%inf
919-%inf];
920y=bessely(alpha,x,1);
921if testresult() then pause,end
922alpha=[-2
923-1
924 0
925 1];
926x=[9.99999999999999955e-08];
927yref=[-127323954473516.594
928 6366197.72367634624
929-10.3349026794206473
930-6366197.72367634624];
931y=bessely(alpha,x);
932if testresult() then pause,end
933yref=[-127323954473516.578
934 6366197.72367634624
935-10.3349026794206473
936-6366197.72367634624];
937y=bessely(alpha,x,1);
938if testresult() then pause,end
939alpha=[-2
940-1
941 0
942 1];
943x=[1];
944yref=[ -1.65068260681625478
945 0.781212821300288907
9460.0882569642156769973
947-0.781212821300288907];
948y=bessely(alpha,x);
949if testresult() then pause,end
950yref=[ -1.65068260681625478
951 0.781212821300288907
9520.0882569642156769696
953-0.781212821300288907];
954y=bessely(alpha,x,1);
955if testresult() then pause,end
956alpha=[-2
957-1
958 0
959 1];
960x=[3 4 5];
961yref=[-0.160400393484923798 -0.324674424791800087 0.376850010012790504 0.324674424791800087
962 0.215903594603614996 -0.397925710557099965 -0.0169407393250650198 0.397925710557099965
963 0.367662882605524655 -0.147863143391226859 -0.308517625249033867 0.147863143391226859];
964y=bessely(alpha,x);
965if testresult() then pause,end
966yref=[-0.160400393484923853 -0.324674424791800031 0.376850010012790504 0.324674424791800031
967 0.215903594603615023 -0.397925710557100021 -0.0169407393250650129 0.397925710557100021
968 0.367662882605524544 -0.147863143391226803 -0.308517625249033811 0.147863143391226803];
969y=bessely(alpha,x,1);
970if testresult() then pause,end
971alpha=[-2
972-1
973 0
974 1];
975x=[-7.29999999999999982 8.40000000000000036];
976yref=[-0.140744947159810951-0.531189823766873737*%i -0.28459437186807196+0.165140860986515869*%i 0.0627738863740378017+0.576433895270028751*%i 0.28459437186807196-0.165140860986515869*%i
977-0.278956274125194581+0*%i 0.0534845083974161278+0*%i 0.266221867363905018+0*%i -0.0534845083974161278+0*%i ];
978y=bessely(alpha,x);
979if testresult() then pause,end
980yref=[-0.140744947159811201-0.531189823766873737*%i -0.28459437186807196+0.165140860986515869*%i 0.0627738863740374964+0.576433895270028751*%i 0.28459437186807196-0.165140860986515869*%i
981-0.278956274125194637+0*%i 0.0534845083974161833+0*%i 0.266221867363905074+0*%i -0.0534845083974161833+0*%i ];
982y=bessely(alpha,x,1);
983if testresult() then pause,end
984alpha=[-2
985-1
986 0
987 1];
988x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
989yref=[-0.134563835842271912-0.526358624867831137*%i -0.273570530695402159+0.163153078769770365*%i 0.0594911047067641363+0.570852419787944876*%i 0.273570530695402159-0.163153078769770365*%i
990-0.279008596242427576+0.000258725570376667352*%i 0.053488563419309812-0.00545213754945528262*%i 0.266276386958715638+0.00106971720159772243*%i -0.053488563419309812+0.00545213754945528262*%i ];
991y=bessely(alpha,x);
992if testresult() then pause,end
993yref=[-0.131899293367662818-0.515936025779016094*%i -0.268153471243458119+0.159922431356041589*%i 0.0583131019071235768+0.55954878453009449*%i 0.268153471243458119-0.159922431356041589*%i
994-0.273483855878007609+0.000253602460833710075*%i 0.0524294189006919362-0.00534417799266206917*%i 0.261003761229849252+0.00104853538181959999*%i -0.0524294189006919362+0.00534417799266206917*%i ];
995y=bessely(alpha,x,1);
996if testresult() then pause,end
997alpha=[ -1.30000000000000004
998-0.300000000000000044
999 0.699999999999999956
1000 1.69999999999999996];
1001x=[];
1002yref=[];
1003y=bessely(alpha,x);
1004if testresult() then pause,end
1005yref=[];
1006y=bessely(alpha,x,1);
1007if testresult() then pause,end
1008alpha=[ -1.30000000000000004
1009-0.300000000000000044
1010 0.699999999999999956
1011 1.69999999999999996];
1012x=[0];
1013yref=[-%inf
1014-%inf
1015-%inf
1016-%inf];
1017y=bessely(alpha,x);
1018if testresult() then pause,end
1019yref=[-%inf
1020-%inf
1021-%inf
1022-%inf];
1023y=bessely(alpha,x,1);
1024if testresult() then pause,end
1025alpha=[ -1.30000000000000004
1026-0.300000000000000044
1027 0.699999999999999956
1028 1.69999999999999996];
1029x=[9.99999999999999955e-08];
1030yref=[ 520508809.337161601
1031-86.7425820843251785
1032-53316.8312102425552
1033-746435636856.653198];
1034y=bessely(alpha,x);
1035if testresult() then pause,end
1036yref=[ 520508809.337161541
1037-86.7425820843251785
1038-53316.8312102425552
1039-746435636856.653198];
1040y=bessely(alpha,x,1);
1041if testresult() then pause,end
1042alpha=[ -1.30000000000000004
1043-0.300000000000000044
1044 0.699999999999999956
1045 1.69999999999999996];
1046x=[1];
1047yref=[ 0.314183689901754548
1048 0.454431262899763211
1049-0.586842447641611864
1050 -1.27601068959802011];
1051y=bessely(alpha,x);
1052if testresult() then pause,end
1053yref=[ 0.314183689901754548
1054 0.454431262899763155
1055-0.586842447641611864
1056 -1.27601068959802033];
1057y=bessely(alpha,x,1);
1058if testresult() then pause,end
1059alpha=[ -1.30000000000000004
1060-0.300000000000000044
1061 0.699999999999999956
1062 1.69999999999999996];
1063x=[3 4 5];
1064yref=[-0.466080611113182774 0.212390924829804428 0.423602426147222433 -0.014709792627767504
1065-0.306380088504755232 -0.199425513451468028 0.336293915522475506 0.317128383884334675
10660.0416545814756713872 -0.355966921924181279 0.00106144915522857208 0.356264127687647869];
1067y=bessely(alpha,x);
1068if testresult() then pause,end
1069yref=[-0.466080611113182774 0.212390924829804484 0.423602426147222433 -0.0147097926277675942
1070-0.306380088504755177 -0.199425513451468028 0.336293915522475562 0.317128383884334619
10710.0416545814756714428 -0.355966921924181279 0.00106144915522855821 0.356264127687647758];
1072y=bessely(alpha,x,1);
1073if testresult() then pause,end
1074alpha=[ -1.30000000000000004
1075-0.300000000000000044
1076 0.699999999999999956
1077 1.69999999999999996];
1078x=[-7.29999999999999982 8.40000000000000036];
1079yref=[ -0.121666717696636423+0.0833864463893779634*%i 0.110649033008904024+0.419447053092179445*%i 0.130761158765862334-0.0489113461352237433*%i -0.135726515511943424-0.4100667949292639*%i
1080-0.0862806191484301321+0*%i 0.268366926287593732+0*%i 0.0671115529850285913+0*%i -0.257181667456753915+0*%i ];
1081y=bessely(alpha,x);
1082if testresult() then pause,end
1083yref=[ -0.12166671769663634+0.0833864463893780328*%i 0.110649033008904024+0.4194470530921795*%i 0.130761158765862362-0.0489113461352239931*%i -0.135726515511943285-0.4100667949292639*%i
1084-0.086280619148430146+0*%i 0.268366926287593677+0*%i 0.0671115529850285913+0*%i -0.257181667456753915+0*%i ];
1085y=bessely(alpha,x,1);
1086if testresult() then pause,end
1087alpha=[ -1.30000000000000004
1088-0.300000000000000044
1089 0.699999999999999956
1090 1.69999999999999996];
1091x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
1092yref=[ -0.113594218202467456+0.080767182276823099*%i 0.109351445375408682+0.417004570698233712*%i 0.122488039151380818-0.046468467367517037*%i -0.132866541475518879-0.408157234754056908*%i
1093-0.0863035346120157726-0.00510058899774245254*%i 0.268418706707374222-0.00153406339209650328*%i 0.0671311394293841734+0.00525581351428117127*%i -0.25722816126732434+0.0023833881553354877*%i ];
1094y=bessely(alpha,x);
1095if testresult() then pause,end
1096yref=[ -0.111344901977376662+0.0791678849144669128*%i 0.107186141681151675+0.408747326961261703*%i 0.120062613472129431-0.0455483300642385985*%i -0.130235607681160659-0.400075180006480513*%i
1097-0.0845946101283816532-0.00499959056867008752*%i 0.263103660205283252-0.00150368690170144784*%i 0.0658018538062531688+0.00515174143384616867*%i -0.252134702411367373+0.00233619390783472436*%i ];
1098y=bessely(alpha,x,1);
1099if testresult() then pause,end
1100alpha=[ -2.29999999999999982
1101 -1.30000000000000004
1102-0.300000000000000044
1103 0.699999999999999956
1104 1.69999999999999996];
1105x=[];
1106yref=[];
1107y=bessely(alpha,x);
1108if testresult() then pause,end
1109yref=[];
1110y=bessely(alpha,x,1);
1111if testresult() then pause,end
1112alpha=[ -2.29999999999999982
1113 -1.30000000000000004
1114-0.300000000000000044
1115 0.699999999999999956
1116 1.69999999999999996];
1117x=[0];
1118yref=[-%inf
1119-%inf
1120-%inf
1121-%inf
1122-%inf];
1123y=bessely(alpha,x);
1124if testresult() then pause,end
1125yref=[-%inf
1126-%inf
1127-%inf
1128-%inf
1129-%inf];
1130y=bessely(alpha,x,1);
1131if testresult() then pause,end
1132alpha=[ -2.29999999999999982
1133 -1.30000000000000004
1134-0.300000000000000044
1135 0.699999999999999956
1136 1.69999999999999996];
1137x=[9.99999999999999955e-08];
1138yref=[ -13533229042766080
1139 520508809.337161601
1140-86.7425820843251785
1141-53316.8312102425552
1142-746435636856.653198];
1143y=bessely(alpha,x);
1144if testresult() then pause,end
1145yref=[ -13533229042766082
1146 520508809.337161541
1147-86.7425820843251785
1148-53316.8312102425552
1149-746435636856.653198];
1150y=bessely(alpha,x,1);
1151if testresult() then pause,end
1152alpha=[ -2.29999999999999982
1153 -1.30000000000000004
1154-0.300000000000000044
1155 0.699999999999999956
1156 1.69999999999999996];
1157x=[1];
1158yref=[ -1.27130885664432558
1159 0.314183689901754548
1160 0.454431262899763211
1161-0.586842447641611864
1162 -1.27601068959802011];
1163y=bessely(alpha,x);
1164if testresult() then pause,end
1165yref=[ -1.2713088566443258
1166 0.314183689901754548
1167 0.454431262899763155
1168-0.586842447641611864
1169 -1.27601068959802033];
1170y=bessely(alpha,x,1);
1171if testresult() then pause,end
1172alpha=[ -2.29999999999999982
1173 -1.30000000000000004
1174-0.300000000000000044
1175 0.699999999999999956
1176 1.69999999999999996];
1177x=[3 4 5];
1178yref=[0.191545604801620084 -0.466080611113182774 0.212390924829804428 0.423602426147222433 -0.014709792627767504
11790.398572570979558172 -0.306380088504755232 -0.199425513451468028 0.336293915522475506 0.317128383884334675
11800.334306539556831872 0.0416545814756713872 -0.355966921924181279 0.00106144915522857208 0.356264127687647869];
1181y=bessely(alpha,x);
1182if testresult() then pause,end
1183yref=[0.191545604801620056 -0.466080611113182774 0.212390924829804484 0.423602426147222433 -0.0147097926277675942
11840.398572570979558116 -0.306380088504755177 -0.199425513451468028 0.336293915522475562 0.317128383884334619
11850.334306539556831872 0.0416545814756714428 -0.355966921924181279 0.00106144915522855821 0.356264127687647758];
1186y=bessely(alpha,x,1);
1187if testresult() then pause,end
1188alpha=[ -2.29999999999999982
1189 -1.30000000000000004
1190-0.300000000000000044
1191 0.699999999999999956
1192 1.69999999999999996];
1193x=[-7.29999999999999982 8.40000000000000036];
1194yref=[-0.153982384517294985-0.389747770816509909*%i -0.121666717696636423+0.0833864463893779634*%i 0.110649033008904024+0.419447053092179445*%i 0.130761158765862334-0.0489113461352237433*%i -0.135726515511943424-0.4100667949292639*%i
1195-0.241661020360698642+0*%i -0.0862806191484301321+0*%i 0.268366926287593732+0*%i 0.0671115529850285913+0*%i -0.257181667456753915+0*%i ];
1196y=bessely(alpha,x);
1197if testresult() then pause,end
1198yref=[-0.153982384517294929-0.38974777081651002*%i -0.12166671769663634+0.0833864463893780328*%i 0.110649033008904024+0.4194470530921795*%i 0.130761158765862362-0.0489113461352239931*%i -0.135726515511943285-0.4100667949292639*%i
1199-0.241661020360698614+0*%i -0.086280619148430146+0*%i 0.268366926287593677+0*%i 0.0671115529850285913+0*%i -0.257181667456753915+0*%i ];
1200y=bessely(alpha,x,1);
1201if testresult() then pause,end
1202alpha=[ -2.29999999999999982
1203 -1.30000000000000004
1204-0.300000000000000044
1205 0.699999999999999956
1206 1.69999999999999996];
1207x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
1208yref=[-0.149888167842517933-0.388349236524206054*%i -0.113594218202467456+0.080767182276823099*%i 0.109351445375408682+0.417004570698233712*%i 0.122488039151380818-0.046468467367517037*%i -0.132866541475518879-0.408157234754056908*%i
1209-0.241702100416107291+0.00304920616213296336*%i -0.0863035346120157726-0.00510058899774245254*%i 0.268418706707374222-0.00153406339209650328*%i 0.0671311394293841734+0.00525581351428117127*%i -0.25722816126732434+0.0023833881553354877*%i ];
1210y=bessely(alpha,x);
1211if testresult() then pause,end
1212yref=[-0.146920183263616416-0.380659406420718049*%i -0.111344901977376662+0.0791678849144669128*%i 0.107186141681151675+0.408747326961261703*%i 0.120062613472129431-0.0455483300642385985*%i -0.130235607681160659-0.400075180006480513*%i
1213-0.236916078163324606+0.00298882783476134631*%i -0.0845946101283816532-0.00499959056867008752*%i 0.263103660205283252-0.00150368690170144784*%i 0.0658018538062531688+0.00515174143384616867*%i -0.252134702411367373+0.00233619390783472436*%i ];
1214y=bessely(alpha,x,1);
1215if testresult() then pause,end
1216