x=0:0.5:2
y=sqrt.(1 .+x.^2)
t=trapz(y)
q=quad(u->sqrt(1+u^2),0,2)
qq=quadl(u->sqrt(1+u^2),0,2)
println(t)
println(q)
println(qq)

x=0:0.05:2
y=sqrt.(1 .+x.^2)
t=trapz(y)
q=quad(u->sqrt(1+u^2),0,2)
qq=quadl(u->sqrt(1+u^2),0,2)
println(t)
println(q)
println(qq)

x=0:0.01:2
y=sqrt.(1 .+x.^2)
t=trapz(y)
q=quad(u->sqrt(1+u^2),0,2)
qq=quadl(u->sqrt(1+u^2),0,2)
println(t)
println(q)
println(qq)

x=0:0.005:2
y=sqrt.(1 .+x.^2)
t=trapz(y)
q=quad(u->sqrt(1+u^2),0,2)
qq=quadl(u->sqrt(1+u^2),0,2)
println(t)
println(q)
println(qq)
为什么trapz这个梯形方法随着步长的减小误差成倍数增加，和另外两种方法相比误差极大，是函数本身的问题还是代码问题呢