Quartic function
The quartic function is the bridge between the cubic function and more advanced functions such as the quintic and sextic.
Objective
edit
|
Lesson
edit
IntroductioneditThe quartic function is the sum of powers of from through :
usually written as:
If the function becomes Within this page we'll say that:
The quartic equation is the quartic function equated to zero: . Roots of the function are values of that satisfy the quartic equation.
Function as product of linear function and cubicedit
Function defined by 5 pointsedit
Function defined by 3 points and 2 slopesedit
Examplesedit
|
First and second derivatives
edit
Points of inflectionedit
Maxima and minimaedit
Quartic with 2 stationary pointsedit
|
The simplest quartic function
edit
The simplest quartic function has coefficients Red line in diagram has equation: First derivative (not shown) of When There is a stationary point on when point Second derivative (not shown) of When There is a point of inflection on when For every non-zero value of is positive. To left and right of point is always concave up. Point is both local minimum and absolute minimum.
Solve:
This is equivalent to finding a root of function If you use Newton's method to find a root of this would be more efficient than solving |
Roots of equal absolute value
edit
The standard quartic function: For in substitute Call this For in substitute Call this Combine and to eliminate and produce an equation in
If is a solution and function has 2 roots of form where
The method works with complex roots of equal absolute value:
|
Equal roots
edit
Equal roots occur when the function and the slope of the function both equal zero.
where:
where:
where:
where:
where:
where:
From
+ 2048aaaaacddeeee - 768aaaaaddddeee - 1536aaaabcdddeee + 576aaaabdddddee
- 1024aaaacccddeee + 1536aaaaccddddee - 648aaaacdddddde + 81aaaadddddddd
+ 1152aaabbccddeee - 480aaabbcddddee + 18aaabbdddddde - 640aaabcccdddee
+ 384aaabccddddde - 54aaabcddddddd + 128aaacccccddee - 80aaaccccdddde
+ 12aaacccdddddd - 216aabbbbcddeee + 81aabbbbddddee + 144aabbbccdddee
- 86aabbbcddddde + 12aabbbddddddd - 32aabbccccddee + 20aabbcccdddde
- 3aabbccdddddd which, by removing values (common to all values), may be reduced to: status = (
+ 2048aaaceeee - 768aaaddeee - 1536aabcdeee + 576aabdddee
- 1024aaccceee + 1536aaccddee - 648aacdddde + 81aadddddd
+ 1152abbcceee - 480abbcddee + 18abbdddde - 640abcccdee
+ 384abccddde - 54abcddddd + 128acccccee - 80accccdde
+ 12acccdddd - 216bbbbceee + 81bbbbddee + 144bbbccdee
- 86bbbcddde + 12bbbddddd - 32bbccccee + 20bbcccdde
- 3bbccdddd
) If there are at least 2 equal roots which may be calculated as shown below.
No equal rootsedit
Exactly 2 equal rootsedit
Exactly 3 equal rootsedit
Four equal rootsedit
Two pairs of equal rootsedit
Summmaryedit
Cautionedit
|
Depressed quartic
edit
A depressed quartic is any quartic function with any one or more of coefficients missing. Within this section a depressed quartic has coefficient missing. To produce the depressed quartic:
Let Substitute in expand and simplify:
where:
Be prepared for the possibility that any 1 or more of may be zero. Coefficient B missingedit
Coefficient C missingedit
Resolvent cubicedit
|
Solving quartic equation
editThis section presents 4 examples that show how to use the depressed quartic and the resolvent cubic to solve the quartic equation.
Four real roots
edit
Calculate roots of: |
Calculate coefficients of depressed quartic:
a,b,c,d,e = 1,-1,-19,-11,30
A = 16*a*c - 6*b*b
B = 64*a*a*d - 32*a*b*c + 8*b*b*b
C = 256*a*a*a*e - 64*a*a*b*d + 16*a*b*b*c - 3*b*b*b*b
print (A,B,C)
-310 -1320 6669
|
Calculate coefficients of resolvent cubic:
P = 64
Q = 32*A
R = 4*A*A - 16*C
S = -B*B
print (P,Q,R,S)
64 -9920 277696 -1742400
|
Calculate roots of cubic function:
There are 3 real, positive roots:
Using 3 roots of calculate 4 roots of
# python code
for U in (9, 25, 121) :
print ('\nU =', U)
sqrtU = U ** 0.5
for u in (sqrtU, -sqrtU) :
V = -(A/2 + U) - B/(4*u)
v = V ** .5
for t in (u+v, u-v) :
x = (-b+t) / (4*a)
y = a*x**4 + b*x**3 + c*x**2 + d*x + e
print ('x:',x, '; y:',y)
U = 9
x: 5.0 ; y: 0.0
x: -3.0 ; y: 0.0
x: 1.0 ; y: 0.0
x: -2.0 ; y: 0.0
U = 25
x: 5.0 ; y: 0.0
x: -2.0 ; y: 0.0
x: 1.0 ; y: 0.0
x: -3.0 ; y: 0.0
U = 121
x: 5.0 ; y: 0.0
x: 1.0 ; y: 0.0
x: -2.0 ; y: 0.0
x: -3.0 ; y: 0.0 Roots of are: |
All 3 values of produce the same results, but not in same sequence.
It is not necessary to calculate all 3 roots of resolvent cubic. Any one non-zero root is sufficient to do the job.
Two real and two complex roots
editExample 1
edit
Calculate roots of: |
Calculate coefficients of depressed quartic:
Calculate coefficients of resolvent cubic:
Calculate one real root of cubic function: is one real root. Choose Calculate roots of
|
Example 2
edit
Calculate roots of: |
Calculate coefficients of depressed quartic:
a,b,c,d,e = 3, -6, -41, 44, -189
A = 16*a*c - 6*b*b
B = 64*a*a*d - 32*a*b*c + 8*b*b*b
C = 256*a*a*a*e - 64*a*a*b*d + 16*a*b*b*c - 3*b*b*b*b
print (A,B,C)
-2184, 0, -1229040 Notice that coefficient |
Calculate coefficients of resolvent cubic:
P = 64
Q = 32*A
R = 4*A*A - 16*C
S = -B*B
print (P,Q,R,S)
1, -1092, 605376, 0 Notice that coefficient |
Calculate roots of cubic function:
Roots are
Value cannot be used because it will cause error Divide by zero
at statement V = -(A/2 + U) - B/(4*u)
.
Calculate roots of
# python code
U = 546+554.3103823671355j
print ('\nU =',U)
sqrtU = U ** 0.5
for u in (sqrtU, -sqrtU) :
V = -(A/2 + U) - B/(4*u)
v = V ** 0.5
s1 = '\nu,v'
print (s1,eval(s1))
for t in (u+v, u-v) :
x = (-b+t)/(4*a)
# Check result. Expecting sum = 0.
sum = a*x**4 + b*x**3 + c*x**2 + d*x + e
print ('x =', x,'; sum =',sum)
U = (546+554.3103823671355j)
u,v ((25.729935131257832+10.771701901683684j), (25.729935131257832-10.771701901683684j))
x = (4.788322521876306 + 0j) ; sum = (1.9895196601282805e-13 + 0j)
x = (0.5 + 1.795283650280614j) ; sum = (5.684341886080802e-14 + 0j)
u,v ((-25.729935131257832-10.771701901683684j), (25.729935131257832-10.771701901683684j))
x = (0.5 - 1.795283650280614j) ; sum = (5.684341886080802e-14 + 0j)
x = (-3.7883225218763052 + 0j) ; sum = (1.7053025658242404e-13 + 0j) Values of are: |
Depressed quartic as quadratic
editIn this example coefficient of depressed quartic
Therefore, resolvent cubic can be ignored and depressed quartic processed as quadratic in
where
Solutions of this quadratic are:
T1,T2 = 2648.1182474349434, -464.11824743494344
t1 = T1 ** 0.5; t2 = ((-T2) ** 0.5) * 1j
for t in (t1,-t1,t2,-t2) :
x = (-b+t)/(4*a)
# Check result. Expecting sum = 0.
sum = a*x**4 + b*x**3 + c*x**2 + d*x + e
print ('x =', x,'; sum =',sum)
x = 4.788322521876305 ; sum = -1.7053025658242404e-13
x = -3.788322521876305 ; sum = -1.7053025658242404e-13
x = (0.5 + 1.7952836502806138j) ; sum = (-2.842170943040401e-14 + 0j)
x = (0.5 - 1.7952836502806138j) ; sum = (-2.842170943040401e-14 + 0j)
or
x = 0.5 ± 4.288322521876305, 0.5 ± 1.7952836502806138j
With precision of 15, values of are same as those shown above.
When roots of quartic function are of form p ± q, p ± r,
coefficient of depressed function
Four complex roots
edit
Calculate roots of: |
Calculate coefficients of depressed quartic:
4128 344064 9683200
|
Calculate coefficients of resolvent cubic:
64 132096 -86769664 -118380036096
|
Calculate one root of cubic function:
There are 3 real roots: Choose
Negative is chosen here to show that any 1 of the roots produces the correct result.
Calculate roots of
# python code
U = -784
u1 = 1j * (-U)**.5
for u in (u1, -u1) :
V = -(A/2 + U) - B/(4*u)
v = V**.5
for t in (u+v, u-v) :
x = (-b+t)/(4*a)
# Check result. Expecting sum = 0.
sum = a*x**4 + b*x**3 + c*x**2 + d*x + e
print ('x =', x,'; sum =',sum)
# python expresses complex numbers with 'j'.
x = (13+19j) ; sum = 0j
x = (-3-5j) ; sum = 0j
x = (13-19j) ; sum = 0j
x = (-3+5j) ; sum = 0j
|
Quartic formula
editThe substitutions made above can be used to produce a formula for the solution of the quartic equation.
See main articles "The general case" or "General formula for roots."
Both links above point to formula for equation |
Given quartic equation: calculate the 4 values of
where:
Coefficients of depressed quartic:
|
Coefficients of resolvent cubic:
|
Coefficients of depressed cubic:
|
One root of resolvent cubic:
may be negative.
|
One root of quartic: may be positive or negative.
may be positive or negative.
|
Formula above produces one value of Python code below utilizes and to produce 4 values of and then, four values of
An example:
edit
Calculate roots of # Python code.
a,b,c,d,e = 4, 4, -75, -776, -1869
values_of_t = [
t
# Coefficients of depressed quartic:
for A in (16*a*c - 6*b*b,)
for B in (64*a*a*d -32*a*b*c + 8*b*b*b,)
for C in (256*a*a*a*e - 64*a*a*b*d + 16*a*b*b*c - 3*b*b*b*b,)
# Coefficients of resolvent cubic:
for a1 in (64,)
for b1 in (32*A,)
for c1 in (4*A*A - 16*C,)
for d1 in (-B*B,)
for U in [
# The resolvent cubic:
(-b1+t1)/(3*a1)
# Coefficients of depressed resolvent cubic:
for A1 in (9*a1*c1 - 3*b1*b1,)
for B1 in (27*a1*a1*d1 - 9*a1*b1*c1 + 2*b1*b1*b1,)
# One root of resolvent cubic:
for C1 in (-A1/3,)
for Δ in (B1*B1 - 4*C1*C1*C1,)
for δ in (Δ**0.5,)
for W in ((-B1 + δ)/2,)
for w in (W**(1/3),)
for t1 in (w + C1/w,) # See note below.
]
# Prepare to calculate 4 values of t.
for u1 in (U**.5,)
for v1 in ( -(A/2 + U) ,)
# Calculate 4 values of t.
for u in (u1, -u1,)
for V in ( v1 - B/(4*u),)
for v in (V**.5,)
for t in (u+v, u-v)
]
print ('values_of_t =', values_of_t)
# Python code.
# Calculate 4 separate roots.
values_of_x = [
(-b + t)/(4*a)
for t in values_of_t
]
print ('values_of_x =', values_of_x)
values_of_x = [7, -3, (-2.5 + 4j), (-2.5 - 4j)]
In python the imaginary part of a complex number is shown with instead of If contains 4 equal roots and If contains 3 or more equal roots, statement Before using this formula, check for equal roots as in "Exactly 3 equal roots" above.
|
In practice
edit
The following Python code implements the quartic formula. However, under statement
Examplesedit
|
Two Conic Sections
editExamples of conic sections include: ellipse, circle, parabola and hyperbola.
This section presents examples of two conic sections, circle and ellipse, and how to calculate the coordinates of the point/s of intersection, if any, of the two sections.
Let one section with name have equation
Let other section with name have equation
Because there can be as many as 4 points of intersection, a special "resolvent" quartic function is used to calculate the coordinates of the point/s of intersection.
Coefficients of associated "resolvent" quartic are calculated as follows:
# python code
def intersection_of_2_conic_sections (abcdef, ABCDEF) :
'''
A_,B_,C_,D_,E_ = intersection_of_2_conic_sections (abcdef, ABCDEF)
where A_,B_,C_,D_,E_ are coefficients of associated resolvent quartic function:
y = f(x) = A_*x**4 + B_*x**3 + C_*x**2 + D_*x + E_
'''
A,B,C,D,E,F = ABCDEF
a,b,c,d,e,f = abcdef
G = ((-1)*(B)*(a) + (1)*(A)*(b))
H = ((-1)*(B)*(d) + (1)*(D)*(b))
I = ((-1)*(B)*(f) + (1)*(F)*(b))
J = ((-1)*(C)*(a) + (1)*(A)*(c))
K = ((-1)*(C)*(d) + (-1)*(E)*(a) + (1)*(A)*(e) + (1)*(D)*(c))
L = ((-1)*(C)*(f) + (-1)*(E)*(d) + (1)*(D)*(e) + (1)*(F)*(c))
M = ((-1)*(E)*(f) + (1)*(F)*(e))
g = ((-1)*(C)*(b) + (1)*(B)*(c))
h = ((-1)*(E)*(b) + (1)*(B)*(e))
i = ((-1)*(A)*(b) + (1)*(B)*(a))
j = ((-1)*(D)*(b) + (1)*(B)*(d))
k = ((-1)*(F)*(b) + (1)*(B)*(f))
A_ = ((-1)*(J)*(g) + (1)*(G)*(i))
B_ = ((-1)*(J)*(h) + (-1)*(K)*(g) + (1)*(G)*(j) + (1)*(H)*(i))
C_ = ((-1)*(K)*(h) + (-1)*(L)*(g) + (1)*(G)*(k) + (1)*(H)*(j) + (1)*(I)*(i))
D_ = ((-1)*(L)*(h) + (-1)*(M)*(g) + (1)*(H)*(k) + (1)*(I)*(j))
E_ = ((-1)*(M)*(h) + (1)*(I)*(k))
str1 = 'y = ({})x^4 + ({})x^3 + ({})x^2 + ({})x + ({}) '.format(A_,B_,C_,D_,E_)
print (str1)
return A_,B_,C_,D_,E_
With no common point
edit
Let ellipse (red curve) have equation: Let circle (blue curve) have equation: Then, resolvent quartic function (black curve) has equation:
has no real roots. Therefore, there is no point of intersection. |
With one common point
edit
Let ellipse (red curve) have equation: Let circle (blue curve) have equation: Then, resolvent quartic function (black curve) has equation:
Roots of are: has 2 equal, real roots at effectively 1 real root where Therefore, there is one point of intersection where |
With two common points
editExample 1
edit
Let ellipse (red curve) have equation: Let circle (blue curve) have equation: Then, resolvent quartic function (black curve) has equation:
Roots of are: has 2 unique, real roots at Therefore, there are two points of intersections where |
Example 2
edit
Let ellipse (red curve) have equation: Let circle (blue curve) have equation: Then, resolvent quartic function (black curve) has equation:
Roots of are: has 2 pairs of equal roots at effectively 2 real roots. Therefore, there are two points of intersection where |
With 3 common points
edit
Let ellipse (red curve) have equation: Let circle (blue curve) have equation: Then, resolvent quartic function (black curve) has equation:
Roots of are: has 1 pair of equal roots at and 2 unique, real roots at effectively 3 real roots. Therefore, there are three points of intersection where |
With 4 common points
edit
Let ellipse (red curve) have equation: Let circle (blue curve) have equation: Then, resolvent quartic function (black curve) has equation:
Roots of are: has 4 real roots as shown above. Therefore, there are four points of intersection where |