作业帮求高手看一下matlab编的龙贝格求积公式
来源：学生作业帮 编辑：作业帮 分类：综合作业 时间：2024/05/02 20:16:56
求高手看一下matlab编的龙贝格求积公式
function [I,step]=Roberg(f,a,b,eps)
M=1;
h=b-a;
err=1
J=0;
R=zeros(4,4);
R(1,1)=h*(feval('f',a)+feval('f',b))/2
while ((err>delta)&(J
function [I,step]=Roberg(f,a,b,eps)
M=1;
h=b-a;
err=1
J=0;
R=zeros(4,4);
R(1,1)=h*(feval('f',a)+feval('f',b))/2
while ((err>delta)&(J
试试我为你编的这个代码
function I=romberg(f,a,b,e)
%%对给定函数f(x),积分区间[a,b]和误差界e,用Romberg积分法计算积分
T1=0;T2=0;S1=0;S2=0;C1=0;C2=0;R1=0;R2=0;
h=b-a;
T1=h/2*(f(a)+f(b));
ee=1;
while ee>e
for k=1:2
R1=R2;
for j=1:2
C1=C2;
for i=1:2
S1=S2;
s=0;
x=a+h/2;
while x
再问： 那请问function y = DELGKT2_mid(f,h,a,b,y0,varvec) format long; N = (b-a)/h; y = zeros(N+1,1); y(1) = y0; x = a:h:b; var = findsym(f); for i=2:N+1 v1 = Funval(f,varvec,[x(i-1) y(i-1)]); t = y(i-1) + h*v1/2; v2 = Funval(f,varvec,[x(i)+h/2 t]); y(i) = y(i-1)+h*v2; end format short; 有没有错误呢？为什么我用不起来。
function I=romberg(f,a,b,e)
%%对给定函数f(x),积分区间[a,b]和误差界e,用Romberg积分法计算积分
T1=0;T2=0;S1=0;S2=0;C1=0;C2=0;R1=0;R2=0;
h=b-a;
T1=h/2*(f(a)+f(b));
ee=1;
while ee>e
for k=1:2
R1=R2;
for j=1:2
C1=C2;
for i=1:2
S1=S2;
s=0;
x=a+h/2;
while x
再问： 那请问function y = DELGKT2_mid(f,h,a,b,y0,varvec) format long; N = (b-a)/h; y = zeros(N+1,1); y(1) = y0; x = a:h:b; var = findsym(f); for i=2:N+1 v1 = Funval(f,varvec,[x(i-1) y(i-1)]); t = y(i-1) + h*v1/2; v2 = Funval(f,varvec,[x(i)+h/2 t]); y(i) = y(i-1)+h*v2; end format short; 有没有错误呢？为什么我用不起来。