From 2d1df4714e2736dbde7855ddcd76b4c1de822fa5 Mon Sep 17 00:00:00 2001
From: David Walter Seikel
Date: Mon, 23 Jan 2012 21:58:02 +1000
Subject: Added a big bunch of example lua scripts for testing the speed of lua
 compiling.

---
 .../testLua/yueliang-0.4.1/test_lua/5.1/bisect.lua | 27 ++++++++++++++++++++++
 1 file changed, 27 insertions(+)
 create mode 100644 LuaSL/testLua/yueliang-0.4.1/test_lua/5.1/bisect.lua

(limited to 'LuaSL/testLua/yueliang-0.4.1/test_lua/5.1/bisect.lua')

diff --git a/LuaSL/testLua/yueliang-0.4.1/test_lua/5.1/bisect.lua b/LuaSL/testLua/yueliang-0.4.1/test_lua/5.1/bisect.lua
new file mode 100644
index 0000000..f91e69b
--- /dev/null
+++ b/LuaSL/testLua/yueliang-0.4.1/test_lua/5.1/bisect.lua
@@ -0,0 +1,27 @@
+-- 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)
-- 
cgit v1.1