MATLAB第四次实验报告.docx

MATLAB第四次实验报告1. 选用复合Simpson公式，计算并用Matlab的符号运算工具箱计算其精确值。比较结果，找出问题原理，提出解决问题的方法。程序：主程序：clcclearf=inline(log(1+x).*log(1-x),x)y=simpson(f,-0.9999,0.9999,100)Simpson函数：function s=simpson(f,a,b,n)h=(b-a)/(2*n);s1=0;s2=0;for k=1:n x=a+h*(2*k-1); s1=s1+feval(f,x);endfor k=1:(n-1) x=a+h*k*2; s2=s2+feval(f,x);ends=h*(feval(f,a)+feval(f,b)+4*s1+2*s2)/3;运行结果：y = -1.11042. 求积分方程的数值解和精确解，分析二者的差异。程序：主程序：clcclear