Untitled

mail@pastecode.io avatar
unknown
plain_text
a year ago
2.6 kB
2
Indexable
Never
from math import *


print('------------------------------------------------------------------------------------------------------------------------------')
print('                                           Лабораторная работа №1                                                             ')
print('                                                Вариант 4                                                                     ')
print('                                                                                   Выполнил: Бикеев Вадим Владимрович         ')
print('                                                                                   Группа: М6О-308С-21                        ')
print('                                                                                   Преподаватель: Котельников Михаил Вадимович')
print('------------------------------------------------------------------------------------------------------------------------------') 

def f (x): #Исходная функция
    return x**3 + x**2 - x - 0.5

def fI (x): #Первая производная
    return 3*x**2 + 2*x - 1


e = 0.001 #Точность 
x_i = 0.5 #Начальное приближение
print('Начальное приближение :', 'x =' , x_i)
print('------------------------------------------------------------------------------------------------------------------------------') 
print()


x_j = round(x_i - f(x_i)/fI(x_i), 4)#Находим х после первой итерации
print('Номер итерации 1 :', 'x =' , x_j)
print('------------------------------------------------------------------------------------------------------------------------------') 
print()

i = 1
while abs(x_j - x_i) >= e: #Делаем новые итерации, пока разность между элементами не станет меньше точности 
    x = x_i
    x_i = x_j
    x_j = round(x_i - f(x_i)/fI(x_i), 4)
    i += 1
    print('Номер итерации ', i, ":", 'x =' , x_j)
    print('Разность между последними итерациями равна', abs(x_j - x_i))
    print('------------------------------------------------------------------------------------------------------------------------------') 
    print()

print('Положительным решением нелинейного уравнения является x =' , x_j)