一元高次方程求解

//本程序基於.NET平臺，編譯環境是Microsoft Visual C++ .NET
//對於VC++6.0平臺的修改，只須把stdafx.h換成iostream.h，然後去掉主函數裏的“using namespace std”
//詳情見blog.csdn.net/chinley。
//#include "stdafx.h"
#include "iostream.h"
#include "math.h"
//using namespace std;
double du,dv;

//求解一元二次方程
void solu_w(double b,double c)
{
    double tem,tem1;
    tem1=b*b-4*c;
    double m,n;
    m=(-1)*b/2;

//判斷是實數解還是複數解。
    if (tem1>=0)
    {
        tem=sqrt(tem1);
        n=tem/2;
        cout<<"/n方程根爲："<<'/n'<<"/t X1 = "<<m+n<<'/n'<<"/t X2 = "<<m-n<<endl;
    }
    else
    {
        tem=sqrt((-1)*tem1);
        n=tem/2;
        cout<<"/n方程根爲："<<'/n'<<"/t X1 = "<<m<<" + "<<n<<" i"<<'/n'<<"/t X2 = "<<m<<" - "<<n<<" i"<<'/n'<<endl;
    }
}

//求解二元一次方程組，其中（a、b、c）、（d、e、f）是方程組兩個的係數，一次是常數項、x係數、y係數。
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<<"詳情見blog.csdn.net/chinley。/n";

    int t1,t2,t3;
    double * a;
    cout<<"輸入階數：";//具體原因待查
    cin>>t1;
    a=new double [t1];
    t3=t2=t1;

    cout<<"輸入各次項的係數：/n";
    for (t1;t1>=0;t1--)
    {
        cout<<"/ta["<<t1<<"]=";
        cin>>a[t1];
    }
    cout<<"方程爲：/n/t";
    for (t2;t2>=1;t2--)
    {
        cout<<a[t2]<<" X^"<<t2<<" + ";
    }
    cout<<a[0]<<" = 0"<<endl;
    double u,v;
    cout<<"輸入因子的常數項和一次項係數："<<endl;
    cout<<"/t u = ";
    cin>>u;
    cout<<"/t v = ";
    cin>>v;
    cout<<"/n因子式爲：/n/t w(X) = "<<"X^2 + "<<u<<" (X) + "<<v<<endl;

    double p[98],r0,r1;

    cout<<"/n輸入修正次數： ";
    int count1;
    cin>>count1;

    for (count1;count1>=0;count1--)
    {
        p[t3]=p[t3-1]=0;
        int t4=t3-2;
        for (t4;t4>=0;t4--)
        {
            p[t4]=a[t4+2]-p[t4+2]*v-p[t4+1]*u;
            r0=a[1]-v*p[1]-u*p[0];
            r1=a[0]-v*p[0];
            //下面一行加在程序中，將導致的問題是：如果修正次數太大，程序運行時間大大延長，主要時間消耗在屏幕顯示上。
            //cout<<"/t p(X) = "<<p[2]<<" (X*X) + "<<p[1]<<" (X) + "<<p[0]<<" = 0"<<endl;
            double r0v,r1v,r0u,r1u;
            r0v=u*p[2]-p[1];
            r1v=v*p[2]-p[0];
            r0u=r1v-u*r0v;
            r1u=(-1)*v*r0v;
            get_du_dv(r0,r0u,r0v,r1,r1u,r1v);
        }

        u=u+du;
        v=v+dv;
    }
    solu_w(u,v);
}