This sems to be homework.
Vectorise the expoinentiation:
y= (2*x.^1.01)+ sin(3*pi/2 *x)-10*cos(x)-3;
↑ ← HERE
Then to find the approximate indices of the zero-crossings:
zxi = find(diff(sign(y)));
Then use those with interp1 to find much more precise estimates of the zero-crossing x-values.