©1jul00/dburkett@infinet.com
©Save polynomials in list bpoly[].
local k,k1,k2,bpdim,ff1,usermode,choice
©Save user's mode; set modes
getmode("all")→usermode
setmode({"Complex Format","Real","Vector
Format","Rectangular","Exact/Approx","Exact","Pretty Print","Off"})
©If bpoly[] exists unarchive it, else create it & initialize it.
if gettype(spfn\bpoly)="NONE" then
newlist(1)→spfn\bpoly
x-1/2→spfn\bpoly[1]
else
unarchiv(spfn\bpoly)
endif
dim(spfn\bpoly)→bpdim
clrio
©Loop to derive Bernoulli polynomials
if nxx>bpdim then
augment(spfn\bpoly,newlist(nxx-bpdim))→spfn\bpoly
for k,bpdim+1,nxx
á(k*spfn\bpoly[k-1],x)→ ff1
ff1-á(ff1,x,0,1)→spfn\bpoly[k]
disp k
disp spfn\bpoly[k]
endfor
endif
©Prompt to archive bpoly[]
1→choice
dialog
title "BNPOLY"
dropdown "Archive polynomial list?",{"yes","no"},choice
enddlog
if ok=1 and choice=1 then
archive spfn\bpoly
endif
©Restore user's modes
setmode(usermode)
clrio
disphome
Endprgm
Then, use bpo() to recall a particular polynomial:
bpo(kk)
func
©(k) Bernoulli polynomial Bk
©1jul00/dburkett@infinet.com
©Generates "Domain error" if kk>dim(bpoly)
if kk<0:return "bpo arg error"
when(kk=0,1,spfn\bpoly[kk])
6 - 42
To Next Page
Loading...
Need help?
General
