数值积分算法：复合梯形公式、复合辛普森公式和龙贝格算法

定义函数：函数名为fun41，输入一个变量x，输出变量y，其中y=sqrt(x) * log(x)。
使用复合梯形公式和复合辛普森公式计算数值积分：

定义步长h，将其分为n个小步长。
对于复合梯形公式，初始化t=0，通过循环累加fun41函数在每个小步长上的值，最后根据公式计算积分值T。
对于复合辛普森公式，初始化s1和s2为0，通过两个循环分别计算fun41函数在每个小步长上的值，最后根据公式计算积分值S。

使用龙贝格算法计算数值积分：

初始化m=16，h=1，计算T表的第一项。
循环迭代计算T表的对角线上的元素，每次将h减半，计算新的积分值t，并根据梯形公式计算出新的T表项。
对T表进行迭代，不断更新T表的对角线上的元素，最终得到积分值。

代码示例：

% 1. 定义函数
function y=fun41(x)
y=sqrt(x).*log(x);
end

% 2. 复合梯形公式和复合辛普森公式
clear;
clc;
h=0.001; % 步长
n=1/h;
t=0;
for i=1:n-1
  t=t+fun41(i*h);
end
T=h/2*(0+2*t+fun41(1));
T=vpa(T,10) % 复合梯形公式

% 复合辛普森公式
s1=0;
s2=0;
for i=0:n-1
  s1=s1+fun41(h/2+i*h);
end
for i=1:n-1
  s2=s2+fun41(i*h);
end
S=h/6*(0+4*s1+2*s2+fun41(1));
S=vpa(S,10)

% 3. 龙贝格算法
clear;
clc;
m=16;
h=1;
T(1)=(0+fun41(1))*h/2
for i=2:m
  h=h/2;
  n=1/h;
  t=0;
  for j=1:2:n-1
    t=t+fun41(j*h);
  end
  T(i)=T(i-1)/2+h*t; % 梯形公式
end
for i=1:m-1
  for j=m:i+1
    T(j)=4^i/(4^i-1)*T(j)-1/(4^i-1)*T(j-1);
  end
end
vpa(T(m),10)

本文详细介绍了三种常用的数值积分算法，并提供了相应的代码示例。读者可以根据自己的需求选择合适的算法，并使用代码示例进行实践。