Fortran, решение нелинейного

[craSu]

Надо решить нелинейное уравнение 0,5^x-(x+2)^2=3
Пробовала методом Ньютона (касательных)
но не выходит, спасибо заранее



module mod
implicit none
contains
function prsv(x) result(prF)
integer, intent(in) :: x !input
integer :: prF !output
prF=(0.5**x)*log(0.5)-2x-4
return
end function prsv

function func(x) result(F)
integer, intent(in) :: x !input
integer :: F ! output
F=(0.5**x)-(x+2)**2-3
return
end function func

end module mod

program nelin
use mod
integer :: x
! integer :: func, prsv
integer :: x1, x2, eps, n

eps=0.0001
x1=0.5
n=0
if ((abs(x2-x1)<eps).or.(abs(func(x1))<eps)) then
x2=x1-func(x1)/prsv(x1)
n=n+1
else
print*,x2,n
end if

[type_Oleg]

Метод касательных не всегда сходится.
Попробуйте метод бисекций ( половинного деления) , он более всеядный.
Или нужно именно секущих ?

[craSu]

Допускается использование метод касательных, метод бисекции, метод хорд, метод простых итераций (если что-то не забыла).
Какой тогда может будет проще в реализации?

[craSu]

да, допускается.
тогда попробую

[type_Oleg]

Да, лучше бисекций.
Посмотрите про https://Метод Ньютона - для него условия применения аж целый академик и нобелевский лауреат Канторович вычисил.

А для бисекций каких-то особых требований к функции нет.

[Vago]

Простите, а почему x и функции - типа integer ?!