菜鸟求助，含参积分那个贴看了，不是很懂，求指点。

潮音洞主

  clear all
k=linspace(1,5);
for i=1:length(k)
    kk=k(i);
    fun=inline(subs('sin(k*x)*x^2','k',kk));
    y(i)=quadl(fun,1,5);
end
hold on;
plot(k,y,'r')

求教，k=linspace(1,5); kk=k(i);和inline()函数的作用不甚清楚。
我想我自己的公式代替'sin(k*x)*x^2'吧？积分限怎么改呢？我的积分限是0到350000.
谢谢！

rocwoods

不要再用inline这种老掉牙又低效的函数结构，效率低还不如匿名函数好使。

1. fun = @(k)@(x) sin(k*x).*x.^2;
   
2. k=linspace(1,5);
   
3. y = zeros(size(k));
   
4. for i = 1:length(k)
   
5. y(i) = quadgk(fun(k(i)),0,350000);
   
6. end
   

复制代码

或者用arrayfun，更简洁，一句话：

1. y = arrayfun(@(k)quadgk(@(x) sin(k*x).*x.^2,0,350000),linspace(1,5));

复制代码

潮音洞主

回复 2 # rocwoods 的帖子

非常感谢，可惜用的matlab7.0，没有quadgk。
请问k=linspace(1,5);是什么作用？1,5何意？