function [I,step]=Roberg(f,a,b,eps)
if(nargin==3)
eps=1.0e-4;
end;
M=1;
tol=10;
k=0;
T=zeros(1,1);
h=b-a;
T(1,1)=(h/2)*(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b));
while tol>eps
k=k+1;
h=h/2;
Q=0;
for i=1:M
x=a+h*(2*i-1);
Q=Q+subs(sym(f),findsym(sym(f)),x);
end
T(k+1,1)=T(k,1)/2+h*Q;
M=2*M;
for j=1:k
T(k+1,j+1)=T(k+1,j)+(T(k+1,j)-T(k,j))/(4^j-1);
end
tol=abs(T(k+1,j+1)-T(k,j));
end
I=T(k+1,k+1);
step=k;
我想显示出T表请问怎么样做
matlab计算过程显示
答案:2 悬赏:30
解决时间 2021-12-29 10:38
- 提问者网友:恋你成殇
- 2021-12-28 14:22
最佳答案
- 二级知识专家网友:魅世女王
- 2021-12-28 14:57
在输出那块多加一个T,不行的话就可在函数最后另起一行写T,不带分号,运行的时候应该就能显示了
全部回答
- 1楼网友:转身→时光静好
- 2021-12-28 16:34
通过现有的程序结构是做不到的,需要改写comet函数才行(因为更新绘图的循环过程是在comet函数内部)。
我要举报
如以上问答内容为低俗、色情、不良、暴力、侵权、涉及违法等信息,可以点下面链接进行举报!
大家都在看
推荐资讯