clear;clc;
ii=sqrt(-1);
r1=1;
nw=1000;
ky=1;
r2=0.5;
mR=0.6;
w=linspace(0.0001,4,nw);
vgt=zeros(1,4);
vg1=zeros(1,nw);