实例解析,【例8-1】复化求积示例（具体题目见书本例8-1）。,解：编写如下语句： v=simple(dsolve(Dv=g-c_d/m*v2,v(0)=0,t) % 微分方程解析解 f=eval(t,g,m,c_d),char(v); % 将符号表达式写成匿名函数 g=9.81;m=68.1;c_d=0.25; v_3=subs(v,t,g,m,c_d,3,g,m,c_d) % 第3秒的下降速度 x1=trape(f,0,3,g,m,c_d) % 复化梯形公式求积 x2=simpson(t)f(t,g,m,c_d),0,3) % 复化辛普森公式求积 x3=Cotes(t)f(t,g,m,c_d),0,3) % 复化Cotes公式求积,运行结果：x1=41.9480 x2 =41.9481 x3 =41.9481,【例8-2】自适应求积方法示例（具体题目见书本例8-9）。,解：首先编写如下被积函数描述文件Dampedsinewave.m，程序如下： function f=Dampedsinewave(t,xi) alpha=atan(-xi/sqrt(1-xi2); %
