Решение нелинейного уравнения методом Ньютона (касательных)

@seco234 · Регистрация: 19.11.2016

Здравствуйте!
Прошу вашей помощи в решении нелинейного уравнения методом Ньютона.
Уравнение: exp^(-cos(x))-x^3=0

Matlab M

x = -2 : 0.01 : 2;
y = exp(-cos(x)) - x.^3;
plot(x, y)
grid;
eps = 0.0001;
y = @(x) exp(-cos(x)) - x.^3;
a = 0;
b = 1;
 
syms x;
diff1 = diff(exp(-cos(x)) - x.^3)
diff2 = diff(exp(-cos(x)) - x.^3,2)
 
if (y(a) * y (b) > 0)
    error('Значения функции на концах интервала имеют не разные знаки.')
end;
 
if (y(a) * diff2(a) > 0)
    x0 = a;
else
    x0 = b;
end;
 
i=0;
xk = x0;
xn = b;
while abs(xk-x)> eps
    xn = xk;
    xk = xn - y(xn) / diff1(xn);
    i = i + 1;
end
xk
i;

У меня выходит ошибка:

diff1 =

exp(-cos(x))*sin(x) - 3*x^2

diff2 =

exp(-cos(x))*sin(x)^2 - 6*x + exp(-cos(x))*cos(x)

Subscript indices must either be real positive integers or logicals.

Error in sym/subsref (line 841)
R_tilde = builtin('subsref',L_tilde,Idx);

Error in Program3 (line 18)
if (y(a) * diff2(a) > 0)

Я в матлабе полнейший ноль и не пойму, что он от меня хочет.
Может, проблема в том, что x символическая, а всё остальное - нет? Без syms x не находятся производные. Короче, помогите, пожалуйста!

@Centurio · 27.10.2017, 21:20

seco234, посмотрите внизу в разделе "Похожие темы"

@Centurio Модератор 1172 / 1068 / 317 Регистрация: 13.09.2015 Сообщений: 3,775
	27.10.2017, 21:20	2
	seco234, посмотрите внизу в разделе "Похожие темы" 0

Опции темы