查看: 136|回复: 0|关注: 0

[已答复] 積分式跑不出來

[复制链接]

新手

9 麦片

财富积分


050


2

主题

5

帖子

0

最佳答案
发表于 2019-11-6 21:41:28 | 显示全部楼层 |阅读模式
本帖最后由 alexliang121 于 2019-11-11 14:18 编辑

請問各路高手..
我想跑附圖的積分式!
我的程式哪裡出錯了!
s=real(10);
sprintf ('how many values of M?');
sprintf ('Potential integral');
sprintf ('Evaluated with pulse functions');
sprintf ('M integral extrapolation');
M=1;
n=8;
k=1,n;
for k=1:n;
    s=0;
    for j=1,M;
        for i=1,M;
            s=s+1/(sqrt((i-0.5)*(i-0.5)+(j-0.5)*(j-0.5)))/M*2;
            Q(k,1)=s;
        end
    end
    c=1;
    for i=2:k;
        c=c*2;
        Q(k,i)=(c*Q(k,i-1)-Q(k-1,i-1))/(c-1);
    end
    print (M,Q(k,1),Q(K,i));
    M=M*2;
end



積分產生的數值

積分產生的數值
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /3 下一条

快速回复 返回顶部 返回列表