function y=SA(fun,p1,p2,numOfLoops)

syms lamda xn x;
fun = subs(fun,x,xn);

lamda = getLamda(fun,p1,p2);


fun = xn+(lamda*fun);

xn = (p1+p2)/2;
x = 0;
tmp = xn;
i = 0;

while((abs(tmp-x) > 0.00001 || abs(subs(diff(fun))) < 1) && i < numOfLoops)

x =subs(fun);
tmp = xn;
xn = x;
i = i+1;
end
y = x;