Limits should be [0, pi/2] & fun should be sin
fun = @(theta) sin(theta);
S = 4*pi*R^2*(integral(fun,theta_1,theta_2/2)); % this will give correct answer
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!