求解一元四次方程

//本程序基于.NET平台，编译环境是Microsoft Visual C++ .NET
//对于VC++6.0平台的修改，只须把stdafx.h换成iostream.h，然后去掉代码中的“using namespace std”
#include "stdafx.h"
#include "math.h"
using namespace std;
double du,dv;
void solu_w(double a,double b,double c)
{
    double tem,tem1;
    tem1=b*b-4*a*c;
    double m,n;
    m=(-1)*b/(2*a);
    if (tem1>=0)
    {
        tem=sqrt(tem1);
        n=tem/(2*a);
        cout<<"/n方程根为："<<'/n'<<"/t X1 = "<<m+n<<'/n'<<"/t X2 = "<<m-n<<endl;
    }
    else
    {
        tem=sqrt((-1)*tem1);
        n=tem/(2*a);
        cout<<"/n方程根为："<<'/n'<<"/t X1 = "<<m<<" + "<<n<<" i"<<'/n'<<"/t X2 = "<<m<<" - "<<n<<" i"<<'/n'<<endl;
    }
}
void get_du_dv(double a,double b,double c,double d,double e,double f)
{
    dv=(b*d-a*e)/(c*e-b*f);
    du=(c*d-a*f)/(b*f-c*e);
}
void main()
{
    cout<<"输入四次方程的各个系数"<<endl;
    double a4,a3,a2,a1,a0;
    cout<<"  a4 = ";
    cin>>a4;
    cout<<"  a3 = ";
    cin>>a3;
    cout<<"  a2 = ";
    cin>>a2;
    cout<<"  a1 = ";
    cin>>a1;
    cout<<"  a0 = ";
    cin>>a0;
    cout<<"/n方程为：/n/t "<<a4<<" (X*X*X*X) + "<<a3<<" (X*X*X) + "<<a2<<" (X*X) + "<<a1<<" (X) + "<<a0<<" = 0"<<'/n'<<endl;
    double u,v;
    cout<<"输入因子的常数项和一次项系数："<<endl;
    cout<<"  u = ";
    cin>>u;
    cout<<"  v = ";
    cin>>v;
    cout<<"/n因子式为：/n/t w(X) = "<<"(X*X) + "<<u<<" (X) + "<<v<<endl;
    double p0,p1,p2,r0,r1;
    p2=a4;
    p1=a3-u*p2;
    p0=a2-v*p2-u*p1;
    r0=a1-v*p1-u*p0;
    r1=a0-v*p0;
    cout<<"/n/t p(X) = "<<p2<<" (X*X) + "<<p1<<" (X) + "<<p0<<" = 0"<<endl;
    double pp3,pp2,pp1,pp0,r0v,r1v,r0u,r1u;
    pp3=(-1)*p2;
    pp2=(-1)*p1;
    pp1=(-1)*p0;
    pp0=0;
    r0v=u*p2-p1;
    r1v=v*p2-p0;
    r0u=r1v-u*r0v;
    r1u=(-1)*v*r0v;
    get_du_dv(r0,r0u,r0v,r1,r1u,r1v);
    solu_w(1,u+du,v+dv);
}