Лаб 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)