Scribd Logo
Upload

Welcome to Scribd!

  • 49 views

    Método de Ferrari, Ejemplo 4

    Uploaded by

    Sergio Sanchez
    This document presents Ferrari's method for solving quartic equations of the form ax^4 + bx^3 + cx^2 + dx + e = 0. It describes shifting the equation and using Cardano's formula to find the roots of a related cubic equation, then back-substituting to find the four roots of the original quartic equation. An example application of the Ferrari's method to solve the equation 2x^4 + x^3 + 2x^2 - 3x - 2 = 0 is shown, returning the four complex roots.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd

    Método de Ferrari, Ejemplo 4

    Uploaded by

    Sergio Sanchez
    49 views2 pages
    This document presents Ferrari's method for solving quartic equations of the form ax^4 + bx^3 + cx^2 + dx + e = 0. It describes shifting the equation and using Cardano's formula to find the roots of a related cubic equation, then back-substituting to find the four roots of the original quartic equation. An example application of the Ferrari's method to solve the equation 2x^4 + x^3 + 2x^2 - 3x - 2 = 0 is shown, returning the four complex roots.

    Original Title

    Método de ferrari, Ejemplo 4

    Copyright

    © © All Rights Reserved

    Available Formats

    PDF, TXT or read online from Scribd

    Did you find this document useful?

    Is this content inappropriate?

    This document presents Ferrari's method for solving quartic equations of the form ax^4 + bx^3 + cx^2 + dx + e = 0. It describes shifting the equation and using Cardano's formula to find the roots of a related cubic equation, then back-substituting to find the four roots of the original quartic equation. An example application of the Ferrari's method to solve the equation 2x^4 + x^3 + 2x^2 - 3x - 2 = 0 is shown, returning the four complex roots.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    49 views2 pages

    Método de Ferrari, Ejemplo 4

    Uploaded by

    Sergio Sanchez
    This document presents Ferrari's method for solving quartic equations of the form ax^4 + bx^3 + cx^2 + dx + e = 0. It describes shifting the equation and using Cardano's formula to find the roots of a related cubic equation, then back-substituting to find the four roots of the original quartic equation. An example application of the Ferrari's method to solve the equation 2x^4 + x^3 + 2x^2 - 3x - 2 = 0 is shown, returning the four complex roots.

    Copyright:

    © All Rights Reserved
    Download as PDF, TXT or read online from Scribd
    Download as pdf or txt
    of 2

    M?

    todo de ferrari, Ejemplo 4

    November 30, 2019

    [1]: import numpy as np


    import cmath

    J=cmath.exp(2j*cmath.pi/3)
    Jc=1/J

    def ferrari(a,b,c,d,e):
    "Ferrarai's Method"
    "resolution of P=ax^4+bx^3+cx^2+dx+e=0, coeffs reals"
    "First shift : x= z-b/4/a => P=z^4+pz^2+qz+r"
    z0=b/4/a
    a2,b2,c2,d2 = a*a,b*b,c*c,d*d
    p = -3*b2/(8*a2)+c/a
    q = b*b2/8/a/a2 - 1/2*b*c/a2 + d/a
    r = -3/256*b2*b2/a2/a2 +c*b2/a2/a/16-b*d/a2/4+e/a
    "Second find y so P2=Ay^3+By^2+Cy+D=0"
    A=8
    B=-4*p
    C=-8*r
    D=4*r*p-q*q
    y0,y1,y2=Cardano(A,B,C,D)
    if abs(y1.imag)<abs(y0.imag): y0=y1
    if abs(y2.imag)<abs(y0.imag): y0=y2
    a0=(-p+2*y0)**.5
    if a0==0 : b0=y0**2-r
    else : b0=-q/2/a0
    r0,r1=Roots_2(1,a0,y0+b0)
    r2,r3=Roots_2(1,-a0,y0-b0)
    return (r0-z0,r1-z0,r2-z0,r3-z0)

    #~ @jit(nopython=True)
    def Cardano(a,b,c,d):
    z0=b/3/a
    a2,b2 = a*a,b*b
    p=-b2/3/a2 +c/a
    q=(b/27*(2*b2/a2-9*c/a)+d)/a

    1
    D=-4*p*p*p-27*q*q
    r=cmath.sqrt(-D/27+0j)
    u=((-q-r)/2)**0.33333333333333333333333
    v=((-q+r)/2)**0.33333333333333333333333
    w=u*v
    w0=abs(w+p/3)
    w1=abs(w*J+p/3)
    w2=abs(w*Jc+p/3)
    if w0<w1:
    if w2<w0 : v*=Jc
    elif w2<w1 : v*=Jc
    else: v*=J
    return u+v-z0, u*J+v*Jc-z0, u*Jc+v*J-z0

    #~ @jit(nopython=True)
    def Roots_2(a,b,c):
    bp=b/2
    delta=bp*bp-a*c
    r1=(-bp-delta**.5)/a
    r2=-r1-b/a
    return r1,r2
    [2]: ferrari(2,1,2,-3,-2)
    [2]: ((-0.7279799860324166+1.1258550469212985j),
    (-0.7249079646793538-1.3282742563830554j),
    (0.19574974703672704+0.528218846563451j),
    (0.7571382036750434-0.325799637101694j))
    [ ]:

    You might also like