作为我的代码的一部分，我通过积分计算2d弧长来构建用于网格划分的几何。 下面是一个示例代码：
fun=@(ss)integral(@(t)sqrt((r0+a1*cos(1*t)+a2*cos(2*t)+a3*cos(3*t)+a4*cos(4*t)+a5*cos(5*t)+a6*cos(6*t)+a7*cos(7*t)+a8*cos(8*t)).^2 ...
+ (-a1*1*sin(1*t)-a2*2*sin(2*t)-a3*3*sin(3*t)-a4*4*sin(4*t)-a5*5*sin(5*t)-a6*6*sin(6*t)-a7*7*sin(7*t)-a8*8*sin(8*t)).^2),0,ss);
整合是在＆＃39;所有其他参数都是常数（即r0，a1等）
但是集成需要很长时间，而且由于我尝试了许多不同的几何形状，如果我能找到一种方法来加速这种需要一段时间的整合，那将是很好的。
如何以更快的方式计算相同的积分？