梯形积分
梯形积分
提问于 2011-11-29 22:06:26
我想用matlab来代替集成命令，一个基本的自己创建的命令。你知道如何修复这个错误吗？如果我使用Matlab命令，我的算法运行良好，但当我尝试使用我自己创建的积分函数时，并不令人惊讶，它不工作：(
M文件：
function y = trapapa(low, up, ints, fun)
y = 0;
step = (up - low) / ints;
for j = low : step : up
y = y + feval(fun,j);
end
y = (y - (feval(fun, low) + feval(fun, up))/2) * step;均值算法：
clear;
x0=linspace(0,4,3);
y=linspace(0,2,3);
for i=1:length(x0)
for j=1:length(y)
x(i,j)=y(j)+x0(i);
alpha=@(rho)((5-2*x(i,j)).*exp(y(j)-rho))./2;
%int(i,j)=quad(alpha,0,y(j))
int(i,j)=trapapa(alpha,0,y(j),10)
end
end回答 1
Stack Overflow用户
回答已采纳
发布于 2011-11-29 23:25:00
您没有遵循脚本中的函数定义。有趣的参数(变量alpha)应该是最后一个。
试试int(i,j)=trapapa(0,y(j),10,alpha)
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/8312070
复制相关文章
相似问题
腾讯云开发者
