I'm getting some odd output from the lgamma function: it seems to be oddly re-setting itself when fetching the log-gamma function values of individual numbers: perldl> p lgamma(60) 184 1 perldl> p lgamma(80) 13 1 perldl> p lgamma(70) 226 1 perldl> p lgamma(75) 247 1 perldl> p lgamma(79) 8 1 perldl> p lgamma(77) 0 1 perldl> p lgamma(76) 251 1
everything looks okay when these values are called as part of a larger array: perldl> $a=sequence(100)+1 perldl> ($b,)=lgamma($a) perldl> p $b [0 0 0.69314718 1.7917595 3.1780538 4.7874917 6.5792512 8.5251614 10.604603 12.801827 15.104413 17.502308 19.987214 22.552164 25.191221 27.899271 30.67186 33.505073 36.395445 39.339884 42.335616 45.380139 48.471181 51.606676 54.784729 58.003605 61.261702 64.557539 67.889743 71.257039 74.658236 78.092224 81.557959 85.054467 88.580828 92.136176 95.719695 99.330612 102.9682 106.63176 110.32064 114.03421 117.77188 121.53308 125.31727 129.12393 132.95258 136.80272 140.67392 144.56574 148.47777 152.40959 156.36084 160.33113 164.32011 168.32745 172.3528 176.39585 180.45629 184.53383 188.62817 192.73905 196.86618 201.00932 205.1682 209.34259 213.53224 217.73693 221.95644 226.19055 230.43904 234.70172 238.97839 243.26885 247.57291 251.8904 256.22114 260.56494 264.92165 269.2911 273.67312 278.06757 282.47429 286.89313 291.32395 295.7666 300.22095 304.68686 309.16419 313.65283 318.15264 322.6635 327.18529 331.71789 336.26118 340.81506 345.37941 349.95412 354.53909 359.13421] any ideas? Steve
_______________________________________________ Perldl mailing list [email protected] http://mailman.jach.hawaii.edu/mailman/listinfo/perldl
