lspoly fails with higher order polynomials

Hi,

I'm trying to run lspoly to get a good polynomial interpolation. Let's
consider this example :

------------------------
begin
 x = (/0.,1.,2.,3.,4.,5.,6./)
 y = (/0.,1.,2.,3.,4.,5.,6./)
 c = lspoly(x,y,1,5)
 print(c)
end
------------------------

which returns obviously the output (y = x) :

------------------------
Variable: c
Type: float
Total Size: 20 bytes
            5 values
Number of Dimensions: 1
Dimensions and sizes: [5]
Coordinates:
(0) 0
(1) 1
(2) 0
(3) 0
(4) 0
------------------------

But, if I just change the order of the requested polynomial from 4 to 5 :

------------------------
 c = lspoly(x,y,1,6)
------------------------

I get the result :

------------------------
Dimensions and sizes: [6]
Coordinates:
(0) 1e+20
(1) 1e+20
(2) 1e+20
(3) 1e+20
(4) 1e+20
(5) 1e+20
------------------------

Is this a bug that could be fixed soon ? I would really need a higher
order polynomial.

Thanks,
François

-- 
International Arctic Research Center
930, Koyukuk Drive
Fairbanks, AK, 99775-7335
USA
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk