Где тут в коде происходит нахождение корня? (метод половинного деления)

#include <stdio.h>
#include <windows.h>
#include <locale.h>
#include <math.h>
 
HMODULE hLib;
 
float *pa = new float;
float *pb = new float;
float *pc = new float;
float *pe = new float;
float *pgran1 = new float;
float *pgran2 = new float;
int *pz = new int;
 
void vvabc();
void vvz();
void vvgran();
void vve();
int vvod();
float af;
 
int main()
{
    setlocale(LC_ALL,"Russian");//Отображение русского текста в консоли
 
    //Создание типа указателя на функцию
    typedef int (APIENTRY *func1)
    (
        float a,float b,float c,float gran1, float gran2,float e,int z
    );
 
    //Инициализация библиотеки и обозначение этой области памяти
    HINSTANCE hLib=LoadLibrary("resh.dll");
    
    //Сообщение об ошибке если библиотека не инициализировалась
    if(!hLib)
    {
        MessageBox
        (
            NULL, "Библиотека resh.dll - не найдена! "," Ошибка",MB_OK | MB_ICONERROR
        );
    }
    
    //Создание указателя на фенкцию в библиотеке
    func1 Mess=( int(APIENTRY*)
    (
        float a,float b,float c,float gran1, float gran2,float e,int z
    ))GetProcAddress(hLib,"ras4et");
 
    //Вывод сообщения об ошибке если в библиотеке нет этой функции
    if(Mess==NULL)
    {
        MessageBox
        (
            NULL, "Функция не найдена!"," Ошибка",MB_OK | MB_ICONERROR
        );
    }
    else
    { 
         //Ввод данных и использование функции из библиотеки 
         int i;
         i=2;       
         while (i)
         {
               vvz();
               vvabc();
               vvgran();
               vve();
    
               i=Mess(*pa,*pb,*pc,*pgran1,*pgran2,*pe,*pz);
       
               if (i==1 || i==5)
               {
                   break;    
               }        
         }  
    }
    
    FreeLibrary(hLib);    
    system("PAUSE");
    return EXIT_SUCCESS;
}
 
//Выбор функции  
void vvz()
{    
    printf("%s\n", "Выберите функцию и введите ее номер");
    printf("1)y=a*x^3+b*x+c 2)y=a*sin(x+b)+c 3)y=a*ln(abs(x+b))+c\n"); 
 
    int zu;
    bool success = false;
    while(!success)
    {
         if ( vvod() )
         {
              zu=int(af);
              if (zu>0 && zu<4)
              {
                   printf("%s\n","Вы ввели номер функции");
                   *pz=zu;
                   success=true;    
              }
              else
              {
                   printf("%s\n","Значение введено неверно. Повторите ввод.1");
              }          
         }
         else
         {
             printf("%s\n","Значение введено неверно. Повторите ввод.1");
         }
    }    
}   
 
//Ввод А В С 
void vvabc()
{
    printf("%s\n","Введите значения коэффициентов А, В и С");     
    int l;
    l=0;
    float au;
    float bu;
    float cu;
    bool success = false;
    while(!success)
    {
        if (l==0)
        {    
             printf("%s\n","Введите А");
             
             if ( vvod() ) 
             {
                 au=af;             
                 l++;
             }
             else
             {
                 printf("%s\n", "Значения введены неверно. Повторите ввод.2");
             }
        }
        
        if (l==1)
        {    
             printf("%s\n","Введите В");
             
             if ( vvod() ) 
             {
                 bu=af;             
                 l++;
             }
             else
             {
                 printf("%s\n", "Значения введены неверно. Повторите ввод.2");
             }
        }
                
        if (l==2)
        {    
             printf("%s\n","Введите C");
             
             if ( vvod() ) 
             {
                 cu=af;             
                 l++;
             }
             else
             {
                 printf("%s\n", "Значения введены неверно. Повторите ввод.2");
             }
        }
        
                        
        if(l==3)
        {
             printf("%s\n","Вы ввели значения А, В и С");  
             success = true;
        }
    }
    *pa=au;
    *pb=bu;
    *pc=cu;     
}       
 
//Ввод границ   
void vvgran()
{
    printf("%s\n","Введите граничные значения х1 и х2");  
    
    float x1u;
    float x2u;
    
    int l;
    l=0;
    bool success = false;
    while(!success)
    {
        if (l==0)
        {           
            printf("%s\n","Введите х1");
            if ( vvod() ) 
            {
                x1u=af; 
                l++;
            }
        }    
        
        if (l==1)
        {           
            printf("%s\n","Введите х2");
            if ( vvod() ) 
            {
                x2u=af; 
                l++;
            }
        }
        
        if (l==2 && x1u==x2u)
        {
            printf("%s\n","Значения х1 и х2 не должны совпадать.");
            l=1;     
        }
        
        if(l==2)
        {
             printf("%s\n","Вы ввели значения х1 и х2"); 
             if (fabs(x1u-x2u)>100 || fabs(x1u-x2u)<0.1)
             {
                  int p;
                  printf("%s\n","Но ваш интервал слишком большой или маленький.");
                  printf("%s\n","Вы уверены что хотите использовать эти значения?");
                  printf("%s\n","1)Да 2)Введу другие");
                  
                  if ( vvod() )
                  {
                       p=int(af);
                  }
                  
                  if (p==1)
                  {
                       success = true;
                  }
                  else
                  {
                       l=0;
                  }                                                 
             }
             else
             {
                  success = true;
             }
        }
        else 
        {
             printf("%s\n","Значения введены неверно. Повторите ввод.3");              
        }                           
    }
    
    if (x2u<x1u)
    {
        float p;
        p=x2u;
        x2u=x1u;
        x1u=p;       
    }
    
    *pgran1=x1u;
    *pgran2=x2u;                 
}  
  
//Ввод Е  
void vve()
{
    float eu;
    bool success = false;
    while(!success)
    {
        printf("%s\n","Введите значение погрешности Е");
        if ( vvod() )
        {
             eu=fabs(af);
             if (eu>0)
             {
                  printf("%s\n","Вы ввели значение Е");                  
                  if (eu<0.000000001 || eu>0.001)
                  {
                      int p;              
                      printf("%s\n","Но ваше значение слишком маленькое или слишком большое.");
                      printf("%s\n","Вы уверены что хотите использовать это значение?");
                      printf("%s\n","1)Да 2)Введу другое");
                      
                      if ( vvod() )
                      {
                          p=int(af);
                      }
 
                      if (p==1)
                      {
                          success = true;
                      }
 
                  }
                  else
                  {
                      success = true;
                  }    
                                                       
             }
             else
             {                                   
                  printf("%s\n","Значение введено неверно. Повторите ввод.4");     
             }
        }
        else
        {   
            printf("%s\n","Значение введено неверно. Повторите ввод.4");
        }
    }    
    *pe=eu;                   
}  
 
//Функция считывания числа для того чтобы избежать неадекватно введенные данные
int vvod()
{
     char s[100];
     int er;
     er=0;
     scanf("%s",s);
     float aff;
     aff=atof(s);
     af=aff;
     if (aff==0)
     {
          int i;
          i=0;
          while (s[i]==' ')
          {
               i++; 
          }
          if (s[i]!='0')
          {
              er=1;          
          }            
     }           
     
     if (er==1)
     {
          return 0;      
     }
     else
     {
          return 1;
     }
     
}

//Экспорт функции
extern "C" __declspec (dllexport) int ras4et
(float a,float b,float c,float gran1, float gran2,float e,int z);
 
#include <windows.h>
#include <graphics.h>
#include <math.h>
#include <locale.h>
#include <sstream>
using std::stringstream;
#include <iostream>
using std::cout;
using std::endl;
using std::cin;
 
//Здесь объявлены указатели на параметры функции
//границы, погрешность и используемое уравнение функции
float *pa = new float;
float *pb = new float;
float *pc = new float;
float *pgran1 = new float;
float *pgran2 = new float;
float *pn1 = new float;
float *pn2 = new float;
float *pm1 = new float;
float *pm2 = new float;
float *pe = new float;        
int *pz = new int;  
//float to char
char *f2c(float a)
{
    char *s = new char;    
    stringstream iostr;
    iostr << a;
    iostr >> s;
    return s;
}
//End of Float to Char
 
//Уравнение графика функции
float func (float l)
{
     if (*pz==1)
     {
         return (*pa)*l*l*l+(*pb)*l+(*pc);       
     }
     if (*pz==2)
     {
         return (*pa)*sin(l+(*pb))+(*pc);       
     }
     if (*pz==3)
     {
         return (*pa)*log(fabs(l+(*pb)))+(*pc);       
     }                
}      
 
//Здесь объявлены коэффициенты для построения графика в графическом окне
float *pkx = new float; 
float *pky = new float;
 
int *pcx = new int;
int *pcy = new int;
 
//Объявление функций построения графика
//и вывода при большом количестве попыток
int zapros(int number);
void graf(float x1, float x2, float zerno);
 
//Корень, возвращаемое значение
float *pxcor = new float;
int *pflag = new int;
 
//Основная функция         
int ras4et (float a,float b,float c,float gran1, float gran2,float e,int z)
{
    setlocale(LC_ALL,"Russian");//Позволяет отображать русский шрифт в консоле нормально
    *pa=a;
    *pb=b;
    *pc=c;
    *pn1=gran1;
    *pn2=gran2;
    *pgran1=gran1;
    *pgran2=gran2;
    *pe=e;
    *pz=z;
    initwindow(640,480);
 
    float *pzerno = new float;     
    int *pt = new int;
    *pt=0;
 
    {
        *pcx=5;
        *pcy=240;       
        float *pminx = new float;
        float *pmaxx = new float;
        float *pmaxy = new float;
        float *pminy = new float;
        float *ptek = new float;
        float *t = new float;
           
        *pminx=(*pgran1);
        *pmaxx=(*pgran2);
        *pminy=func(*pgran1);
        *pmaxy=*pminy;
        *pzerno=((*pgran2)-(*pgran1))/640;
        if (*pzerno<0.0000001)
        {
            *pzerno=0.0000001;                  
        }        
        
        for (*t=(*pgran1); *t<=(*pgran2); *t+=(*pzerno))
        {
            *ptek=func(*t);
            if (*ptek>*pmaxy) 
            {
               *pmaxy=*ptek;                       
            }
            if (*ptek<*pminy)
            {
               *pminy=*ptek;        
            }            
        }
        *pm1=*pmaxy;
        *pm2=*pminy;
        
        *pmaxx=fabs(*pmaxx-*pminx);
        *pmaxy=fabs(*pmaxy-*pminy);
 
        *pky=470/(*pmaxy);
        *pkx=630/(*pmaxx);
        *pcy=5-int((*pminy)*(*pky));        
        *pcx=(*pcx)-int((*pgran1)*(*pkx));
        delete pminx;
        delete pmaxx;
        delete pmaxy;
        delete pminy;
        delete ptek;
        delete t;
    }
    
    int *pst = new int;
    *pst=1;
    *pflag=0;
 
    printf("\n\n\n\n\n\n\n\n%s\n", "На этом интервале указана середина, которая будет проверена на");  
    printf("%s\n", "условие корня, и если условие не выполнится, то будет взят"); 
    printf("%s\n", "новый интервал, который ближе к этому корню и находится между");
    printf("%s\n", "вертикальными линиями.");
 
//Построение графика и принятие результатов подсчета               
    do 
    {
        (*pt)++;
        graf(*pgran1,*pgran2,*pzerno);
 
        printf("%s %d \n","Текущая итерация",*pt);   
        if (fabs((*pxcor))<0.000001 || fabs(func(*pxcor))<0.000001)
        {
            printf("Текущий х=%e\tу=%e\n",*pxcor,func(*pxcor));                 
        }
        else
        {
            printf("Текущий х=%f\tу=%f\n",*pxcor,func(*pxcor));
        }
                             
        system("PAUSE");
        if (*pflag==2) 
        {
            *pgran1=*pxcor;
        }           
        if (*pflag==3)
        {
            *pgran2=*pxcor;
        } 
           
        if ((*pt)>=25)
        {
            *pst=zapros(*pt);  
        }
        
        if (*pflag==1)
        {
            break;
        }
    }while (*pst==1);
    delete pst;    
 
//Если действие завершается с найденным корнем то он выводится      
    if (*pflag==1)
    {
        printf("\n\n\n\n\n%s", "Программа нашла корень х=");
        if (fabs((*pxcor))<0.000001)
        {
            printf("%e\n", *pxcor);
        }
        else
        {
            printf("%f\n", *pxcor);            
        }
    }
 
    closegraph();
    return(*pflag);
    delete pa;
    delete pb;
    delete pc;
    delete pgran1;
    delete pgran2;
    delete pe;
    delete pz;
    delete pkx;
    delete pky;
    delete pcx;
    delete pcy;
    delete pxcor;
    delete pflag;
    delete pt;
}
 
//Если слишком много действий предлагается завершить программу
int zapros(int number)
{   
    printf("\n\n\n\n\n\n%s\n", "Программа не нашла решение этой функции с введенными вами");
    printf("%s","параметрами за ");
    printf("%d",number);
    printf("%s \n", " попыток.");
    printf("%s\n", "Если вы хотите продолжить расчеты введите 1");
    printf("%s\n", "Также вы можете прервать работу программы");
    printf("%s\n", "и ввести другие данные, для этого введите 2");
    printf("%s\n", "Для завершения работы программы введите 3");
    
    int *pref = new int;
    *pref=3;
    scanf("%d",&*pref);
    if (*pref==1) 
    {
        return 1;
    }
    else
    {
        if (*pref==2) 
        {
            *pflag=4;         
            return 2;
        }
        else
        {
            *pflag=5;
            return 3;
        }    
    }
    delete pref;
}
 
//Построение графика и расчет коэффициентов
void graf(float x1, float x2, float zerno)
{
    setcolor (2);//Цвет для штрихов
    line(128, 450, 128, 480);//Штрихи ОХ
    line(256, 450, 256, 480);
    line(384, 450, 384, 480);
    line(512, 450, 512, 480);
    
    //Шкала делений для Х
    setcolor (4);
    stringstream iostr;
    char *s = new char;
    iostr << ((*pn2-*pn1)/5+*pn1);
    iostr >> s;
    outtextxy(130,450,s);
    
    char *s2 = new char;    
    stringstream iostr2;
    iostr2 << ((*pn2-*pn1)*0.4+*pn1);
    iostr2 >> s2;
    outtextxy(258,450,s2);
    
    char *s3 = new char;
    stringstream iostr3;
    iostr3 << ((*pn2-*pn1)*0.6+*pn1);
    iostr3 >> s3;
    outtextxy(386,450,s3);
    
    char *s4 = new char;
    stringstream iostr4;
    iostr4 << ((*pn2-*pn1)*0.8+*pn1);
    iostr4 >> s4;
    outtextxy(514,450,s4);
    line(0,475,640,475);
    outtextxy(620,455,"X");
    // Конец шкалы для OХ
    
    //Начало шкалы OY
    setcolor (2);//Цвет для штрихов
    line(0, 96, 20, 96);//Штрихи ОY
    line(0, 192, 20, 192);
    line(0, 288, 20, 288);
    line(0, 384, 20, 384);
    
    setcolor (4);
    stringstream iostrY;
    char *ys = new char;
    iostrY << ((*pm2-*pm1)/5+*pm1);
    iostrY >> ys;
    outtextxy(10,98,ys);
    
    char *ys2 = new char;    
    stringstream iostrY2;
    iostrY2 << ((*pm2-*pm1)*0.4+*pm1);
    iostrY2 >> ys2;
    outtextxy(10,194,ys2);
    
    char *ys3 = new char;
    stringstream iostrY3;
    iostrY3 << ((*pm2-*pm1)*0.6+*pm1);
    iostrY3 >> ys3;
    outtextxy(10,290,ys3);
    
    char *ys4 = new char;
    stringstream iostrY4;
    iostrY4 << ((*pm2-*pm1)*0.8+*pm1);
    iostrY4 >> ys4;
    outtextxy(10,386,ys4);
    outtextxy(10,0,"Y");
    //Конец OY
    
    setcolor (15);
    moveto((*pcx)+int((*pn1)*(*pkx)),480-((*pcy)+int(func(*pn1)*(*pky))));
    lineto((*pcx)+int((*pn1)*(*pkx)),480-((*pcy)+int(func(*pn1)*(*pky))));
    
    float *t = new float;
    for (*t=(*pn1)+zerno; *t<(*pn2); *t+=zerno)
    {
        lineto((*pcx)+int((*t)*(*pkx)),480-((*pcy)+int(func(*t)*(*pky))));
    }
 
    setcolor (8);
    moveto((*pcx)+int(x1*(*pkx)),1);
    lineto((*pcx)+int(x1*(*pkx)),480);
    moveto((*pcx)+int(x2*(*pkx)),1);
    lineto((*pcx)+int(x2*(*pkx)),480);  
    
    setcolor (8);//Меняем Координатную линию
    moveto(1,480-(*pcy));
    lineto(640,480-(*pcy));
 
    setcolor (6);//Отвечает за текущую линию
    moveto((*pcx)+int(((x2-x1)/2+x1)*(*pkx)),1);
    lineto((*pcx)+int(((x2-x1)/2+x1)*(*pkx)),480);
    
    *pxcor=x1+(x2-x1)/2;
//    setcolor(1);
    
//    outtextxy(60,20,xcount);
    
    if (fabs(func(*pxcor))<=(*pe))
    {
        *pflag=1;        
    }
    else 
    {
         if (func(x2)*func(x1)>0)
         {
            if (   fabs( func((x2-x1)/4+x1) )   >   fabs( func((x2-x1)*3/4+x1) )   )
            {
               *pflag=2;
            }
            else
            {
               *pflag=3;
            }
         }
         else
         {
             if (( func(*pxcor)>0&&func(x1)>0 ) || ( func(*pxcor)<0&&func(x1)<0 ))
             {
                 *pflag=2;
             }
             else
             {
                 *pflag=3;
             }
         }    
    }    
}