clcclearclose?all;
tic
V=2.4000;
U=1.6453;
W=1.7473;

Npoint=501;
Rx=2;
Ry=2;

x=linspace（-RxRxNpoint）;
y=linspace（-RyRyNpoint）;
X=meshgrid（xy）;
Y=meshgrid（yx）;
Y=Y‘;
R=sqrt（X.^2+Y.^2）