function [xphi, phi, xpsi, psi] = Plot_Multiwavelet(H, r, m, n_refine) Compute the point values of the scaling function and wavelet(s) at all points of the form n*dx which fall inside the support of phi.