#include "stdafx.h"
#include "stdio.h"
#include "math.h"
static float d[20],f[20],e1[10],e2[10];
void main()
{float b[10],c[10]; float delta; float eps;
int print=0; float fru0,fru1,s0,s1,frv0,frv1;
 float deltau,deltav;  float u,v;  float a[10];
 int n,m,j; float r0,r1; float r,q;
  int k,i,l1=0; d[0]=0; f[0]=0;
		cout<<”請輸入方程的最高次數n=："<<endl; cin>>"n";
        cout<<"請輸入方程的系數:"<<endl;
	    for(m=0;m<=n;m++)
	        cin>>a[m];
        if(n==1)
           printf("Found root:%.3f\n",-a[1]/a[0]); 
	 else
	 {Cout<<"請輸入精度：eps=";cin>>eps; }
	Cout<<"方程為:"<<endl;
	for(j=0;j<=n;j++)
		printf("+%.3fx^%d",a[j],n-j);                        //輸出方程的表達式
	  cout<<endl;
	if(n%2==0)  m=n/2;
	else 	m=(n/2)+1;
//求方程的高階因子
    for(k=0;k<m;k++)
	{ u=a[n]/a[n-1]-2*k; v=a[n-1]/a[n-2]+4*k;
	do
	{ b[0]=a[0]; b[1]=a[1]-u*b[0];
        for(j=2;j<=n;j++)
		    b[j]=a[j]-u*b[j-1]-v*b[j-2];
	    r0=b[n-1];  r1=b[n]+u*b[n-1];
	    c[0]=b[0];  c[1]=b[1]-u*b[0];
	    for(j=2;j<=n-2;j++)
		  c[j]=b[j]-u*c[j-1]-v*c[j-2];
		s0=c[n-3]; s1=c[n-2]+u*c[n-3];
		fru0=u*s0-s1; fru1=v*s0; frv0=-s0; frv1=-s1;

		deltau=-(frv1*r0-frv0*r1)/(fru0*frv1-fru1*frv0);
		deltav=-(fru1*r0-fru0*r1)/(frv0*fru1-frv1*fru0);
		u=u+deltau;  v=v+deltav; delta=u*u-4*v; 
}   
       while((fabs(deltau)>eps)||(fabs(deltav)>eps));
	    r=-u/2; q=sqrt(fabs(delta))/2;
		cout<<"得到方程的根為:";
        if(delta<0)
		{ printf("%.3f+%.3fi\t",r,fabs(q)); printf("%.3f-%.3fi\n",r,fabs(q)); }
	    else
		{ printf("%.3f\t\t",r+q);  printf("%.3f\n",r-q); }
	}
}