Shevchenko / 1sem / KursvCodeV11

clc
clear
close all
%Y1=cos(x) %Y2=sqrt(x)-4
disp('Корепанов Денис Максимович 3206')
disp('Вариант 11')
format long
n=input('Необходимая точка пересечения:');
eps=input('Погрешность:');
xq=0; %начальное значение
dx=0.001; %изменение аргумента
q=1;
r=1; %переменные для счетчиков
xt(q)=xq;

while r<=n
  dy0=cos(xt(q))-(sqrt(xt(q))-4); %разность значений функций в точке 0
  q=q+1;
xt(q)=xt(q-1)+dx;
dy1=cos(xt(q))-(sqrt(xt(q))-4);
  while dy1*dy0>0 %цикл нахождения точки пересечения
    dy1=cos(xt(q))-(sqrt(xt(q))-4);
    xt(q)=xt(q)+dx;
    xt(q-1)=xt(q-1)+dx;
  end
w=1;
c=(xt(q)+xt(q-1))/2;
    while (xt(q)-xt(q-1))>eps %точное нахождение точки пересечения путем половинного деления отрезка
      c=(xt(q)+xt(q-1))/2;
      w=w+1;
      dyc=cos(c)-(sqrt(c)-4);
      if dyc*dy0>0 %определение приближенного значения x в зависимости от знака
       xt(q)=xt(q);
       xt(q-1)==xt(q)-(c-xt(q-1));
      elseif dyc*dy0<0
       xt(q-1)=xt(q-1);
       xt(q)=xt(q-1)+(c-xt(q-1));
      end
end
q=q+1;
xt(q)=xt(q-1)+dx;
r=r+1;
end
fplot('cos(x)',[0 c+100*dx],'r'); %Вывод графиков функций y1(x) и y2(x)
hold on
fplot('sqrt(x)-4',[0 c+100*dx]);
grid on
legend('cos(x)','sqrt(x)-4');
disp('Координата пересечения x='); disp(c); %вывод значения х
disp('Количество итераций');  disp(w);