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

t_imba

问题少女 发表于 2014-4-30 07:14
可以用do while语句的，当满足什么条件时就执行！

主要是有两个未知变量，一个N，一个I，一个的我还懂。2个用DO WHILE无从下手，求指教

t_imba

program ex15
    implicit none
    real::h,f1,f2,s1=0.0,si,s2
        integer n,i,j
        n=1
                do
                h=2./n
        f1=7.0
        si=0.0
                s1=0.0
                s2=0.0
        i=1
        
        do while(i<=n)
        f2=(i*h)**2+4*(i*h)+7
        si=(f1+f2)*h/2.0
        f1=f2
        s1=s1+si
        s2=s1-si
                i=i+1
                end do
                print*,n
                n=n+1
               
                if(abs(s1-s2)<=1e-3) exit
            enddo
                print*,i,n,s1,s2
      
               
                end
经过大家指教。总算想出不用GOTO怎么写了。谢谢以上回帖的各位

石破天惊

诶，你这个写法挺聪明的

无言吖啊

我是来加油的