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