一、实验目的 许多工程技术和数学研究中要用到定积分,如果无法直接算不出精确值(如含在积分方程中的积分)或计算困难但可用近似值近似时,就用数值积分法方法加以解决。常用的算法有:复化梯形、辛甫生(Simpson)、柯特斯(Cotes)求积法; 龙贝格(Romberg)算法;高斯(Gauss)算法。 二、 ...
分类:
编程语言 时间:
2020-05-12 22:02:32
阅读次数:
140
function Romberg(a,b,n)h=b-a;R(1,1)=h*(f(a)+f(b))/2;for i=2:n m=0; for k=1:2^(i-2) m=m+f(a+(k-0.5)*h); end R(2,1)=1/2*(R(1,1)+h*m);...
分类:
编程语言 时间:
2016-01-11 20:18:18
阅读次数:
303
hdu 1724
题意:
求椭圆面积, x^2/a^2+y^2/b^2=1, 给出a,b,l,r, 求椭圆被直线x=l, x=r所夹面积。思路:
1. 自适应simpson,精度不是很高,所以eps要设小一点。
2. romberg,精度很好,但效率不是很高。自适应simpson:
{CSDN:CODE:617152}romberg:
{CSDN:CODE:617154}...
分类:
其他好文 时间:
2015-03-11 19:42:58
阅读次数:
138