Метод бисекции для нахождения корней в нелинейном уравнении

#include <iostream>
#include <conio.h>
#include <math.h>
 
 
 
double sign(double mak){
    double x=mak;
    if (x > 0) return 1;
if (x < 0) return -1;
return 0;}
 
using namespace std;
 
double bisection_count(double end , double start , double eps ){
    double res;
    double temp;
    double temp_midle;
    double f_st , f_end;
    f_st=start; 
    f_end= end;
    double f_mid_i;
 
    f_st = (pow(start,4) +3*pow(start,3)-2*start+4);
    f_end = ( pow(end,4) +3*pow(end,3)-2*end+4);
    if(f_st==0) {
        cout<<"Root equels"<<endl;
        return f_st;}
 
    else if (f_end==0) {
        cout<< "Root equels 0"<<endl; 
         return f_end;
    }
 
 
    temp=end-start;
    temp_midle=(end-start)/2;
 
    double f_midle =    ( pow(temp_midle,4) +3*pow(temp_midle,3)-2*temp_midle+4);
 
    while (fabs(f_midle)>eps){
 
        temp=temp/2;
        temp_midle=end+temp;
 
        f_mid_i= ( pow(temp_midle,4) +3*pow(temp_midle,3)-2*temp_midle+4);
 
      if(sign(f_st)!=sign(f_mid_i))
      {end=temp_midle;}
      cout<<"Root is found:";
    }
 
    return temp_midle;
}
 
 
int main(){
double start;
double end;
double eps;
cin>>eps>>start>>end;
double z;
z=bisection_count(end,start,eps);
cout<<"res:"<<z;
getch();
 
    return 0;}