On 7/20/2020 7:56 AM, Jeong Dal wrote:
> Dear Fabrice,
>
> You may split Binom(n,k) function into two functions as following:
>
> see original mail
>
> Dalyoung
Best stay in a protected namespace ...
\startluacode
local function fact (n)
if n <= 0 then
return 1
else
return n * fact(n-1)
end
end
local function ncr(n,r)
return fact(n)/(fact(r)*fact(n-r))
end
userdata.P = {
fact = fact,
ncr = ncr,
}
function MP.pascal_ncr(n, r)
mp.print(ncr(n,r))
end
\stopluacode
Watch the last definition. This permits
% tt := lua("mp.print(userdata.P.ncr(" & decimal n & "," &
decimal r & " ))");
replaced by
tt := lua.MP.pascal_ncr(n,r);
which looks nicer.
\startbuffer[pt1]
numeric n, r, s, u, dx, dy, tt;
path p, q;
pair A, B, start, now;
u := 1.8cm;
A := dir(210)*u;
B := dir(-30)*u;
dy := sind(30)*u;
dx := 2*cosd(30)*u;
for n=0 upto 4:
start := n*dir(210)*u;
for r=0 upto n:
s := n-r;
% tt := lua("mp.print(userdata.P.ncr(" & decimal n & "," &
decimal r & " ))");
tt := lua.MP.pascal_ncr(n,r);
now := start+r*right*dx;
dotlabel.top(textext("$\displaystyle {" & decimal n &
"\choose" & decimal r & "} = "& decimal tt & "$"),now);
draw (now+A) -- now -- (now+B);
endfor;
endfor;
\stopbuffer
Now, in context lmtx we can have a different kind of abstraction. We can
do this:
function MP.pascal_ncr_x()
mp.print(ncr(mp.scan.pair()))
end
and then use:
tt := runscript("MP.pascal_ncr_x()") (n,r) ;
Of course one can decide to pick to two numerics instead, like
tt := runscript("MP.pascal_ncr_x()") n r ;
but i leave that as exercise.
% tt := runscript mp_pascal_ncr (n,r) ;
tt := pascal_ncr (n,r) ;
However, we still have the rather verbose runscript here, so we go
further, we register pascal as script:
\startluacode
metapost.registerscript("pascal_ncr",MP.pascal_ncr_x)
\stopluacode
And then define an alias at the metafun end:
\startMPextensions
newinternal mp_pascal_ncr ; mp_pascal_ncr := scriptindex "pascal_ncr" ;
def pascal_ncr =
runscript mp_pascal_ncr
enddef ;
\stopMPextensions
The internal permits this:
tt := runscript mp_pascal_ncr (n,r) ;
while the additional def permits
tt := pascal_ncr (n,r) ;
Now watch out, because we define pascal_ncr here, something
lua.MP.pascal_ncr(n,r) won't work because the last part gets expanded
because that is what mp does (i'll probably cook something for that some
day).
Now, to come back to
"I couldn’t wikify it at that time because I don’t know
how to. I’ll do it soon."
looks like you suddenly have an additional challenge,
Hans
-----------------------------------------------------------------
Hans Hagen | PRAGMA ADE
Ridderstraat 27 | 8061 GH Hasselt | The Netherlands
tel: 038 477 53 69 | www.pragma-ade.nl | www.pragma-pod.nl
-----------------------------------------------------------------
___________________________________________________________________________________
If your question is of interest to others as well, please add an entry to the Wiki!
maillist : ntg-context@ntg.nl / http://www.ntg.nl/mailman/listinfo/ntg-context
webpage : http://www.pragma-ade.nl / http://context.aanhet.net
archive : https://bitbucket.org/phg/context-mirror/commits/
wiki : http://contextgarden.net
___________________________________________________________________________________