program hilbert matrix define `1'=J(`2',`2',0) plugin call hilbertmat, `1' `2' end program hilbertmat, plugin