__global__ void Initalize(double* u, double* un, int nx, int nt, double dx)

{

int i = blockIdx.x * blockDim.x + threadIdx.x;

int j = blockIdx.y * blockDim.y + threadIdx.y;