Lax–Friedrichs method