pulse.c

利用不均一分割的FDTD解析程序

/*
 *
 *  Incident Wave  (Gaussian Pulse)
 *
 */

double pulse(int n)
{
	double V;
	double alpha = 1.0/64;

	if(n < 65){
		V = exp(-alpha*(pow(32.0-n,2.0)));
	}
	else{
		V = 0.0;
	}
	return V;
}