有限差分,beam propagation method,matlab