C语言实现兰姆波频散特性分析代码

#include 
#include 
#include 
#include 
#include 
#include 
#include 

#define PI 3.1415926
//  double cd=5.8e3,ct=3.1e3;                //钢
// double cd=6.42e3,ct=3.04e3;                   //铝
double cd=2.68e3,ct=1.1e3;    //有机玻璃   

int BisectRoot(double a,double b,double h,double eps,double fb,double* x,int n,int *m);
double Func(double cp,double fb);

void main()
{
        int i,j,n=20,m,N=300;
        double a,b,h,eps,*x,*cp,*dcp,*cg,fb,D=1e-3,step=0.02e3;
        FILE *fp;
        char str[10],str1[10],cgfile[50][10];
    printf("Input the cg file name,please:\n";
        scanf("%s",str1);

        for(j=0;j         { 
                sprintf(str,"%d.dat",j);
                strcpy(cgfile[j],str1);
                strcat(cgfile[j],str);
        }

    a=1.0e0;
        b=15.0e3;
        h=0.7;
        eps=1e-6;        
       
        for (j=0;j<10;j++)
        {
                if((fp=fopen(cgfile[j],"w")==NULL)
                { printf("cg file %s can't open.\n",cgfile[j]); exit(0);}
                else printf("file %s is written.\n",cgfile[j]);

        cp=(double*)calloc(N,sizeof(double));
            if(cp==NULL)        exit(1);
                dcp=(double*)calloc(N,sizeof(double));
            if(dcp==NULL)        exit(1);
        cg=(double*)calloc(N,sizeof(double));
            if(cg==NULL)        exit(1);

                m=0; i=0; fb=0.0;   
                do
                {
                        x=(double*)calloc(n,sizeof(double));
                if(x==NULL)                exit(1);
                        
                BisectRoot(a,b,h,eps,fb,x,n,&m);
                        cp=x[j];

                        if(cp>0 && i>2)
                        {