// Здесь ещё забыли инициализацию нулевых элементов массива: wx[0] = 0.0; // (например) wt[0] = 0.0; // (например) for (int i = 1; i < Nx+1; i++) { wx[i+1] = wx[i] + hx; } for (int j = 1; j < Nt+1; j++) { wt[j+1] = wt[j] + ht; }