1.主程序(fortran77 写的):
program vol
implicit double precision (a-h,o-z)
pi=3.14159
read (*,*) r1,r2,h
v=pi/3*h*(r1**2+r2**2+r1*r2)
write(*,*) v
end
2,参数输入
如果主机的系统是linux的话,例如可输入 echo 5,4,3 > par.txt
(其中,5,4,3是R,r,h;par.txt是重复擦写的,不是累加的)
3.参数读出计算
如果主程序的可执行文件是vol.x,那么直接输入vol.x < par.txt ,体积结果就输入输出在屏幕上了。
wjisah ok
已经过点了。
program circle
implicit none
real*8,parameter::pi=3.1415926
real*8 h,r1,r2,v ! fortran 不能区分大下写,因此将R,r定义为r1,r2
write(*,*) 'please input R,r,h.'
read(*,*) r1,r2,h
v=pi*(r1*r1+r2*r2+r1*r2)/h/3.d0
end program