Лаб 2 / z2

syms x

# Функция №7
function Fx = v1 (x)
  Fx = sin(sin(x.^2)) - 0.5
endfunction

# Производная функции №7
function dFx = v2 (x)
  dFx = 2.*x.*cos(x.^2).*cos(sin(x.^2))
endfunction

x0 = 1.5
# Метод Ньютона
for i = 1:15
  x0 = x0 - v1(x0) / v2(x0)
 end

 % Итоговое значение, невязка
 zeroX = x0
 zeroN = v1(x0)

 # Погрешность
 diff = abs(fzero('@(x) sin(sin(x.^2)) - 0.5', [1, 2]) - zeroX)