Gauss积分matlab的程序?

% 高斯-勒让德数值积分

fun=@(x)sin(x).*x.^3 % fun：积分表达式

% a,b：积分上下限

a = 0

b = pi

tol=1e-8 % tol：积分精度，默认1e-6

% 计算求积节点

syms x

p=sym2poly(diff((x^2-1)^(n+1),n+1))/(2^n*factorial(n))

tk=roots(p)% 求积节点

% 计算求积系数

Ak=zeros(n+1,1)

for i=1:n+1

xkt=tk

xkt(i)=[]

pn=poly(xkt)

fp=@(x)polyval(pn,x)/polyval(pn,tk(i))

Ak(i)=quadl(fp,-1,1,tol)% 求积系数

end

% 积分变量代换，将[a,b]变换到[-1,1]

xk=(b-a)/2*tk+(b+a)/2

% 检验积分函数fun有效性

fun=fcnchk(fun,'vectorize')

% 计算变量代换之后积分函数的值

fx=fun(xk)*(b-a)/2

% 计算积分值

ql=sum(Ak.*fx)

quadl(fun,0,pi) % 调用MATLAB内部积分函数检验

1、matlab是一个功能强大的软件，不仅仅在数据处理方面很优秀，在界面编程方面同样优秀，这里简单介绍下matlab界面编程的基础步骤。

2、在打开的matlab程序中，点击new---graphic user interface，打开创建gui向导--我们选择blank gui，创建空白的gui界面---选择左侧我们需要的控件，如下图，我们选择一个button---将控件拖入到gui界面的合适的位置，双击打开设置属性的界面---设计好界面后，我们先不要编写函数内容，先运行界面---他会提醒我们激活界面将保存界面和代码，我们选择yes---输入文件名，点击保存---我们回到界面编辑界面，点击button右键打开右键菜单，点击view callbacks---callback，来跳转到该控件的回调函数---我们在该函数中输入代码---这时，我们运行程序，点击按钮，即可以在命令窗口中看到button执行的效果。

3、先运行界面，使得matlab给我们创建界面的代码，然后在view callback。

---