求【0,2】 定积分 ∫（x2+4X+7）dx ，用梯形法，计算精度为10-3.怎么做

t_imba

program ex15
    implicit none
    real::h,f1,f2,s1=0.0,si,s2=0.0
        integer n,i,j
        read*,n
    h=2./n
        f1=7.0
    si=0.0
        
        do i=1,n
    f2=(i*h)**2+4*(I*h)+7
        si=(f1+f2)*h/2.0
    f1=f2
        s1=s1+si
        s2=s1-si
        
        if (abs(s1-s2)<1e-3) exit
        end do
    print *,s1,i,s1-s2

End
不懂如何把求精度的弄进来，请问N必须自己输入么？如果不自己输入怎么做。求指教

石破天惊

abs(s1-s2)<1e-3这个就是对精度的判定了，你不需要自己输入n，只要在外层加一个循环，如果精度不够的话就让n+1，然后再次进行累加求和

t_imba

石破天惊 发表于 2014-4-30 01:11
abs(s1-s2)

能帮改改么？好像直接再套一个,但是N未定义求不出来

石破天惊

program ex15
    implicit none
    real::h,f1,f2,s1=0.0,si,s2=0.0
        integer n,i,j
        n=10
100 continue
    h=2./n
        f1=7.0
    si=0.0
      
        do i=1,n
    f2=(i*h)**2+4*(I*h)+7
        si=(f1+f2)*h/2.0
    f1=f2
        s1=s1+si
        s2=s1-si
        
        if (abs(s1-s2)<=1e-3) then
        goto 200
        elsefi(abs(s1-s2)>1e-3)then
        n=n+1
        goto 100
        endif
        end do
200 continue
    print *,s1,i,s1-s2

End

t_imba

石破天惊 发表于 2014-4-30 01:16
program ex15
    implicit none
    real::h,f1,f2,s1=0.0,si,s2=0.0

额，可以了。我研究研究，谢谢了

石破天惊

诶，不对，不好意思，应该是这样


program ex15
    implicit none
    real::h,f1,f2,s1=0.0,si,s2=0.0
        integer n,i,j
        n=10
100 continue
    h=2./n
        f1=7.0
    si=0.0
      
        do i=1,n
        f2=(i*h)**2+4*(I*h)+7
        si=(f1+f2)*h/2.0
        f1=f2
        s1=s1+si
        s2=s1-si
        enddo
        if (abs(s1-s2)<=1e-3) then
        goto 200
        elsefi(abs(s1-s2)>1e-3)then
        n=n+1
        goto 100
        endif
200 continue
    print *,s1,i,s1-s2

t_imba

石破天惊 发表于 2014-4-30 01:19
诶，不对，不好意思，应该是这样

额，好像两个值都不对，应该定积分的 值应该是在24.666667这样，S1的值都太离谱了

t_imba

石破天惊 发表于 2014-4-30 01:19
诶，不对，不好意思，应该是这样

请问如果不用GOTO语句还有更好的方法么？

石破天惊

f1为什么是7啊，这个不对吧

t_imba

石破天惊 发表于 2014-4-30 01:26
f1为什么是7啊，这个不对吧

哦。程序写错了。是x2+4X+7