计算物理实验报告数值求积.docVIP

数值求积

实验目的

掌握基本数学运算，学会利用梯形法则、Simpson法则等求积方法，包括了解这些公式得到的结果的精度和误差

熟练运用fortran来解决实际问题

比较梯形法则和Simpson法则的精度

实验问题

写一个程序，用梯形法则和Simpson法则求积分

并且考察它对不同h值的精度

实验方法

梯形法则

梯形法则是采用梯形来估计曲线下方面积，这等同将被积函数近似为直线函数，被积的部分近似为梯形，要求得较准确的数值，可以将要求积的区间分成多个小区间。

Simpson法则

用过三点f-1,f0,f1的抛物线逼近f(__)，对f(__)作泰勒展开

变量代换

令,则积分形式作以下变换

四、实验程序

梯形法则

programt__1

implicitnone

integer::i

real(8)::h,d,s,f,result1

real(8)::sum1,p,n,l,e__act,diff1！准备基本的变量

n=.0

h=1.0/n！取n值，确定格子间隔h

p=h