function[m1]=i1sin(a,b,x0,v)
a1=(a-x0)/v;b1=(b-x0)/v;
m1=(x0/pi)*(asin(b1)-asin(a1))+(v/pi)*(sqrt(1-a1^2)-sqrt(1-b1^2));