程序给的是数据范围,没有坐标限制范围,所以是软件以适应的,基本一致
xlim([-5 100])ylim([-15 100])
已解决,20行修改为contour(x,y,EP,150); 24行改为quiver(x,y,dx,dy);