netlib rpoly.f frustration

zwang5_at_hotmail.com
Date: 03/26/05 

Date: 26 Mar 2005 10:48:36 -0800

Dear all,

I am having trouble to use rpoly subroutine found from
http://www.netlib.org/port/prop.upd/rpoly.f

Here is part of the example and I don't know why I got the error
message "leading coeeficient is zero"?

Thanks!

*******
       double precision alpha1(3)
       print *,'alpha1= ',alpha1
       print *,'alpha1(1)=',alpha1(1)
       print *,'alpha1(2)=',alpha1(2)
       print *,'alpha1(3)=',alpha1(3)
       call rpoly(2,alpha1,rootr,rooti)
*******
Output

 alpha1= 1. 0.01 0.03
 alpha1(1)= 1.
 alpha1(2)= 0.01
 alpha1(3)= 0.03
 ERROR 2 IN
    RPOLY - LEADING COEFFICIENT IS ZERO

*****
      SUBROUTINE RPOLY(DEGREE,COEFF,ZEROR,ZEROI)
C
C FINDS THE ZEROS OF A POLYNOMIAL WITH REAL COEFFICIENTS.
C
C INPUTS -
C
C COEFF - REAL VECTOR OF THE COEFFICIENTS
C IN ORDER OF DECREASING POWERS.
C
C DEGREE - INTEGER DEGREE OF POLYNOMIAL.
C
C OUTPUTS -
C
C ZEROR, ZEROI - REAL VECTORS OF THE REAL
C AND IMAGINARY PARTS OF THE ZEROS.
C
C
C ERROR STATES -
C
C 1 - DEGREE IS LESS THAN 1
C 2 - LEADING COEFFICIENT IS ZERO

Relevant Pages

  • Re: interaction term in linear regression with a dummy coded predictor
    ... of the response vs. X1 or response vs. X2 is not a reasonable thing to ... model coefficient is very large relative to the value of the ... are "zero" with a large interaction. ... Their analysis had been presented to the plant manager and his direct ...
    (sci.stat.math)
  • Re: interaction term in linear regression with a dummy coded predictor
    ... of the response vs. X1 or response vs. X2 is not a reasonable thing to ... model coefficient is very large relative to the value of the ... are "zero" with a large interaction. ... Their analysis had been presented to the plant manager and his direct ...
    (sci.stat.math)
  • Re: How to set up coefficients constraints in my regression?
    ... statistical inference or reasoning behind such that the weighted ... coefficient is zero. ... That's probably why you want to constrain them. ... This is the coefficient vector that solves my equations, which are ...
    (comp.soft-sys.matlab)
  • Regression and p-values
    ... Remember the hypothesis test is whether a coefficient is significantly different from zero as you have set it up. ... If you had a small p-value, then you would reject the null hypothesis and conclude that the slope is different from zero, which is not want you want to do in this case. ...
    (comp.soft-sys.matlab)
  • Re: Multivariable regression (t-stat)
    ... coefficient is actually zero. ... coefficient is considered "statistically significant". ... statistics before assuming that you know much. ...
    (microsoft.public.excel.misc)