回复： 【求助】用FORTRAN语言写出他的不定积分! 急急急急
拜託，你也去看一下fortran的書，這是最基本的東西，不要隨便上來丟個問題，開平方都不會，回去做點功課吧！
回复： 【求助】用FORTRAN语言写出他的不定积分! 急急急急
呵呵，我连题的意思都没看懂，“（1，2）”是干什么的？元旦愉快回复： 【求助】用FORTRAN语言写出他的不定积分! 急急急急
还有，指定用哪种积分方法了吗？回复： 【求助】用FORTRAN语言写出他的不定积分! 急急急急
puti0772 ，看得出你很着急，但我不知道你是想要fortran内部的系统函数还是要自己编的程序，我这有一个龙贝格积分的作业，不知对你是否会有帮助，但不要输入太高的精度，因为它的系统开销很大，我的机器在迭代超过27次时将发生溢出。我将“（1，2）”设为积分的上下限，在VF中运行时，a,b，E分别对应的是积分得上下限和精度。program suanfa implicit none
real(8), allocatable::T(:)
real(8) TT !save the value of T(0)
integer m
real a,b,E
real y1,y2,yy
integer i,j,k,Num,number
print*,'please input the value of a,b,E:'
read *,a,b,E
yy=0
y1=3*tan(sqrt(a)+2)
y2=3*tan(sqrt(b)+2)
m=0
TT=0
!main process
do while(1)
allocate(T(0:m))
T(0)=(a-b)*(y1+y2)/2
if(m>0)then
do 50 k=0,m-1
do 200 i=0 ,2**k-1
yy=yy+3*tan(sqrt(a+(b-a)*(i+1/2)/(2**k))+2)
200 end do
T(k+1)=T(k)/2+(b-a)*yy/(2**(k+1))
yy=0
50 end do
!**********
do 300 j=m-1,0,-1
T(j)=((4.0**m)*T(j+1)-T(j))/(4.0**m-1)
300 end do
end if
if (abs(T(0)-TT)<E)then
print*,'I(f)=',T(0)
stop 'ok'
end if
m=m+1
TT=T(0)
deallocate(T)
end do
end program
回复： 【求助】用FORTRAN语言写出他的不定积分! 急急急急
上文中program suanfa时程序的第一行，那个笑脸标记是T(:)
回复： 【求助】用FORTRAN语言写出他的不定积分! 急急急急
怎么又这样了，我要晕了，笑脸是冒号和右括号回复： 【求助】用FORTRAN语言写出他的不定积分! 急急急急
禁止在这个帖子中使用笑脸标记可以在编辑帖子时取消：）
还有，请楼主注意斑竹的置顶贴，不要重复发帖
回复： 【求助】用FORTRAN语言写出他的不定积分! 急急急急
不好意思喔！他要的是不定積分，這在fortran中是不可能做到的，除非是定積分................
页:
[1]