Add RK6 with residual viscosity update and reduce computational effort of spatial scheme - Add RK6 with residual updates - Change the D2 operator for upwind schemes to one less computationally expensive.

function z = gaussian(x,x0,d)
    z = exp(-sum((x-x0).^2,2)/d^2);
end