#include <math.h>
#include <stdlib.h>
#include <stdio.h>
#include <vector>
#define IE 40
#define JE 40
#define KE 40
#define ia 7
#define ja 7
#define ka 7
#define NFREQS 3
using namespace std;
main()
float l,m,n,i,j,k,ic,jc,kc,nsteps,n_pml;
float ddx,dt,T,epsz,muz,pi,eaf,npml;
int ib,jb,kb;
.
.
.
.
ic=IE*0.5 ;
jc=JE*0.5 ;
....'ic'dose not name a type.
....'jc'dose not name a type.