基于FDTD_1D平板间电磁波仿真(C语言）

//一维FDTD仿真
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define pi 3.1415
#define C 3.0*pow(10,8)
struct parameter
{

double T;

double epsl;
double mu;
double t0;
double dz;
double dt;

};


writedata(double *pe,double*ph,int N)
{
FILE *fp1,*fp2;
int i;
fp1=fopen("outE.dat","w");
for(i=0;i<N;i++)
{
fprintf(fp1,"%f  ",pe[i]);
/***********test*****************************/
printf("Ex:  %f  \n",pe[i]);