声波方程有限差分数值模拟程序C++

#include 
#include 
#include 
using namespace std;

int main（）
{
	int kijt;
	float ***p;
	float dt=0.02h=0.5v=1;
	float a=（v*dt）/h;
	p=new float **[3];
	for（k=0;k<3;k++）
	{
		p[k]=new float *[101];
		for（i=0;i<101;i++）
		{
			p[k][i]=new float [101];
				for（j=0;j<101;j++）
					p[k][i][j]=0;
		}
	}

	p[1][50][50]=1;

	for（t=0;t<1000;t++）
	{
		for（i=1;i<100;i++）
		{
			for（j=1;j<100;j++）
			{
				p[2][i][j]=2*p[1][i][j]-p[0][i][j]+