Richardson extrapolation for Simpson's rule -code matlab
조회 수: 2 (최근 30일)
이전 댓글 표시
I need you to help me with the code for the Richardson extrapolation for Simpson's rule
Does anyone could help me with this code please?
댓글 수: 0
답변 (1개)
bym
2012년 11월 23일
Simpson's rule is a Richardson extrapolation of the trapezoidal rule, the following should get you started:
clc;clear
yLow = sin(0:pi/10:pi);
yHigh = sin(0:pi/20:pi);
iL = pi/10*trapz(yLow);
iH = pi/20*trapz(yHigh);
format long
re = (4*iH-iL)/3; % Richardson extrapolation
err = [iL (iL-2)/2; ...
iH (iH-2)/2; ...
re (re-2)/2];
fprintf('Estimate\t\terror\n')
fprintf('%1.10f\t%1.10f\n',err')
Estimate error
1.9835235375 -0.0082382312
1.9958859727 -0.0020570136
2.0000067844 0.0000033922
댓글 수: 0
참고 항목
카테고리
Help Center 및 File Exchange에서 Numerical Integration and Differential Equations에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!