![]() ![]() int gsl_poly_complex_solve_quadratic ( double a, double b, double c, gsl_complex * z0, gsl_complex * z1 ) ¶ Polynomials with small integer coefficients the discriminant can alwaysīe computed exactly. May cause a discrete change in the number of roots. This will be subject to rounding and cancellationĮrrors when computed in double precision, and will also be subject toĮrrors if the coefficients of the polynomial are inexact. The number of roots found depends on the sign of the discriminant For example will have two roots, which happen ![]() The case of coincident roots is not considered Real roots are found they are stored in x0 and x1 inĪscending order. ![]() If no real rootsĪre found then x0 and x1 are not modified. The number of real roots (either zero, one or two) is returned, and Using the notation above,Ĭan be evaluated by calling gsl_poly_dd_eval() and using On output theĭivided-differences of the Hermite representation are stored in the arrayĭd, also of length 2* size. Incorported into the usual divided-difference algorithm by forming a new Provided in the array dya also of length size. Hermite interpolationĬonstructs polynomials which also match first derivatives which are Interpolating Hermite polynomial for the points stored in ![]() int gsl_poly_dd_hermite_init ( double dd, double za, const double xa, const double ya, const double dya, const size_t size ) ¶ A workspace of length size must be provided in theĪrray w. Point xp are stored in the array c also of length On output the Taylor coefficients of the polynomial expanded about the Is supplied in the arrays dd and xa of length size. This function converts the divided-difference representation of a int gsl_poly_dd_taylor ( double c, double xp, const double dd, const double xa, size_t size, double w ) ¶ An inline version of this function is used when HAVE_INLINE is defined. In the arrays dd and xa of length size at the point This function evaluates the polynomial stored in divided-difference form double gsl_poly_dd_eval ( const double dd, const double xa, const size_t size, const double x ) ¶ On output theĭivided-differences of ( xa, ya) are stored in the arrayĭd, also of length size. Interpolating polynomial for the points stored in This function computes a divided-difference representation of the int gsl_poly_dd_init ( double dd, const double xa, const double ya, size_t size ) ¶ Are discussed in Burden and Faires, section 3.4. ![]()
0 Comments
Leave a Reply. |
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |