高等数学--定积分的应用，补充符号方程组解法（八）

元素法

①U是与一个变量x的变化区间[a,b]有关的量
②U对于区间[a,b]具有可加性，就是说如果把[a,b]分成若干份，则U相应的分成许多分量，U等于所有分量之和。
③部分量Δu的近似值可以用f(ε)Δx表示
du=f(ε)dx
那么U=∫ f(x)dx
这个方法就是元素法。

例1 计算两条抛物线y^2=x，y=x ^2所围成的图形面积。

syms x y
eq1=y^2==x;
eq2=y==x^2;
solve(eq1,eq2)
hold on
fplot(sqrt(x))
fplot(x^2)

axis([0,2,0,3])

解方程组可知x1=0，x2=1，取横坐标[0,1]为积分区间，把它分成n等分，每一份都可看成是宽dx，高度x^(1/2)-x ^2 窄矩形组成，从而得到：
A = ∫（x^(1/2)-x ^2）dx

int((x^(1/2)-x ^2),x,0,1)

ans=1/3

利用这种思想，还能不规则物体的体积，曲线的弧长度，变化力所做的功等等。

查询官方文档，补充另外一种解方程的写法：

 syms x f(x) 
f(x)=[sqrt(x),x^2];

通过syms设定一个symfun型变量，名字f(x)，把所有的方程都写入f(x)，用[]号与逗号隔开。

这种写法的意思是，f(x)的表达式可以取其中之一，但是无法直接访问单个函数，需要

fbody = formula(f);
fbody(1)
fbody(2)

要求解这个方程组，不能直接solve(f(x))，必须写成

solve(fbody(1)==fbody(2))

这种方法比起前一种写起来更麻烦，它的优点在于自变量变化时，可以把所有可能的结果都直接计算出来，比如x=2时

 f(2)

当x=4时

有多个自变量，就可以写成syms f(x,y) f(x,y)=[g(x,y),h(x,y),T(x,y)] 形式。