u[k * nt + 0] = 1;

}

for (int it = 1; it <= nx - 1; it++)

{

u[m * nx + it] = un[m * nx + it - 1] - c*dt / dx*(un[m * nx + it - 1] - un[(m - 1) * nx + it - 1]);

}