function [m0]=i0sin(a,b,x0,v)
if a<=x0-v
a1=-1;
else
a1=(a-x0)/v;
end
if b>=x0+v
b1=1;
else
b1=(b-x0)/v;
end
m0=(1/pi)*(asin(b1)-asin(a1));