ymin(i,1)=sqrt(rs(i,1)^2-(rg(i,1)-h)^2);
ymax(i,1)=((z0*ampl(i,1)*h)/(1.5*Cfoa(t,1)))*(1+.7*(V(i,1)/V0)*(1/sqrt(1-.5*(V(i,1)/V0)^2)));
if ((y(i,1)==0 | abs(y(i,1))<=ymin(i,1)) & x(i,1)<=0.05*rs(i,1)) %voltage emali bar ruye tire bedune barghgir
Ig=(Eg*ro)/(2*pi*((Ro)^2));
ir=(ampl(i,1))*((Ro*Z)/(Ro+Z))*(1/Ro);
Ri=(Ro/(sqrt(1+(ir/Ig))));
Vpole(m,1)=Varr+ampl(i,1)*Ri;
% Ipolearr=(1.5*Cfoa-Varr)/5;
%ppole=1/(1+(Ipolearr/31)^2.6);
I_pole(m,1)=ampl(i,1);
m=m+1;