c++程序哪里出错了

来源:百度知道 编辑:UC知道 时间:2024/05/21 07:14:41
#include<iostream>
#include<fstream>
#include<cmath>
using namespace std;

const double PI=3.14;
const double c0=3e8;//光速
const double e0=8.85*pow(10,-12);//真空中的电导率
const double u0=4*PI*pow(10,-7);//真空中的磁导率

double L1=300,L3=300; //总长,宽
double T=L1/c0;//总时间

double w=pow(10,7);
double dy=PI*c0/(w*10);//!!!!!!!!!!
double dx=dy;
double dt=dy/(2*c0);
int n=int(T/dt);//时间总格
int Ni=int(L3/dy);//横空间总格100
int Nj=int(L1/dy);//200

int i=0,j=0;//全局变量
double t=0;

double ez[300][300];
double hx[300][300];
double hy[300][300];

double *ezi1=new double[Ni-1];//i,j
double *ezi0=new double[Ni-1];
double *ez1j=new double[Nj-1];
double *ez0j=new double[Nj-1];
double *eznj=new double[Nj-1];//n代表减1
double *ezNj=new double[Nj-1];

double Ei=0;
double Ei1=0;

void FDTD();

void main()
{

扩号位置:

#include<iostream>
#include<fstream>
#include<cmath>
using namespace std;

const double PI=3.14;
const double c0=3e8;//光速
const double e0=8.85*pow(10,-12);//真空中的??率
const double u0=4*PI*pow(10,-7);//真空中的磁?率

double L1=300,L3=300; //??,?
double T=L1/c0;//???

double w=pow(10,7);
double dy=PI*c0/(w*10);//!!!!!!!!!!
double dx=dy;
double dt=dy/(2*c0);
int n=int(T/dt);//???格
int Ni=int(L3/dy);//?空??格100
int Nj=int(L1/dy);//200

int i=0,j=0;//全局?量
double t=0;

double ez[300][300];
double hx[300][300];
double hy[300][300];

double *ezi1=new double[Ni-1];//i,j
double *ezi0=new double[Ni-1];
double *ez1j=new double[Nj-1];
double *ez0j=new double[Nj-1];
double *eznj=new double[Nj-1];//n代表?1
double *ezNj=new double[Nj-1];

double Ei=0;
double Ei1=0;