d10r10.cpp

工程计算中经常用到的数值vc算法

# include<math.h>
# include<iomanip.h>
# include<iostream.h>
# include<stdlib.h>
# include<process.h>
# include<stdio.h>

void main()
{
    //program d10r10
    //driver for routine laguer
	int m,mp1,ntry,j,i,n,iflag,polish;
	double eps,aaa,bbb ,x[3];
    double  a[3][6],y[3][22];
	m = 4;
    mp1 = m + 1;
    ntry = 21;
    eps = 0.000001;
    for (j = 1; j<=mp1; j++)
        for (i = 1; i<=2; i++)
            a[i][j] = 0.0;
    a[2][1] = 2.0;
    a[1][3] = -1.0;
    a[2][3] = -2.0;
    a[1][5] = 1.0;
    cout<<endl;
    cout<<"Roots of polynomial x^4-[1+2i]*x^2+2i"<<endl;
    cout<<endl;
    cout<<"            Real     Complex"<<endl;
    n = 0;
    polish = 0;
    for (i = 1; i<=ntry; i++)
	{
        x[1] = (i - 11.0) / 10.0;
        x[2] = (i - 11.0) / 10.0;
        laguer(a,m,x,eps,polish);
        if (n == 0)
		{
            n = 1;
            y[1][1] = x[1];
            y[2][1] = x[2];
			printf("%6.1d %10.4f %10.4f\n",n,x[1],x[2]);
            
		}
        else
		{
            iflag = 0;
            for (j = 1; j<=n; j++)
			{
                aaa = cabs(x[1] - y[1][j],x[2] - y[2][j]);
                bbb = eps * cabs(x[1],x[2]);
                if (aaa <= bbb) iflag = 1;
            }
            if (iflag == 0)
			{
                n = n + 1;
                y[1][n] = x[1];
                y[2][n] = x[2];
				printf("%6.1d %10.4f %10.4f\n",n,x[1],x[2]);
            }
        }
    }
}