[問題] 用梯形法做積分有warning看不懂
各位大大
小弟剛入門fortran
有些經驗還是不夠
找不出來哪邊是錯的
請提點一下
程式如下：
program trapezoid
integer i
real*8 h,sum,n,a,b,ans
read(*,*)a
read(*,*)b
read(*,*)n
h=(b-a)/n
sum=0
x=0
do i=1,n-1
sum=sum+f(x)
x=x+h
enddo
ans=((f(a)+f(b))/2+sum)*h
write(*,*) ans
end
function f(x)
real x
f=x*x
return
end
而跑出來的warning是說
inconsistent data type for argument 1 in call to F
我試一個簡單的積分跑出來都沒有值
看了很久不知道哪邊有問題
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.116.78.69
推
03/19 17:05, , 1F
03/19 17:05, 1F
推
03/19 17:14, , 2F
03/19 17:14, 2F
→
03/19 17:15, , 3F
03/19 17:15, 3F
推
03/19 17:22, , 4F
03/19 17:22, 4F
→
03/19 17:23, , 5F
03/19 17:23, 5F
→
03/19 17:23, , 6F
03/19 17:23, 6F
→
03/19 17:24, , 7F
03/19 17:24, 7F
推
03/19 17:31, , 8F
03/19 17:31, 8F
推
03/19 17:33, , 9F
03/19 17:33, 9F
→
03/19 17:43, , 10F
03/19 17:43, 10F
→
03/19 17:58, , 11F
03/19 17:58, 11F
→
03/21 01:50, , 12F
03/21 01:50, 12F
Fortran 近期熱門文章
PTT數位生活區 即時熱門文章
