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

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