请教个C++的问题

#include<iostream.h>
#include"math.h"
double fum(double a,double b,double c)
{
if(b*b-4*a*c<0)
return 0;
else
return (-b-(sqrt)(b*b-4*a*c))/2*a;
}
double fun(double m,double n,double p)
{
if(n*n-4*m*p<0)
return 0;
else
return (-n+(sqrt)(n*n-4*m*p))/2*m;
}
void main()
{
cout<<"这是求ax平方+bx+c=0方程的根的公式："<<endl;
double x,y,z;
cin>>x>>y>>z;
cout<<"两根为:"<<fum(x,y,z)<<"或"<<fun(x,y,z)<<endl;
}

自己百度相关的算法,牛顿法是其中一个