Reduce test precision as the test fail on Linux 32bits
[scilab.git] / scilab / modules / special_functions / tests / unit_tests / erf.tst
1 //<-- CLI SHELL MODE -->
2 // =============================================================================
3 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
4 // Copyright (C) ????-2008 - INRIA Serge Steer
5 //
6 //  This file is distributed under the same license as the Scilab package.
7 // =============================================================================
8
9
10
11 function [x,y]=erf_ref()
12     //reference Yr=erf(X) generee par Maple
13     //pm:=array(1..1001):x:=-5:
14     //for k from 1  to 1001 do
15     //pm[k]:=evalf[24](erf(x)):
16     //x:=x+0.01:end do;
17
18
19     x=-5:0.01:5;
20
21     y=[-0.999999999998462540205572
22     -0.999999999998297733894900
23     -0.999999999998115631480261
24     -0.999999999997914458007001
25     -0.999999999997692260661541
26     -0.999999999997446891397149
27     -0.999999999997175987908301
28     -0.999999999996876952801260
29     -0.999999999996546930794861
30     -0.999999999996182783770717
31     -0.999999999995781063475994
32     -0.999999999995337981664472
33     -0.999999999994849377442675
34     -0.999999999994310681567322
35     -0.999999999993716877418075
36     -0.999999999993062458345374
37     -0.999999999992341381066975
38     -0.999999999991547014758366
39     -0.999999999990672085451468
40     -0.999999999989708615322670
41     -0.999999999988647856415078
42     -0.999999999987480218300757
43     -0.999999999986195189146334
44     -0.999999999984781249599477
45     -0.999999999983225778864110
46     -0.999999999981514952278515
47     -0.999999999979633629652374
48     -0.999999999977565233555991
49     -0.999999999975291616686970
50     -0.999999999972792917366275
51     -0.999999999970047402136203
52     -0.999999999967031294347156
53     -0.999999999963718587527495
54     -0.999999999960080842230859
55     -0.999999999956086964947412
56     -0.999999999951702967549089
57     -0.999999999946891705613252
58     -0.999999999941612593833703
59     -0.999999999935821296581860
60     -0.999999999929469391523360
61     -0.999999999922504004025582
62     -0.999999999914867409908580
63     -0.999999999906496603894840
64     -0.999999999897322830900966
65     -0.999999999887271077085832
66     -0.999999999876259517323697
67     -0.999999999864198915505960
68     -0.999999999850991973790366
69     -0.999999999836532626610029
70     -0.999999999820705274925103
71     -0.999999999803383955845711
72     -0.999999999784431442373954
73     -0.999999999763698267603729
74     -0.999999999741021667277623
75     -0.999999999716224434128131
76     -0.999999999689113676923764
77     -0.999999999659479476596644
78     -0.999999999627093431244556
79     -0.999999999591707081174197
80     -0.999999999553050204480830
81     -0.999999999510828972939411
82     -0.999999999464723957210342
83     -0.999999999414387969535759
84     -0.999999999359443731215952
85     -0.999999999299481351206210
86     -0.999999999234055601157792
87     -0.999999999162682971138542
88     -0.999999999084838489103913
89     -0.999999998999952285943071
90     -0.999999998907405886591731
91     -0.999999998806528206277959
92     -0.999999998696591229443244
93     -0.999999998576805347252390
94     -0.999999998446314327865434
95     -0.999999998304189891785759
96     -0.999999998149425862613257
97     -0.999999997980931861411821
98     -0.999999997797526510638154
99     -0.999999997597930111164948
100     -0.999999997380756753356393
101     -0.999999997144505820407811
102     -0.999999996887552839233372
103     -0.999999996608139631065173
104     -0.999999996304363710601733
105     -0.999999995974166879001692
106     -0.999999995615322952246180
107     -0.999999995225424562376994
108     -0.999999994801868964843021
109     -0.999999994341842780638679
110     -0.999999993842305597079599
111     -0.999999993299972345915102
112     -0.999999992711294372006456
113     -0.999999992072439099985541
114     -0.999999991379268200130508
115     -0.999999990627314148132540
116     -0.999999989811755066458970
117     -0.999999988927387727619661
118     -0.999999987968598591791738
119     -0.999999986929332742926881
120     -0.999999985803060578628974
121     -0.999999984582742099719981
122     -0.999999983260788635479157
123     -0.999999981829021830014495
124     -0.999999980278629704073472
125     -0.999999978600119594788978
126     -0.999999976783267763350693
127     -0.999999974817065447355071
128     -0.999999972689661120580137
129     -0.999999970388298708114035
130     -0.999999967899251489096719
131     -0.999999965207751402768258
132     -0.999999962297913456009131
133     -0.999999959152654912059634
134     -0.999999955753608920567185
135     -0.999999952081032228479841
136     -0.999999948113706589527419
137     -0.999999943828833467051722
138     -0.999999939201921600705774
139     -0.999999934206666981977303
140     -0.999999928814824756540487
141     -0.999999922996072543035871
142     -0.999999916717864627952817
143     -0.999999909945276464770342
144     -0.999999902640838872326778
145     -0.999999894764361292459493
146     -0.999999886272743430203347
147     -0.999999877119774561178120
148     -0.999999867255919750145400
149     -0.999999856628092181986012
150     -0.999999845179410761448597
151     -0.999999832848942090853798
152     -0.999999819571425885409310
153     -0.999999805276982833797809
154     -0.999999789890803857138683
155     -0.999999773332819662188394
156     -0.999999755517349424622676
157     -0.999999736352727375323111
158     -0.999999715740905996653838
159     -0.999999693577034466641037
160     -0.999999669749010916634964
161     -0.999999644137006992314701
162     -0.999999616612963128659530
163     -0.999999587040052866624547
164     -0.999999555272114452585247
165     -0.999999521153047871016725
166     -0.999999484516175366204968
167     -0.999999445183563409904602
168     -0.999999402965303968610622
169     -0.999999357658752816349255
170     -0.999999309047722526460751
171     -0.999999256901627658587254
172     -0.999999200974579534832405
173     -0.999999141004427871663589
174     -0.999999076711746401418634
175     -0.999999007798759479089598
176     -0.999998933948206526218906
177     -0.999998854822141014087435
178     -0.999998770060660532728832
179     -0.999998679280564330496742
180     -0.999998582073934540768369
181     -0.999998478006637137714638
182     -0.999998366616738481729930
183     -0.999998247412833126918517
184     -0.999998119872278367806604
185     -0.999997983439330800015059
186     -0.999997837523179959816750
187     -0.999997681495873890143584
188     -0.999997514690131255533768
189     -0.999997336397034395554135
190     -0.999997145863597465233101
191     -0.999996942290203561838538
192     -0.999996724827904479777092
193     -0.999996492575576469327953
194     -0.999996244576925100210857
195     -0.999995979817332047487635
196     -0.999995697220536324878170
197     -0.999995395645142189112255
198     -0.999995073880945628324028
199     -0.999994730645071027620088
200     -0.999994364577909275721057
201     -0.999993974238848237905028
202     -0.999993558101786172297790
203     -0.999993114550418308799710
204     -0.999992641873286442566995
205     -0.999992138258581016945169
206     -0.999991601788684784070191
207     -0.999991030434446735009890
208     -0.999990422049174585335565
209     -0.999989774362333686430315
210     -0.999989084972939807716761
211     -0.999988351342632800403966
212     -0.999987570788417709415870
213     -0.999986740475059447000171
214     -0.999985857407116679281056
215     -0.999984918420600105892450
216     -0.999983920174239833019164
217     -0.999982859140346051919802
218     -0.999981731595246738576625
219     -0.999980533609285585815045
220     -0.999979261036362867394040
221     -0.999977909503001414558627
222     -0.999976474396919359773184
223     -0.999974950855090770264459
224     -0.999973333751274756079605
225     -0.999971617682993094132439
226     -0.999969796957935861737626
227     -0.999967865579774021028779
228     -0.999965817233357340078774
229     -0.999963645269275478192029
230     -0.999961342687759502470105
231     -0.999958902121900541164316
232     -0.999956315820161717373451
233     -0.999953575628158945224767
234     -0.999950672969685610748981
235     -0.999947598826955602236751
236     -0.999944343720038601012541
237     -0.999940897685460994406321
238     -0.999937250253945229426789
239     -0.999933390427259889483441
240     -0.999929306654152248770470
241     -0.999924986805334540975777
242     -0.999920418147494672237463
243     -0.999915587316301614224693
244     -0.999910480288375233419033
245     -0.999905082352189848729210
246     -0.999899378077880363163096
247     -0.999893351285919388145981
248     -0.999886985014633373018311
249     -0.999880261486525369137321
250     -0.999873162073371699777589
251     -0.999865667260059475670860
252     -0.999857756607131593604098
253     -0.999849408712005584127087
254     -0.999840601168832436289035
255     -0.999831310526961324669729
256     -0.999821512247975999096039
257     -0.999811180661268472697860
258     -0.999800288918115561773568
259     -0.999788808944223793772524
260     -0.999776711390708210081766
261     -0.999763965583470650796009
262     -0.999750539470943221871744
263     -0.999736399570162813678909
264     -0.999721510911142766670001
265     -0.999705836979508067426713
266     -0.999689339657360809492885
267     -0.999671979162343070964282
268     -0.999653713984864847611602
269     -0.999634500823465239215587
270     -0.999614294518275720661921
271     -0.999593047982555041060436
272     -0.999570712132266086606323
273     -0.999547235813665918980236
274     -0.999522565728881163674247
275     -0.999496646359441974605568
276     -0.999469419887748945596554
277     -0.999440826116448578590456
278     -0.999410802385694255639232
279     -0.999379283488271099505731
280     -0.999346201582564648884856
281     -0.999311486103354921430255
282     -0.999275063670419193575791
283     -0.999236857994928693095374
284     -0.999196789783626380902110
285     -0.999154776640775095111632
286     -0.999110732967867545150889
287     -0.999064569861091978842036
288     -0.999016195006549802974103
289     -0.998965512573224019809696
290     -0.998912423103700050040183
291     -0.998856823402643348534653
292     -0.998798606423041184302075
293     -0.998737661150219051712021
294     -0.998673872483645407328131
295     -0.998607121116541786642636
296     -0.998537283413318848302089
297     -0.998464231284862520615845
298     -0.998387832061698186549948
299     -0.998307948365064739107359
300     -0.998224437975934368814823
301     -0.998137153702018108556548
302     -0.998045943242801457529124
303     -0.997950649052658834678793
304     -0.997851108202100171372239
305     -0.997747152237207641699493
306     -0.997638607037325344858751
307     -0.997525292671069695356701
308     -0.997407023250733340776017
309     -0.997283606785160610778704
310     -0.997154845031177801648948
311     -0.997020533343667014496115
312     -0.996880460524377788338765
313     -0.996734408669576397404104
314     -0.996582153016638410456033
315     -0.996423461789695933775340
316     -0.996258096044456873132539
317     -0.996085809512319547895015
318     -0.995906348443912066070672
319     -0.995719451452192015978668
320     -0.995524849355248241312063
321     -0.995322265018952734162069
322     -0.995111413199616997238347
323     -0.994892000386813583757223
324     -0.994663724646529912575485
325     -0.994426275464827868000969
326     -0.994179333592189118776836
327     -0.993922570888732519066544
328     -0.993655650170496375504198
329     -0.993378225056984767759469
330     -0.993089939820183484438217
331     -0.992790429235257469948358
332     -0.992479318433147959270642
333     -0.992156222755293694076721
334     -0.991820747610706752658708
335     -0.991472488335639574649461
336     -0.991111030056085706155210
337     -0.990735947553362618000491
338     -0.990346805133030645297762
339     -0.989943156497407646205095
340     -0.989524544621944366953582
341     -0.989090501635730714183733
342     -0.988640548706408159263284
343     -0.988174195929768317289394
344     -0.987690942224322340437041
345     -0.987190275231130125566791
346     -0.986671671219182443772211
347     -0.986134594996632938939331
348     -0.985578499828180497431188
349     -0.985002827358905745666562
350     -0.984407007544868370574912
351     -0.983790458590774563626243
352     -0.983152586895026146341339
353     -0.982492787002464827943096
354     -0.981810441565126558280198
355     -0.981104921311322055650339
356     -0.980375585023360294162452
357     -0.979621779524232013515387
358     -0.978842839673570150500483
359     -0.978038088373203471420329
360     -0.977206836582618593554779
361     -0.976348383344644007774283
362     -0.975462015821667639794134
363     -0.974547009342696901235481
364     -0.973602627461567070965916
365     -0.972628122026600200033799
366     -0.971622733262012538372566
367     -0.970585689861363727448927
368     -0.969516209093335679948252
369     -0.968413496920123165931304
370     -0.967276748128711635913826
371     -0.966105146475310727066976
372     -0.964897864843204212102907
373     -0.963654065414268855166142
374     -0.962372899854405729886811
375     -0.961053509513118027313721
376     -0.959695025637459232377443
377     -0.958296569600564774513511
378     -0.956857253144968859850497
379     -0.955376178640896168788291
380     -0.953852439359705454847226
381     -0.952285119762648810516483
382     -0.950673295805096476827758
383     -0.949016035256362570111835
384     -0.947312398035251987699776
385     -0.945561436561433041323207
386     -0.943762196122724061065830
387     -0.941913715258365323657307
388     -0.940015026158330197754595
389     -0.938065155078711378050909
390     -0.936063122773199513379509
391     -0.934007944940652436603893
392     -0.931898632688733592554243
393     -0.929734193013578152514172
394     -0.927513629295424719100128
395     -0.925235941810129484510736
396     -0.922900128256458230136523
397     -0.920505184299029669863949
398     -0.918050104126761367892733
399     -0.915533881026646830667193
400     -0.912955507972669409081356
401     -0.910313978229635380238406
402     -0.907608285971685037790485
403     -0.904837426915216837567545
404     -0.902000398965935653353894
405     -0.899096202879712030000041
406     -0.896123842936915012419506
407     -0.893082327629856715369185
408     -0.889970670362962317254718
409     -0.886787890165254649549837
410     -0.883533012414718050732089
411     -0.880205069574081699771868
412     -0.876803101937538279856830
413     -0.873326158387889589649984
414     -0.869773297163586659132690
415     -0.866143586635108082774537
416     -0.862436106090096697765226
417     -0.858649946526651453191925
418     -0.854784211454148381612389
419     -0.850838017700942042052444
420     -0.846810496228276697886304
421     -0.842700792949714869341221
422     -0.838508069555369803579790
423     -0.834231504340207880514292
424     -0.829870293035667064651115
425     -0.825423649643818267103596
426     -0.820890807273277941907934
427     -0.816271018976062450999077
428     -0.811563558584557736234438
429     -0.806767721547761677691787
430     -0.801882825765941237429626
431     -0.796908212422832128518725
432     -0.791843246814495351087002
433     -0.786687319173932540042657
434     -0.781439845490550711651092
435     -0.776100268323556715663465
436     -0.770668057608352532380083
437     -0.765142711454994534663544
438     -0.759523756937772996814901
439     -0.753810750874962507155164
440     -0.748003280597789559230994
441     -0.742100964707660486167111
442     -0.736103453820691090374864
443     -0.730010431298578831244359
444     -0.723821613964859289254273
445     -0.717536752805590846718625
446     -0.711155633653515131598938
447     -0.704678077854745777673723
448     -0.698103942917044475821896
449     -0.691433123138751138842000
450     -0.684665550217444285287115
451     -0.677801193837418472975629
452     -0.670840062235077784356887
453     -0.663782202741357985424868
454     -0.656627702300305046439084
455     -0.649376687962954222750275
456     -0.642029327355671841292946
457     -0.634585829122141313859706
458     -0.627046443338195690464736
459     -0.619411461898721260562292
460     -0.611681218875880289504687
461     -0.603856090847925922562622
462     -0.595936497197908576518189
463     -0.587922900381600743987037
464     -0.579815806163996030236782
465     -0.571615763823768395668684
466     -0.563323366325108955991688
467     -0.554939250456390260497539
468     -0.546464096935141688222533
469     -0.537898630478854431106614
470     -0.529243619841170429010070
471     -0.520499877813046537682747
472     -0.511668261188523099211697
473     -0.502749670694764894533266
474     -0.493745050886082136517092
475     -0.484655390001679655026703
476     -0.475481719786923675319682
477     -0.466225115277957539319705
478     -0.456886694549540304987442
479     -0.447467618426025319819411
480     -0.437969090155439536373130
481     -0.428392355046668455103604
482     -0.418738700069796075759558
483     -0.409009453419694044918239
484     -0.399205984042999234582144
485     -0.389329701128664204753628
486     -0.379382053562310320297300
487     -0.369364529344658637204124
488     -0.359278654974358972988254
489     -0.349125994795582757315695
490     -0.338908150310790248699499
491     -0.328626759459127427638914
492     -0.318283495860952262338992
493     -0.307880068029034016649000
494     -0.297418218547012757893882
495     -0.286899723215749152048332
496     -0.276326390168236932985068
497     -0.265700058953792029795712
498     -0.255022599592273160770179
499     -0.244295911599128684869262
500     -0.233521922982103572925263
501     -0.222702589210478454140139
502     -0.211839892157749742936947
503     -0.200935839018695792648413
504     -0.189992461201808793395082
505     -0.179011813198105671155184
506     -0.167995971427363494770210
507     -0.156947033062855800871902
508     -0.145867114835695749082958
509     -0.134758351819920069295014
510     -0.123622896199474309769839
511     -0.112462916018284892203275
512     -0.101280593914626883352498
513     -0.0900781258410181607233922
514     -0.0788577197708907433569970
515     -0.0676215943933084420794315
516     -0.0563719777970166238312711
517     -0.0451111061451247520897491
518     -0.0338412223417354333022167
519     -0.0225645746918449442243659
520     -0.0112834155558496169159095
521     0.
522     0.0112834155558496169159095
523     0.0225645746918449442243659
524     0.0338412223417354333022167
525     0.0451111061451247520897491
526     0.0563719777970166238312711
527     0.0676215943933084420794315
528     0.0788577197708907433569970
529     0.0900781258410181607233922
530     0.101280593914626883352498
531     0.112462916018284892203275
532     0.123622896199474309769839
533     0.134758351819920069295014
534     0.145867114835695749082958
535     0.156947033062855800871902
536     0.167995971427363494770210
537     0.179011813198105671155184
538     0.189992461201808793395082
539     0.200935839018695792648413
540     0.211839892157749742936947
541     0.222702589210478454140139
542     0.233521922982103572925263
543     0.244295911599128684869262
544     0.255022599592273160770179
545     0.265700058953792029795712
546     0.276326390168236932985068
547     0.286899723215749152048332
548     0.297418218547012757893882
549     0.307880068029034016649000
550     0.318283495860952262338992
551     0.328626759459127427638914
552     0.338908150310790248699499
553     0.349125994795582757315695
554     0.359278654974358972988254
555     0.369364529344658637204124
556     0.379382053562310320297300
557     0.389329701128664204753628
558     0.399205984042999234582144
559     0.409009453419694044918239
560     0.418738700069796075759558
561     0.428392355046668455103604
562     0.437969090155439536373130
563     0.447467618426025319819411
564     0.456886694549540304987442
565     0.466225115277957539319705
566     0.475481719786923675319682
567     0.484655390001679655026703
568     0.493745050886082136517092
569     0.502749670694764894533266
570     0.511668261188523099211697
571     0.520499877813046537682747
572     0.529243619841170429010070
573     0.537898630478854431106614
574     0.546464096935141688222533
575     0.554939250456390260497539
576     0.563323366325108955991688
577     0.571615763823768395668684
578     0.579815806163996030236782
579     0.587922900381600743987037
580     0.595936497197908576518189
581     0.603856090847925922562622
582     0.611681218875880289504687
583     0.619411461898721260562292
584     0.627046443338195690464736
585     0.634585829122141313859706
586     0.642029327355671841292946
587     0.649376687962954222750275
588     0.656627702300305046439084
589     0.663782202741357985424868
590     0.670840062235077784356887
591     0.677801193837418472975629
592     0.684665550217444285287115
593     0.691433123138751138842000
594     0.698103942917044475821896
595     0.704678077854745777673723
596     0.711155633653515131598938
597     0.717536752805590846718625
598     0.723821613964859289254273
599     0.730010431298578831244359
600     0.736103453820691090374864
601     0.742100964707660486167111
602     0.748003280597789559230994
603     0.753810750874962507155164
604     0.759523756937772996814901
605     0.765142711454994534663544
606     0.770668057608352532380083
607     0.776100268323556715663465
608     0.781439845490550711651092
609     0.786687319173932540042657
610     0.791843246814495351087002
611     0.796908212422832128518725
612     0.801882825765941237429626
613     0.806767721547761677691787
614     0.811563558584557736234438
615     0.816271018976062450999077
616     0.820890807273277941907934
617     0.825423649643818267103596
618     0.829870293035667064651115
619     0.834231504340207880514292
620     0.838508069555369803579790
621     0.842700792949714869341221
622     0.846810496228276697886304
623     0.850838017700942042052444
624     0.854784211454148381612389
625     0.858649946526651453191925
626     0.862436106090096697765226
627     0.866143586635108082774537
628     0.869773297163586659132690
629     0.873326158387889589649984
630     0.876803101937538279856830
631     0.880205069574081699771868
632     0.883533012414718050732089
633     0.886787890165254649549837
634     0.889970670362962317254718
635     0.893082327629856715369185
636     0.896123842936915012419506
637     0.899096202879712030000041
638     0.902000398965935653353894
639     0.904837426915216837567545
640     0.907608285971685037790485
641     0.910313978229635380238406
642     0.912955507972669409081356
643     0.915533881026646830667193
644     0.918050104126761367892733
645     0.920505184299029669863949
646     0.922900128256458230136523
647     0.925235941810129484510736
648     0.927513629295424719100128
649     0.929734193013578152514172
650     0.931898632688733592554243
651     0.934007944940652436603893
652     0.936063122773199513379509
653     0.938065155078711378050909
654     0.940015026158330197754595
655     0.941913715258365323657307
656     0.943762196122724061065830
657     0.945561436561433041323207
658     0.947312398035251987699776
659     0.949016035256362570111835
660     0.950673295805096476827758
661     0.952285119762648810516483
662     0.953852439359705454847226
663     0.955376178640896168788291
664     0.956857253144968859850497
665     0.958296569600564774513511
666     0.959695025637459232377443
667     0.961053509513118027313721
668     0.962372899854405729886811
669     0.963654065414268855166142
670     0.964897864843204212102907
671     0.966105146475310727066976
672     0.967276748128711635913826
673     0.968413496920123165931304
674     0.969516209093335679948252
675     0.970585689861363727448927
676     0.971622733262012538372566
677     0.972628122026600200033799
678     0.973602627461567070965916
679     0.974547009342696901235481
680     0.975462015821667639794134
681     0.976348383344644007774283
682     0.977206836582618593554779
683     0.978038088373203471420329
684     0.978842839673570150500483
685     0.979621779524232013515387
686     0.980375585023360294162452
687     0.981104921311322055650339
688     0.981810441565126558280198
689     0.982492787002464827943096
690     0.983152586895026146341339
691     0.983790458590774563626243
692     0.984407007544868370574912
693     0.985002827358905745666562
694     0.985578499828180497431188
695     0.986134594996632938939331
696     0.986671671219182443772211
697     0.987190275231130125566791
698     0.987690942224322340437041
699     0.988174195929768317289394
700     0.988640548706408159263284
701     0.989090501635730714183733
702     0.989524544621944366953582
703     0.989943156497407646205095
704     0.990346805133030645297762
705     0.990735947553362618000491
706     0.991111030056085706155210
707     0.991472488335639574649461
708     0.991820747610706752658708
709     0.992156222755293694076721
710     0.992479318433147959270642
711     0.992790429235257469948358
712     0.993089939820183484438217
713     0.993378225056984767759469
714     0.993655650170496375504198
715     0.993922570888732519066544
716     0.994179333592189118776836
717     0.994426275464827868000969
718     0.994663724646529912575485
719     0.994892000386813583757223
720     0.995111413199616997238347
721     0.995322265018952734162069
722     0.995524849355248241312063
723     0.995719451452192015978668
724     0.995906348443912066070672
725     0.996085809512319547895015
726     0.996258096044456873132539
727     0.996423461789695933775340
728     0.996582153016638410456033
729     0.996734408669576397404104
730     0.996880460524377788338765
731     0.997020533343667014496115
732     0.997154845031177801648948
733     0.997283606785160610778704
734     0.997407023250733340776017
735     0.997525292671069695356701
736     0.997638607037325344858751
737     0.997747152237207641699493
738     0.997851108202100171372239
739     0.997950649052658834678793
740     0.998045943242801457529124
741     0.998137153702018108556548
742     0.998224437975934368814823
743     0.998307948365064739107359
744     0.998387832061698186549948
745     0.998464231284862520615845
746     0.998537283413318848302089
747     0.998607121116541786642636
748     0.998673872483645407328131
749     0.998737661150219051712021
750     0.998798606423041184302075
751     0.998856823402643348534653
752     0.998912423103700050040183
753     0.998965512573224019809696
754     0.999016195006549802974103
755     0.999064569861091978842036
756     0.999110732967867545150889
757     0.999154776640775095111632
758     0.999196789783626380902110
759     0.999236857994928693095374
760     0.999275063670419193575791
761     0.999311486103354921430255
762     0.999346201582564648884856
763     0.999379283488271099505731
764     0.999410802385694255639232
765     0.999440826116448578590456
766     0.999469419887748945596554
767     0.999496646359441974605568
768     0.999522565728881163674247
769     0.999547235813665918980236
770     0.999570712132266086606323
771     0.999593047982555041060436
772     0.999614294518275720661921
773     0.999634500823465239215587
774     0.999653713984864847611602
775     0.999671979162343070964282
776     0.999689339657360809492885
777     0.999705836979508067426713
778     0.999721510911142766670001
779     0.999736399570162813678909
780     0.999750539470943221871744
781     0.999763965583470650796009
782     0.999776711390708210081766
783     0.999788808944223793772524
784     0.999800288918115561773568
785     0.999811180661268472697860
786     0.999821512247975999096039
787     0.999831310526961324669729
788     0.999840601168832436289035
789     0.999849408712005584127087
790     0.999857756607131593604098
791     0.999865667260059475670860
792     0.999873162073371699777589
793     0.999880261486525369137321
794     0.999886985014633373018311
795     0.999893351285919388145981
796     0.999899378077880363163096
797     0.999905082352189848729210
798     0.999910480288375233419033
799     0.999915587316301614224693
800     0.999920418147494672237463
801     0.999924986805334540975777
802     0.999929306654152248770470
803     0.999933390427259889483441
804     0.999937250253945229426789
805     0.999940897685460994406321
806     0.999944343720038601012541
807     0.999947598826955602236751
808     0.999950672969685610748981
809     0.999953575628158945224767
810     0.999956315820161717373451
811     0.999958902121900541164316
812     0.999961342687759502470105
813     0.999963645269275478192029
814     0.999965817233357340078774
815     0.999967865579774021028779
816     0.999969796957935861737626
817     0.999971617682993094132439
818     0.999973333751274756079605
819     0.999974950855090770264459
820     0.999976474396919359773184
821     0.999977909503001414558627
822     0.999979261036362867394040
823     0.999980533609285585815045
824     0.999981731595246738576625
825     0.999982859140346051919802
826     0.999983920174239833019164
827     0.999984918420600105892450
828     0.999985857407116679281056
829     0.999986740475059447000171
830     0.999987570788417709415870
831     0.999988351342632800403966
832     0.999989084972939807716761
833     0.999989774362333686430315
834     0.999990422049174585335565
835     0.999991030434446735009890
836     0.999991601788684784070191
837     0.999992138258581016945169
838     0.999992641873286442566995
839     0.999993114550418308799710
840     0.999993558101786172297790
841     0.999993974238848237905028
842     0.999994364577909275721057
843     0.999994730645071027620088
844     0.999995073880945628324028
845     0.999995395645142189112255
846     0.999995697220536324878170
847     0.999995979817332047487635
848     0.999996244576925100210857
849     0.999996492575576469327953
850     0.999996724827904479777092
851     0.999996942290203561838538
852     0.999997145863597465233101
853     0.999997336397034395554135
854     0.999997514690131255533768
855     0.999997681495873890143584
856     0.999997837523179959816750
857     0.999997983439330800015059
858     0.999998119872278367806604
859     0.999998247412833126918517
860     0.999998366616738481729930
861     0.999998478006637137714638
862     0.999998582073934540768369
863     0.999998679280564330496742
864     0.999998770060660532728832
865     0.999998854822141014087435
866     0.999998933948206526218906
867     0.999999007798759479089598
868     0.999999076711746401418634
869     0.999999141004427871663589
870     0.999999200974579534832405
871     0.999999256901627658587254
872     0.999999309047722526460751
873     0.999999357658752816349255
874     0.999999402965303968610622
875     0.999999445183563409904602
876     0.999999484516175366204968
877     0.999999521153047871016725
878     0.999999555272114452585247
879     0.999999587040052866624547
880     0.999999616612963128659530
881     0.999999644137006992314701
882     0.999999669749010916634964
883     0.999999693577034466641037
884     0.999999715740905996653838
885     0.999999736352727375323111
886     0.999999755517349424622676
887     0.999999773332819662188394
888     0.999999789890803857138683
889     0.999999805276982833797809
890     0.999999819571425885409310
891     0.999999832848942090853798
892     0.999999845179410761448597
893     0.999999856628092181986012
894     0.999999867255919750145400
895     0.999999877119774561178120
896     0.999999886272743430203347
897     0.999999894764361292459493
898     0.999999902640838872326778
899     0.999999909945276464770342
900     0.999999916717864627952817
901     0.999999922996072543035871
902     0.999999928814824756540487
903     0.999999934206666981977303
904     0.999999939201921600705774
905     0.999999943828833467051722
906     0.999999948113706589527419
907     0.999999952081032228479841
908     0.999999955753608920567185
909     0.999999959152654912059634
910     0.999999962297913456009131
911     0.999999965207751402768258
912     0.999999967899251489096719
913     0.999999970388298708114035
914     0.999999972689661120580137
915     0.999999974817065447355071
916     0.999999976783267763350693
917     0.999999978600119594788978
918     0.999999980278629704073472
919     0.999999981829021830014495
920     0.999999983260788635479157
921     0.999999984582742099719981
922     0.999999985803060578628974
923     0.999999986929332742926881
924     0.999999987968598591791738
925     0.999999988927387727619661
926     0.999999989811755066458970
927     0.999999990627314148132540
928     0.999999991379268200130508
929     0.999999992072439099985541
930     0.999999992711294372006456
931     0.999999993299972345915102
932     0.999999993842305597079599
933     0.999999994341842780638679
934     0.999999994801868964843021
935     0.999999995225424562376994
936     0.999999995615322952246180
937     0.999999995974166879001692
938     0.999999996304363710601733
939     0.999999996608139631065173
940     0.999999996887552839233372
941     0.999999997144505820407811
942     0.999999997380756753356393
943     0.999999997597930111164948
944     0.999999997797526510638154
945     0.999999997980931861411821
946     0.999999998149425862613257
947     0.999999998304189891785759
948     0.999999998446314327865434
949     0.999999998576805347252390
950     0.999999998696591229443244
951     0.999999998806528206277959
952     0.999999998907405886591731
953     0.999999998999952285943071
954     0.999999999084838489103913
955     0.999999999162682971138542
956     0.999999999234055601157792
957     0.999999999299481351206210
958     0.999999999359443731215952
959     0.999999999414387969535759
960     0.999999999464723957210342
961     0.999999999510828972939411
962     0.999999999553050204480830
963     0.999999999591707081174197
964     0.999999999627093431244556
965     0.999999999659479476596644
966     0.999999999689113676923764
967     0.999999999716224434128131
968     0.999999999741021667277623
969     0.999999999763698267603729
970     0.999999999784431442373954
971     0.999999999803383955845711
972     0.999999999820705274925103
973     0.999999999836532626610029
974     0.999999999850991973790366
975     0.999999999864198915505960
976     0.999999999876259517323697
977     0.999999999887271077085832
978     0.999999999897322830900966
979     0.999999999906496603894840
980     0.999999999914867409908580
981     0.999999999922504004025582
982     0.999999999929469391523360
983     0.999999999935821296581860
984     0.999999999941612593833703
985     0.999999999946891705613252
986     0.999999999951702967549089
987     0.999999999956086964947412
988     0.999999999960080842230859
989     0.999999999963718587527495
990     0.999999999967031294347156
991     0.999999999970047402136203
992     0.999999999972792917366275
993     0.999999999975291616686970
994     0.999999999977565233555991
995     0.999999999979633629652374
996     0.999999999981514952278515
997     0.999999999983225778864110
998     0.999999999984781249599477
999     0.999999999986195189146334
1000     0.999999999987480218300757
1001     0.999999999988647856415078
1002     0.999999999989708615322670
1003     0.999999999990672085451468
1004     0.999999999991547014758366
1005     0.999999999992341381066975
1006     0.999999999993062458345374
1007     0.999999999993716877418075
1008     0.999999999994310681567322
1009     0.999999999994849377442675
1010     0.999999999995337981664472
1011     0.999999999995781063475994
1012     0.999999999996182783770717
1013     0.999999999996546930794861
1014     0.999999999996876952801260
1015     0.999999999997175987908301
1016     0.999999999997446891397149
1017     0.999999999997692260661541
1018     0.999999999997914458007001
1019     0.999999999998115631480261
1020     0.999999999998297733894900
1021     0.999999999998462540205572]';
1022 endfunction
1023
1024 [Xr,Yr]=erf_ref();
1025 y=erf(Xr);
1026 assert_checkalmostequal(Yr,y,sqrt(%eps),sqrt(%eps));
1027
1028 // Tests are from faddeeva code:
1029
1030 assert_checkalmostequal(erf(1 + 2*%i), -0.5366435657785650339917955593141927494421 + -5.049143703447034669543036958614140565553*%i, 1.e-10);
1031 assert_checkalmostequal(erf(-1 + 2*%i), 0.5366435657785650339917955593141927494421 + -5.049143703447034669543036958614140565553*%i, 1.e-10);
1032 assert_checkalmostequal(erf(1 + -2*%i), -0.5366435657785650339917955593141927494421 + 5.049143703447034669543036958614140565553*%i, 1.e-10);
1033 assert_checkalmostequal(erf(-1 + -2*%i), 0.5366435657785650339917955593141927494421 + 5.049143703447034669543036958614140565553*%i, 1.e-10);
1034 assert_checkalmostequal(erf(9 + -28*%i), 0.3359473673830576996788000505817956637777e304 + -0.1999896139679880888755589794455069208455e304*%i, 1.e-10);
1035 assert_checkalmostequal(erf(21 + -33*%i), 0.3584459971462946066523939204836760283645e278 + 0.3818954885257184373734213077678011282505e280*%i, 1.e-10);
1036 assert_checkalmostequal(erf(1e3 + 1e3*%i), 0.9996020422657148639102150147542224526887 + 0.00002801044116908227889681753993542916894856*%i, 1.e-10);
1037 assert_checkalmostequal(erf(-3001 + -1000*%i), -1 + 0*%i, 1.e-10);
1038 assert_checkalmostequal(erf(1e160 + -1e159*%i), 1 + 0*%i, 1.e-10);
1039 assert_checkalmostequal(erf(5.1e-3 + 1e-8*%i), 0.005754683859034800134412990541076554934877 + 0.1128349818335058741511924929801267822634e-7*%i, 1.e-10);
1040 assert_checkalmostequal(erf(-4.9e-3 + 4.95e-3*%i), -0.005529149142341821193633460286828381876955 + 0.005585388387864706679609092447916333443570*%i, 1.e-10);
1041 assert_checkalmostequal(erf(4.9e-3 + 0.5*%i), 0.007099365669981359632319829148438283865814 + 0.6149347012854211635026981277569074001219*%i, 1.e-10);
1042 assert_checkalmostequal(erf(4.9e-4 + -0.5e1*%i), 0.3981176338702323417718189922039863062440e8 + -0.8298176341665249121085423917575122140650e10*%i, 1.e-10);
1043 assert_checkalmostequal(erf(-4.9e-5 + -0.5e2*%i), complex(-%inf, -%inf), 1.e-10);
1044 assert_checkalmostequal(erf(5.1e-3 + 0.5*%i), 0.007389128308257135427153919483147229573895 + 0.6149332524601658796226417164791221815139*%i, 1.e-10);
1045 assert_checkalmostequal(erf(5.1e-4 + -0.5e1*%i), 0.4143671923267934479245651547534414976991e8 + -0.8298168216818314211557046346850921446950e10*%i, 1.e-10);
1046 assert_checkalmostequal(erf(-5.1e-5 + -0.5e2*%i), complex(-%inf, -%inf), 1.e-10);
1047 assert_checkalmostequal(erf(1e-6 + 2e-6*%i), 0.1128379167099649964175513742247082845155e-5 + 0.2256758334191777400570377193451519478895e-5*%i, 1.e-10);
1048 assert_checkalmostequal(erf(0 + 2e-6*%i), 0 + 0.2256758334194034158904576117253481476197e-5*%i, 1.e-10);
1049 assert_checkalmostequal(erf(0 + 2*%i), 0 + 18.56480241457555259870429191324101719886*%i, 1.e-10);
1050 assert_checkalmostequal(erf(0 + 20*%i), 0 + 0.1474797539628786202447733153131835124599e173*%i, 1.e-10);
1051 assert_checkalmostequal(erf(0 + 200*%i), complex(0, %inf), 1.e-10);
1052 assert_checkalmostequal(erf(%inf + 0*%i), 1 + 0*%i, 1.e-10);
1053 assert_checkalmostequal(erf(-%inf + 0*%i), -1 + 0*%i, 1.e-10);
1054 assert_checkalmostequal(erf(0 + %inf*%i), complex(%nan, %nan), 1.e-10);
1055 assert_checkalmostequal(erf(0 + -%inf*%i), complex(%nan, %nan), 1.e-10);
1056 assert_checkalmostequal(erf(%inf + %inf*%i), %nan + %nan*%i, 1.e-10);
1057 assert_checkalmostequal(erf(%inf + -%inf*%i), %nan + %nan*%i, 1.e-10);
1058 assert_checkalmostequal(erf(%nan + %nan*%i), %nan + %nan*%i, 1.e-10);
1059 assert_checkalmostequal(erf(%nan + 0*%i), %nan + 0*%i, 1.e-10);
1060 assert_checkalmostequal(erf(0 + %nan*%i), 0 + %nan*%i, 1.e-10);
1061 assert_checkalmostequal(erf(%nan + %inf*%i), %nan + %nan*%i, 1.e-10);
1062 assert_checkalmostequal(erf(%inf + %nan*%i), %nan + %nan*%i, 1.e-10);
1063 assert_checkalmostequal(erf(1e-3 + %nan*%i), %nan + %nan*%i, 1.e-10);
1064 assert_checkalmostequal(erf(7e-2 + 7e-2*%i), 0.07924380404615782687930591956705225541145 + 0.07872776218046681145537914954027729115247*%i, 1.e-10);
1065 assert_checkalmostequal(erf(7e-2 + -7e-4*%i), 0.07885775828512276968931773651224684454495 + -0.0007860046704118224342390725280161272277506*%i, 1.e-10);
1066 assert_checkalmostequal(erf(-9e-2 + 7e-4*%i), -0.1012806432747198859687963080684978759881 + 0.0007834934747022035607566216654982820299469*%i, 1.e-10);
1067 assert_checkalmostequal(erf(-9e-2 + 9e-2*%i), -0.1020998418798097910247132140051062512527 + 0.1010030778892310851309082083238896270340*%i, 1.e-10);
1068 assert_checkalmostequal(erf(-7e-4 + 9e-2*%i), -0.0007962891763147907785684591823889484764272 + 0.1018289385936278171741809237435404896152*%i, 1.e-10);
1069 assert_checkalmostequal(erf(7e-2 + 0.9e-2*%i), 0.07886408666470478681566329888615410479530 + 0.01010604288780868961492224347707949372245*%i, 1.e-10);
1070 assert_checkalmostequal(erf(7e-2 + 1.1e-2*%i), 0.07886723099940260286824654364807981336591 + 0.01235199327873258197931147306290916629654*%i, 1.e-10);
1071
1072 assert_checkerror("erf(''foo'')",[],999);
1073 assert_checkerror("erf(sparse([1,2;4,5;3,10],[1,2,3]))",[],999);