Метод касательных Ньютона

clc;
clear;
e = 0.000001; % Погрешность
b = -1;
Xkm = b;
Xkp = 0;
Yo = f(b);
k = 1;
while(abs(Yo) > e)
    Xk = Xkm - f(Xkm)/fp(Xkm);
    Xkm = Xk;
    Xkp = fp(Xk);
    Yo = f(Xkp);
    k = k + 1;
    fprintf('Xk = %f\n',Xk)
end

function[F]=f(t)
    EXP = 2.71;
    F = 7*(EXP ^ -(2*t))*sin(7*t+0.3)-1; 
end

function[F]=fp(t)
    EXP = 2.71;
    F = 49*cos(7*t+0.3)*(EXP ^ (-2*t)) - 14*sin(7*t+0.3)*(EXP ^ (-2*t));
end

clc;
clear;
e = 0.001; % Погрешность
b = -1;
Xkm = b;
Xkp = 0;
Yo = f(b);
k = 1;
while (abs(Yo) > e)
    Xk = Xkm - f(Xkm)/fp(Xkm);
    Xkm = Xk;
    Xkp = fp(Xk);
    Yo = f(Xk);
    k = k + 1;
end
fprintf('Xk = %f\n',Xk)
ezplot(@f,[-1 0]), hold on
plot(Xk,0,'or'), grid on

function [F]=f(t)
    F = 7*exp(-2*t)*sin(7*t+0.3)-1; 
end

function [F]=fp(t)
    F = 49*cos(7*t + 0.3).*exp(-2*t) - 14*sin(7*t + 0.3).*exp(-2*t);
end

fun=@(t)7*exp(-2*t).*sin(7*t+0.3)-1
t=-5:0.1:5;
plot(t,fun(t))
grid on

fzero(fun,-2)