gsl-shell.git - gsl-shell

diff --git a/lua-lmtest.lua b/lua-lmtest.lua
deleted file mode 100644
index 46c8944d..00000000
--- a/lua-lmtest.lua
+++ /dev/null

@@ -1,37 +0,0 @@

-

-local sin, cos, exp, sqrt = math.sin, math.cos, math.exp, math.sqrt

-

-local n = 40

-

-local yrf, sigrf

-

-local fdf = function(x, f, J)

- for i=1, n do

- local A, lambda, b = x[1], x[2], x[3]

- local t, y, sig = i-1, yrf[i], sigrf[i]

- local e = exp(- lambda * t)

- if f then f[i] = (A*e+b - y)/sig end

- if J then

- J:set(i, 1, e / sig)

- J:set(i, 2, - t * A * e / sig)

- J:set(i, 3, 1 / sig)

- end

- end

- end

-

-local A, lambda, b = 5, 0.1, 1

-local r = gsl.rng('mt19937')

-r:set(0)

-

-yrf = gsl.new(n, 1, function(i) return A * exp(-lambda*(i-1)) + b + gsl.rnd.gaussian(r, 0.1) end)

-sigrf = gsl.new(n, 1, function() return 0.1 end)

-

-local s = gsl.nlinfit {n= n, p= 3}

-

-s:set(fdf, gsl.vector {1, 0, 0})

-print(gsl.tr(s.x), s.chisq)

-

-for i=1, 10 do

- s:iterate()

- print('ITER=', i, ': ', gsl.tr(s.x), s.chisq)

-end