Решение уравнения методом Ньютона. Локализация корня

[arktest]

Собственно, программа имеется, все довольно таки просто, но необходимо локализовать корень и найти именно положительный корень уравнения. Я так понял, нужно найти начальное приближение Xn и Xm. Приближение необходимо искать методом деления пополам. Как мне объяснили, необходимо выставить пределы для деления пополам от 0 до 100 и потом, когда метод будет сходится уже что-то принимать за Xn и Xm. Я в этом не силен. Необходима помощь, желательно программно.

#include <cstdlib>
#include <iostream>
#include <math.h>
#include <iomanip>
 
using namespace std;
double e = 0.001; 
 
 
double funct(double x)
{
    double f = x*x+4*sin(5*x)+exp(x+1)-3.2;
    return f;
}
 
 
double diff(double x)
{
    double f = 2*x+exp(x+1)+20*cos(5*x);
    return f;
}
 
 
int task_2()
{
    int n = 0;
    double xm = 0.95;
    double xn = 0;
    double xnprev;
    double acc = fabs(xn - xm);
    //int a=0,b=2;
    /*
    double c;
    while ((xn-xm)>=e)
          {
               c=(xn+xm)/2;
               if (funct(xm)*funct(c)<0)
               xn=c;
               else xm=c;                         
           }*/
    cout << "Решение методом Ньютона:" << endl << endl;
    while (acc > e)
    {
        n++;
        xnprev=xn;
        xn = xm - funct(xm)/diff(xm);
        if(xn<0){
            //xn=xnprev;
            //break;
        }
        acc = fabs(xn - xm);
        cout << setw(2) << n << " " << xm << " " << xn << " " << acc << endl;
        xm = xn;
    }
    cout << endl << "Корень уравнения x = " << xn << endl;
    return 0;
}
 
int main()
{
    
          
    setlocale(LC_ALL,"Rus");
    int ch;
   
    cout.precision(4);
    cout.setf(ios::fixed);
    
    task_2();
        
    system("PAUSE");
    return EXIT_SUCCESS;
}

________
Код нужно оформлять по правилам:
тегом [CODE]..[/СODE]
(это кнопочка на панели форматирования с решёточкой #)
Не забывайте об этом!

Модератор.