if x=0:return 0
if x<0:return
⁻
spfn\si(
⁻
x)
if x>Œ then
Œ/2-(((polyeval({
1
,0,38.027264,0,265.
1
87033,0,335.677320,0,38.
1
02495},x))/(polye
val({
1
,0,40.02
1
433,0,322.6249
11
,0,570.23628,0,
1
57.
1
05423},x)))/(x)*cos(x)+(polye
val({
1
,0,42.242855,0,302.757865,0,352.0
1
8498,0,2
1
.82
1
899},x)/polyeval({
1
,0,48.
1
9
6927,0,482.485984,0,
111
4.978885,0,449.690326},x))/(x^2)*sin(x))
elseif x>.0
1
then
nint(sin(z)/z,z,0,x)
else
polyeval({
⁻
2.2774643986765
⁻
9,0,3.06
1
9243582207
⁻
7,0,
⁻
2.834467
1
20
1
8
1
4
⁻
5,0,
1
.666
6666666667
⁻
3,0,
⁻
.055555555555556,0,
1
,0},x)
endif
EndFunc
Code listing for ci(x):
ci(x)
Func
©(x) return ci(x), x real, |x|<
11
2
©Must be in folder \spfn
©23jan0
1
if x=0:return undef
if x<0:return spfn\ci(
⁻
x)+
*Œ
if x>Œ then
(((polyeval({
1
,0,38.027264,0,265.
1
87033,0,335.677320,0,38.
1
02495},x))/(polyeval(
{
1
,0,40.02
1
433,0,322.6249
11
,0,570.23628,0,
1
57.
1
05423},x)))/(x))*sin(x)-((polyeva
l({
1
,0,42.242855,0,302.757865,0,352.0
1
8498,0,2
1
.82
1
899},x)/polyeval({
1
,0,48.
1
969
27,0,482.485984,0,
111
4.978885,0,449.690326},x))/(x^2))*cos(x)
elseif x>.0
1
then
nint((cos(t)-
1
)/t,t,0,x)+ln(x)+.5772
1
566490
1
53
else
ln(x)+.5772
1
566490
1
53+polyeval({3.
1
00
1
984
1
2698
⁻
6,0,
⁻
2.3
1
48
1
48
1
48
1
48
⁻
4,0,.0
1
04
1
6666666667,0,
⁻
.25,0,0},x)
endif
EndFunc
(Thanks to Bhuvanesh for providing test data.)
[6.34] Error function for real arguments
The error function is one of many so-called 'special functions'. It is defined as
6 - 64
Summary of Contents for TI-92+
Page 52: ...Component side of PCB GraphLink I O connector detail 1 41...
Page 53: ...LCD connector detail PCB switch side 1 42...
Page 54: ...Key pad sheet contact side Key pad sheet key side 1 43...
Page 55: ...Key cap detail 1 44...
Page 57: ...Component side of PCB with shield removed A detail view of the intergrated circuits 1 46...
Page 410: ...void extensionroutine2 void Credit to Bhuvanesh Bhatt 10 4...