Raul Miller wrote:
> That seems different from what I know as the "vandermonde matrix".
OK, you're right: it's a Vandermonde matrix multiplied by a geometric series.
> I am not quite sure how to use your approach with this kind
> of polynomial.
>
I have returned to just evaluating polynomials in one variable.
ev1=:(>@{. p."1 0 {:@>@{:) ; }:@>@{:
eval1=:([: >@{. ev1^:([EMAIL PROTECTED]@>@{:)^:_) @ ;
arr&eval1"1 |: i. 3 9
27 49 69 81 79 57 9 _71 _189
Of course, with the boxing I am using, sparese arrsys will not work.
Best wishes,
John
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm