matlab 计算阴影面积,MATLAB求解阴影面积

Figure 1

Matlab solution：

STEP 1. Find the possible intersection points on two circles.

![][01]

[01]: http://latex.codecogs.com/svg.latex?\left{\begin{array}{ll}x2+(y-a)2&=a2\(x-\frac{a}{2})2+(y-\frac{a}{2})2&=(\frac{a}{2})2\end{array}\right.

code:

>> syms x x0 y0 a positive;

>> [x0, y0]=solve('x^2+(y-a)^2=a^2', '(x-a/2)^2+(y-a/2)^2=(a/2)^2');

output:

![][02]

[02]: http://latex.codecogs.com/svg.latex?\left{\begin{array}{ll}x_0&=\frac{5\pm\sqrt{7}}{8}a\y_0&=\frac{3\pm\sqrt{7}}{8}a\end{array}\right.

and the smaller one of x0 is the first cross point.

STEP 2. find the area s0 by integrating function f on (0, x0), where f is