Octave / BisHord

f = @(x)(x.^2 - log(x));
f1 = @(x)(x.^2 + log(x)-4);
m = [-0.5:0.001:0.5];
m2 = [-0.5:0.001:0.5];
subplot(1,2,1)
plot(m,f(m))
xlabel('x')
ylabel('y')
title('Variant 13')
grid on
subplot(1,2,2)
plot(m2,f1(m2))
xlabel('x')
ylabel('y')
title('Variant 17')
grid on
disp ('Метод бисекции')
z = bisection(f,0,1, 0.001)
z2 = bisection(f1,1,2,0.001)
disp('Метод хорд')
z = hord(f,1,2,0.002)
z2 = hord(f1,1,2,0.001)
disp('Метод fzero')
%[X,FVAL,INFO,OUTPUT] = fzero(f,[13;15])
[X,FVAL,INFO,OUTPUT] = fzero(f1,[-0.5;0.5])
disp('Метод fsolve')
[X,FVAL,INFO,OUTPUT] = fsolve(f,[13;15])
[X,FVAL,INFO,OUTPUT] = fsolve(f1,[-0.5;0.5])