Program do liczenia pierwiastków równania sześciennego o współczynnikach rzeczywistych Poniżej podano przykład kodu w języku Python, znajdujący pierwiastki równania sześciennego wg wzorów wyżej podanych dla równania sześciennego o współczynnikach rzeczywistych. W programie obliczane są pierwiastki dla 5 zestawów danych (tj. dane są współczynniki a, b, c, d równania sześciennego) i rysowane są wykresy funkcji sześciennych. Użytkownik może dodać własne zestawy współczynników (linia 55 kodu). W linii 77 kodu: x = linespace(-1, 4, 1000) liczby -1, 4 określają zakres zmiennej x, w jakim rysowany jest wykres - można dostosować do własnych przykładów.# JK 2024.07.19 v.2 Program oblicza pierwiastki równania sześciennego # o współczynnikach rzeczywistych wg wzorów podanych w: # W. Zakowski, Mały poradnik matematyczny, Warszawa, 1977, str. 255-256.
import numpy as np, matplotlib.pyplot as plt, cmath
# GłóWNA FUNKCJA, OBLICZAJĄCA PIERWIASTKI RÓWNANIA SZEŚCIENNEGO #
def solve_cubic(a, b, c, d): p = c/a - b**2 / (3*a**2) q = 2 * b**3 / (27*a**3) + d/a - b * c / (3*a**2) D = (q / 2)**2 + (p / 3)**3 s= b / (3*a)