设为首页收藏本站

EPS数据狗论坛

 找回密码
 立即注册

QQ登录

只需一步,快速开始

查看: 25|回复: 0

关于Matlab中龙贝格代码的实现

[复制链接]

1

主题

0

金钱

9

积分

新手用户

发表于 2021-7-12 11:18:03 | 显示全部楼层 |阅读模式
function [s,k]=Romberg_li(a,b,eps)
if nargin==2
    eps=1.0e-6;
elseif nargin<2
    error
    return
end
t1=10^4;
t2=-10^4;n=2;
t(1,1)=0.5*(b-a)*(f(a)+f(b));area=t(1,1)/2;
while abs(t2-t1)>=eps
    h=(b-a)/2^(n-1);
    for i=1:22^(n-1))
        area=area+h*f(h*i+a);
    end
    t(n,1)=area ;
    for j=2:n
        for i=1:(n-j+1)
            t(i,j)=(4^(j-1)*t(i+1,j-1)-t(i,j-1))/(4^(j-1)-1);
        end
    end
    t1=t(i,n);
    t2=t(i+1,n-1);
    n=n+1;
    area=area/2;
end
s=t1;k=n-1;

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

关闭

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

客服中心
关闭
在线时间:
周一~周五
8:30-17:30
QQ群:
653541906
联系电话:
010-85786021-8017
在线咨询
客服中心

意见反馈|网站地图|手机版|小黑屋|EPS数据狗论坛 ( 京ICP备09019565号-3 )   

Powered by BFIT! X3.4

© 2008-2028 BFIT Inc.

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