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

[Raz0r]

Помогите с задачкой, нужно найти корень уравнения методом касательных, впринципе само решение, насколько я понял, в следующем коде

while (x>e)
    x=x-f(x)/f1(x);

где е - точность
f(x)-значение функции в точке х
f1(x)-значение х в производной от функции
х -точка из которой провели касательную
Но у меня почему то программа зацикливается после первого прохода... В чём ошибка?

[_Dmitry]

Замени while (abs(f(x))>e)

[Raz0r]

Не помогло, условие сразу выполняется и он выводит точку, откуда касательная проведена...

[_Dmitry]

Тогда приведи код полностью

[Raz0r]

#include<stdio.h>
#include<conio.h>
#include<math.h>

float f (float xf)
{
 return exp(-xf*xf+1)+xf-3;
}

float f1 (float xf1)
{
 return -2*xf1*exp(-xf1*xf1+1)+1;
}

float f2 (float xf2)
{
 return (4*xf2*xf2-2)*exp(-xf2*xf2+1);
}

void main()
{
float a, b, e, x, y;
clrscr();
printf("Input interval (a,b) and accuracy E:\n");
scanf("%f%f%f", &a, &b, &e);
if (f(a)*f(b)>0)
 {printf("There are no roots.\n");
  getch();
  return;
 }
 else if (f2(a)*f(a)>0) x=a;
  else x=b;
  while (abs(f(x))>e)
    x=x-f(x)/f1(x);
printf("Root of a equation is %2.3f", x);
getch();
}

[_Dmitry]

Что такое f2 я так и не понял. Убери эту функцию исправь
else if (f(a)*f1(a)>0) x=a;

[Raz0r]

Это функция проверяет мою функцию на выпуклость/вогнутость, там мне заданы координаты отрезка, и это условие определяет из какой точки касательную проводить. А вот дальше проблемы начинаются...

[_Dmitry]

Это же самое делается через производную

if (f(a)*f1(a)>0) x=a;
  else x=b;

теперь ведь работает, результат 3

[Raz0r]

А не через вторую производную? f2 - это она и есть. Попробовал как ты написал, результат тот же, а что с while делать, ума не приложу...

[_Dmitry]

Окончательный код вот

#include<stdio.h>
#include<conio.h>
#include<math.h>

float f (float xf)
{
 return exp(-xf*xf+1)+xf-3;
}

float f1 (float xf1)
{
 return -2*xf1*exp(-xf1*xf1+1)+1;
}

void main()
{
float a, b, e, x, y;
clrscr();
printf("Input interval (a,b) and accuracy E:\n");
scanf("%f%f%f", &a, &b, &e);
if (f(a)*f(b)>0)
 {printf("There are no roots.\n");
  getch();
  return;
 }
 else if (f(a)*f1(a)>0) x=a;
  else x=b;
  while (abs(f(x))>e)
    x=x-f(x)/f1(x);
printf("Root of a equation is %2.3f", x);
getch();
}

по условию, решается задача поиска корня уравнения f(x)=0 на отрезке [a,b], цикл

  while (abs(f(x))>e)
    x=x-f(x)/f1(x);

выполняется, пока f(x)>e, когда становится f(x)<=e расчет заканчивается