diff options
| author | lneto <lneto@NetBSD.org> | 2014-07-19 18:38:33 +0000 |
|---|---|---|
| committer | lneto <lneto@NetBSD.org> | 2014-07-19 18:38:33 +0000 |
| commit | c4b1eb35e827c43a0dcf41027c98aaedf46be036 (patch) | |
| tree | f80015954e9c51266af4f2e1c57f2134cbcc1d62 /external/mit/lua/dist/test/bisect.lua | |
| parent | c6d74635a1d93d95a3770bd621dcd83ed33100c2 (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.lua | 27 |
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) |
