@Adri van der Meer 

thank you

@Christopher2222 

tnx

@gkokovidis 

tnx

@acer 

tnx my friend, i got it. 

@acer 

tnx my friend, i have one problem with ur methods, i try roots : 22.371016133870070862, 50.405610019621053775, 61.672103588567668205, but i can't earn lhs(eq)=0 with these roots.why?

for more explanation if we take x=22.371016133870070862 we earn lhs(eq)=8.2424325*10^8 whlie it must be zero. tnx again. 

@acer 

thank u so much... good luck.

@gaurav_rs 

tnx for ur answer my dear friend, but i must find the other roots. (assuming x>0)

@vv 

thank u so much my dear friend for ur excellent answer, can i solve too these equtions with ur method ?

eq1 := diff(u1(x), x, x)+2*x*(diff(u1(x), x))-u1(x)-(int(u1(s)-u2(s), s = 0 .. x)) = 2+x-exp(x)+2*x*exp(x)-cos(x);
eq2 := diff(u2(x), x, x)+diff(u2(x), x)-2*x*u2(x)-(int(u1(s)+u2(s), s = 0 .. x)) = 2*cos(x)-3*x-(1+2*x)*sin(x)-exp(x);
bcs := u1(0) = 1, (D(u1))(0) = 1, u2(0) = 1, (D(u2))(0) = 1

@Kitonum 

Thank u so much my friend. tnx.

@Preben Alsholm 

Hi, tnx for ur answer, the parameters M, R and aplha are constant, u can assume them M=R=alpha=1 or every value else, i know that must solve it with numerical method but how?

@Carl Love 

thanks my friend, You assume theta=0..Pi, but for theta=0..2Pi how? my plot for theta=0..2*Pi is like theta=0..Pi and this is not true. 

@mehran rajabi 

For more explanation my good friend, the range of Z for this problem is between -0.5 to 0.5 and I will write for loop so first it create values of Z between -0.5 to 0.5 and second the values of function f(z) in this interval([-0.5,0.5]) are closer together. I don't know what some of the values of Z create value closer together for function f(z) and I will obtain this values with maple. 

@ThU 

thank you so much my friend, the range of Z for this problem is between -0.5 to 0.5 and I will write for loop so first it create values of Z between -0.5 to 0.5 and second the values of function f(z) in this interval([-0.5,0.5]) are closer together. I don't know what some of the values of Z create value closer together for function f(z) and I will obtain this values with maple. 

@Mac Dude 

thank you so much my friend, my mean f(z), excuse me for this wrong.

@Carl Love 

thank you so much...