require 'mathlib'
function getw(q, n) return kvad(1, 1/q^2-2, 1-n) end
local q = tonumber(arg[1])
print(q, 'pyco')
print(getw(q, 2))
term% draws -A 1 -Q 2 1 A/Q 1 / 1 Q/A 1 | gview
term% draws -A 1 -Q 2 1 A/Q 1 / 1 Q/A 1 | gview
term% draws -A 1/4 -Q 4 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 4 -Q 4 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 1 -Q 4 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 1/2 -Q 4 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 1/2 -Q 4 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 1/4 -Q 4 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 1/8 -Q 4 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 1/8 -Q 4 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 1/8 -Q 8 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 1/8 -Q 8 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 1/8 -Q 8 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 1/8 -Q 8 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 1/8 -Q 256 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 1/2 -Q 256 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 1/4 -Q 256 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 4 -Q 256 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 4 -Q 4256 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 1/4 -Q 4256 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 1/4 -Q 4256 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 1/4 -Q 4 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 1/4 -Q 4 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% ""
	 term% draws -A 1/4 -Q 4 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
rc: line 2: token '( [SUB]': syntax error
term% draws -A 1/4 -Q 4 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draws -A 1/4 -Q 4 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
term% draw 0 4 'db((1/4)^2*(x^4+x^2*(1/4*1/16-2*1/4)+1/16)/(1/16*x^4+(1/4*1/16-2*1/4)*x^2+1))' | gview
term% { draws -A 1/4 -Q 4 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1; draw 0 4 'db((1/4)^2*(x^4+x^2*(1/4*1/16-2*1/4)+1/16)/(1/16*x^4+(1/4*1/16-2*1/4)*x^2+1))' } | gview
term% { draws -A 1/4 -Q 4 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1; draw 0 4 'db((1/4)^2*(x^4+x^2*(1/4*1/16-2*1/4)+1/16)/(1/16*x^4+(1/4*1/16-2*1/4)*x^2+1))' } | gview
term% draws -A 1/4 -Q 4 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1 | gview
{ draws -A 1/4 -Q 4 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1; draw 0 4 'db(1/4^2*(x^2-1/4)^2/(x^2/4-1)^2)' } | gview
{ draws -A 1/4 -Q 2 -c 'sqrt(A)/Q' A c 1 / 1/A c/A 1; draw 0 4 'db(1/4^2*(x^2-1/4)^2/(x^2/4-1)^2)' } | trim -128 128 | gview
drawz -w 'pi/4' -c '1/tan(w/2)' -d '1/(1+c)' d d / 1 '(1-c)*d' | trim -40 40 | gview