I am interested in increasing "MaxPoints" in NDSolve's "MethodOfLines" in attempt to increase the resolution of the plot of the solution of linearly damped wave equation with transparent boundary conditions and square pulse initial condition. Here is my code: interpolatingFunctLinear[initalPulseFunction_, xBoundLow_, xBoundHigh_, timeBound_] := First[ pdeY = D[Y[x, t], t, t] + .04 D[Y[x, t], t] == D[Y[x, t], x, x]; solnDerivativeY = NDSolve[{pdeY, Y[x, 0] == initalPulseFunction, Derivative[0, 1][Y][x, 0] == 0, Derivative[1, 0][Y][xBoundLow, t] == Derivative[0, 1][Y][xBoundLow, t], Derivative[1, 0][Y][xBoundHigh, t] == -Derivative[0, 1][Y][xBoundHigh, t]}, Y, {x, xBoundLow, xBoundHigh}, {t, 0, timeBound }, Method -> {"MethodOfLines", "SpatialDiscretization" -> {"TensorProductGrid", "MaxPoints" -> 2000}}]] Here is my square wave: Piecewise[{{1 , Abs[x] 8.40749}}] Now running the whole thing together, we get the...