#include "iostream.h"
#include "math.h"
#include "stdio.h"

void Initialize(double U[4][N+3][M+3])
{
	int i,j;
	for(i=0;i<N+3;i++)
		for(j=0;j<M+3;j++)
		{
			U[0][i][j]=rou0;
			U[1][i][j]=rou0*u0;
			U[2][i][j]=rou0*v0;
			U[3][i][j]=p0/(k-1)+rou0*(u0*u0+v0*v0)/2.0;
		}
}