summaryrefslogtreecommitdiff
path: root/external/mit/lua/dist/test/bisect.lua
diff options
context:
space:
mode:
authorlneto <lneto@NetBSD.org>2014-07-19 18:38:33 +0000
committerlneto <lneto@NetBSD.org>2014-07-19 18:38:33 +0000
commitc4b1eb35e827c43a0dcf41027c98aaedf46be036 (patch)
treef80015954e9c51266af4f2e1c57f2134cbcc1d62 /external/mit/lua/dist/test/bisect.lua
parentc6d74635a1d93d95a3770bd621dcd83ed33100c2 (diff)
lua: updated from 5.1 to 5.3 work3
* lua(1): - changed lua_Integer to intmax_t - updated distrib/sets/lists and etc/mtree - updated bsd.lua.mk - fixed bozohttpd (lua-bozo.c) - compatibilized bindings: gpio, sqlite * lua(4): - removed floating-point and libc dependencies using '#ifndef _KERNEL' - fixed division by zero and exponentiation - libkern: added isalnum(), iscntrl(), isgraph(), isprint() and ispunct() - acpica: removed isprint() from acnetbsd.h - libc: moved strcspn.c, strpbrk.c and strspn.c to common - removed stub headers - compatibilized bindings: luapmf, luasystm * reorganized luaconf.h * updated doc/CHANGES and doc/RESPONSIBLE
Diffstat (limited to 'external/mit/lua/dist/test/bisect.lua')
-rw-r--r--external/mit/lua/dist/test/bisect.lua27
1 files changed, 0 insertions, 27 deletions
diff --git a/external/mit/lua/dist/test/bisect.lua b/external/mit/lua/dist/test/bisect.lua
deleted file mode 100644
index f91e69bfbaf..00000000000
--- a/external/mit/lua/dist/test/bisect.lua
+++ /dev/null
@@ -1,27 +0,0 @@
--- bisection method for solving non-linear equations
-
-delta=1e-6 -- tolerance
-
-function bisect(f,a,b,fa,fb)
- local c=(a+b)/2
- io.write(n," c=",c," a=",a," b=",b,"\n")
- if c==a or c==b or math.abs(a-b)<delta then return c,b-a end
- n=n+1
- local fc=f(c)
- if fa*fc<0 then return bisect(f,a,c,fa,fc) else return bisect(f,c,b,fc,fb) end
-end
-
--- find root of f in the inverval [a,b]. needs f(a)*f(b)<0
-function solve(f,a,b)
- n=0
- local z,e=bisect(f,a,b,f(a),f(b))
- io.write(string.format("after %d steps, root is %.17g with error %.1e, f=%.1e\n",n,z,e,f(z)))
-end
-
--- our function
-function f(x)
- return x*x*x-x-1
-end
-
--- find zero in [1,2]
-solve(f,1,2)