fortran中write什么意思

F90程序中的数据都有三个特征：类型、种别、属性,由类型说明语句来定义说明.其一般形式是：
类型说明[(种别说明)][,属性说明表] ::变量名表[=初值]
例：REAL(KIND=2),DIMENSION(1:10)::X,Y
说明变量X,Y都是实型,种别是2,属性是一维数组.这样X、Y实际上是两个种别参数为2的一维实型数组,各具有10个元素.在某些场合下,种别说明与属性说明可以省略,此时它的一般形式是最基本形式：
类型说明::变量名表
这个符号是在定义数据类型时所用的符号,我想你应该明白了吧

PROGRAM ex
INTEGER::m(6),i,j,k
call random_seed ()
DO k=1,6,1
call sub1(m(k))
ENDDO
END PROGRAM
SUBROUTINE sub1(n)
REAL::x
INTEGER::n
call random_number (x)
n=nint(x*100)
write(*,*) n
RETURN
END SUBROUTINE
random_seed 只需调用一次既可

我虽然不懂fortran,但我也看出来了：cannot find module TEXTGRAPHLIB
你第二行里用到了TextGraphLib,要引入相关文件.

如果您说的e是那个常数,那请用exp(x)
：）
如果它只是一个变量名,那请用e**x
：）

(x-y)**4既可

你是计算机1091的吧,我知道你了,你准备重修吧.

status
Character*(*).Describes the nature of the file to be opened.Can be one of the following:
●'OLD' The file must already exist.If the file exists,it is opened.If it does not exist,an I/O error occurs.
●'NEW' The file must not already exist.If the file does not exist,it is created.If it exists,an I/O error occurs.
●'SCRATCH' If you omit the file parameter when opening a file,the value of status defaults to 'SCRATCH'.Scratch files are temporary.They are deleted when the unit is closed or when the program terminates.
●'REPLACE' The file opened replaces one of the same name.If there is no existing file of the same name,a new file is created.
●'UNKNOWN' (Default) The run-time system first attempts to open the file with status equal to 'OLD',and then with status equal to 'NEW'.If the file exists,it is opened; if it does not exist,it is created.
●Using STATUS= 'UNKNOWN' avoids any run-time errors associated with opening an existing file using STATUS='NEW' or opening a nonexistent file using STATUS='OLD'.
x05Values of status affect only disk files,and are ignored for devices such as the keyboard or the printer.
没有发现'readonly'这个选项.

cpu_time返回的是一个浮点型,表示当前CPU运行时间（以秒计）
示例：
program test_cpu_time
real ::start,finish
call cpu_time(start)
!put code to test here
call cpu_time(finish)
print '("Time = ",f6.3," seconds.")',finish-start
end program test_cpu_time
date_and_time稍微复杂,但包括很多信息,包括日期,时间,时区,最后一个参数是一个大小为8的整型数组,记录了年,月,日,时区差（以分钟计）,小时,分钟,秒,毫秒.
program test_time_and_date
character(8) ::date
character(10) ::time
character(5) ::zone
integer,dimension(8) ::values
!using keyword arguments
call date_and_time(date,time,zone,values)
call date_and_time(DATE=date,ZONE=zone)
call date_and_time(TIME=time)
call date_and_time(VALUES=values)
print '(a,2x,a,2x,a)',date,time,zone
print '(8i5))',values
end program test_time_and_date
system_clock用于决定处理器时钟周期.count_rate用于决定每秒CPU的时钟计数,
PROGRAM test_system_clock
INTEGER ::count,count_rate,count_max
CALL SYSTEM_CLOCK(count,count_rate,count_max)
WRITE(*,*) count,count_rate,count_max
END PROGRAM

错误原因如ls
至于“但当我把所有的i和x交换,其他不变,就没有问题,能够的到正确的结果.”是因为fortran的I-N规则,以字母IJKMLN开头的变量,如果没有定义,默认为整形,交换以后,虚参I是整形,所以就没有问题了

ec是用在write语句里的
格式为rec=r
表示当前访问文件中的第r个记录.

program randPoint
implicit none
integer i
real ::x,y
call random_seed () !随机种子
do i=1,100
call random_number (x)
call random_number (y)
write(*,"(a,f6.3,f6.3,a)") "(",x,y,")"
end do
end program random

integer,allocatable::rex(*),rey(*),rez(*),ref(*)
改为
integer,allocatable::rex(:),rey(:),rez(:),ref(:)
如果数组的大小出现在虚参中,可以不用 allocatable,直接 integer rex(I),前提是 I 是虚参.

不知道你这程序是干什么的
但感觉你里边的一个 is(K)=j 应该改成 js(K)=j
改了后可以运行

lkind 在前面定义了,肯定是一个 integer 的 parameter 常量.
它用来表示整数的 kind.
所谓 kind ,也就是整型的长度.具体值因不同编译器而不同.大多数主流编译器,如果 kind 是 2,表示短整型,如果是 4,则表示整型（32位）.如果 kind 是 8 则表示长整型.
如果写成
integer(2)
integer(4)
integer(8)
或许你会更熟悉它们一些.
之所以用 lkind 代替,是为了以后修改程序方便.试想,如果你要把整个程序的 integer(4) 改为 integer(8),是具有较大工作量的,而用 lkind代替,以后只需要修改 lkind 的值既可.
另外,如果你的代码更换为其他编译器,正如我上面所说,不同编译器,对 kind 的具体值规定不同,也需要对 lkind 进行修改.
L=int(PJ(1,I),Ikind) 也是处于同样的目的.int 函数接收第二个参数来确定转换的 kind
integer(Ikind),intent(out) ::JC(6) 这句话,我相信你能理解.定义一个kind为 lkind 的整型数组,名字叫 JC,大小是 6.
intent(OUT) 可能你理解有困难,他表示,JC 这个数组,出现在函数虚参里,而且作为该函数的一个输出参数.在该函数返回时,JC 必须赋值（因为它是一个输出参数）

F77的用法,一般用标号指向一个format（格式语句）.这里合并起来写了.

(a9,i3,2(a11,i5))表示的含义应该是：
输出宽度为9的实数、宽度为3的整数、然后输出2个(宽度为11的实数、宽度为5的整数）

从一个名为fort.1或已被打开且通道号为1的文本文件中读入6个数付给数组tout

NRND1 如果是函数,请书写它.
NRND1 如果是数组,请声明它.

eal::d,e,x1,x2 要放在所有可执行语句之前.
write(*,*)"a=" 引号要用单引号.

allocate(a(n),b(n),c(n),d(n),g(n),q(n),o(n))
这一句要放到read(*,*) n 之后
还有一个错误,fortran的数组是从1开始的,当你改正上面的错误了之后,接下来会抱错关于o(0)的问题,所以建议修改.

program ex01
implicit none
integer , parameter :: n=50 , l=1 , p=1
real , parameter :: h= 1.0* l / (n-1)
Real , parameter :: pi=3.1415926535
Real :: a(n-2,n),b(n-2,1),y(n,1),q(n,1)
integer :: i
y(:,:)=0
do i=1,n
q(i,1)=p*sin((i-1)*h*pi/l)
end do
do i=1,n
write(*,*) q(i,1)
end do
end program ex01h=l/(n-1) = 0
注意 l 和 n 是整数.
另外你的编程风格非常不好,建议这样写.

write(*,*) ' 请输入三角形三个顶点的坐标 x1,y1,x2,y2,x3,y3：'read(*,*) x1,y1,x2,y2,x3,y3a=sqrt((x2-x1)^2+(y2-y1)^2)b=sqrt((x3-x1)^2+(y3-y1)^2)c=sqrt((x3-x2)^2+(y3-y2)^2)L=(a+b+c)/2s=sqrt(L*(L-a)*(L-b)*...

你可能是按错工具栏的按钮了,要按那个感叹号的,或者快捷键 Ctrl+F5
你按的是调试的执行按钮,没有设置断点,所以直接正常退出了.

L4运算结果为FALSE

intent(IN） 表示这个参数是输入的
intent(OUT) 表示参数是输出的
intent(INOUT)表示这个参数同时用于两个方向的数据传递

我估计,由于你随便输入值赋给x,不会恰好是pi/2或者-pi/2的倍数,所以运算后,-1

主要问题：调用gaussj时,实参与虚参类型不一致.您在主程序中定义数组a和b是双精度8个字节的,而在子程序中却是默认的4个字节.我输入x(1:3)=(5.5 ,5.5, 5.5)时输出结果如下图：
SUBROUTINE gaussj(a,n,np,b,m,mp) !引入GaussJ消元法INTEGER m,mp,n,np,NMAXREAL*8 a(np,np),b(np,mp)PARAMETER (NMAX=50)INTEGER i,icol,irow,j,k,l,ll,indxc(NMAX),indxr(NMAX),ipiv(NMAX)REAL big,dum,pivinvdo 11 j=1,nipiv(j)=011continuedo 22 i=1,nbig=0.do 13 j=1,nif(ipiv(j).ne.1)thendo 12 k=1,nif (ipiv(k).eq.0) thenif (abs(a(j,k)).ge.big)then big=abs(a(j,k)) irow=j icol=kendifelse if (ipiv(k).gt.1) thenpause 'singular matrix in gaussj'endif12continueendif13continueipiv(icol)=ipiv(icol)+1if (irow.ne.icol) thendo 14 l=1,ndum=a(irow,l)a(irow,l)=a(icol,l)a(icol,l)=dum14continuedo 15 l=1,mdum=b(irow,l)b(irow,l)=b(icol,l)b(icol,l)=dum15continueendifindxr(i)=irowindxc(i)=icolif (a(icol,icol).eq.0.) pause 'singular matrix in gaussj'pivinv=1./a(icol,icol)a(icol,icol)=1.do 16 l=1,na(icol,l)=a(icol,l)*pivinv16continuedo 17 l=1,mb(icol,l)=b(icol,l)*pivinv17continuedo 21 ll=1,nif(ll.ne.icol)thendum=a(ll,icol)a(ll,icol)=0.do 18 l=1,na(ll,l)=a(ll,l)-a(icol,l)*dum18continuedo 19 l=1,mb(ll,l)=b(ll,l)-b(icol,l)*dum19continueendif21continue22continuedo 24 l=n,1,-1if(indxr(l).ne.indxc(l))thendo 23 k=1,ndum=a(k,indxr(l))a(k,indxr(l))=a(k,indxc(l))a(k,indxc(l))=dum23continueendif24continuereturn end
program mainimplicit none
integer,parameter :: row=3integer,parameter :: col=3integer nn!关于xinteger mm!关于binteger iiinteger jj integer kkreal*8:: HI(row,col) real*8 AI(row,col) !储存矩阵HI的值real*8 x(row)real*8 b(row)real*4 iii !整数型与浮点型的转换real*4 jjj
do ii=1,row,1!给H矩阵赋值do jj=1,col,1iii=iijjj=jjHI(ii,jj)=1/(iii+jjj-1)AI(ii,jj)=HI(ii,jj)!将初始矩阵H的值赋给替换矩阵Awrite(*,"(A3,I1,I1,A3,F9.4)")"H(",ii,jj,")=",HI(ii,jj)enddoenddo
write(*,*) "The substitude matrix is: " !显示替换矩阵Ado ii=1,row,1do jj=1,col,1write(*,"(A3,I1,I1,A3,F9.4)")"A(",ii,jj,")=",AI(ii,jj)enddoenddowrite(*,*) "Enter the x: " !给x赋值do nn=1,row,1read(*,*) x(nn)enddo
do mm=1,row,1!求b的值do jj=1,col,1b(mm)=b(mm)+HI(mm,jj)*x(jj)enddowrite(*,"(A3,I1,A3,F9.4)")"b(",mm,")=",b(mm)enddo!以上代码将矩阵H、解X、常数b均已确定
call gaussj(HI,row,row,b,row,row)write(*,*) ""write(*,"(A20)") "The result is : " !显示经gaussj消元法得到的矩阵Hdo ii=1,row,1do jj=1,col,1write(*,"(A3,I1,I1,A3,F9.4)")"H(",ii,jj,")=",HI(ii,jj)enddoenddo
do mm=1,row,1write(*,*)"b(",mm,")=",b(mm)enddo
do ii=1,row,1x(ii)=0enddo
write(*,*) "The solution vectors are: "!显示解do ii=1,row,1do kk=1,col,1x(ii)=x(ii)+AI(ii,kk)*b(kk)enddowrite(*,*)"x(",ii,")=",x(ii)enddo stop end

写2就是跳到括号里第二个标号,220

补充：l 微分方程求解问题意大利生物学家Ancona曾致力于鱼类种群相互制约关系的研究,他从第一次世界大战期间,地中海各港口捕获的几种鱼类捕获量百分比的资料中,发现鲨鱼的比例有明显增加.为了反映在没有人工捕获的自然环境中食饵与捕食者之间的制约关系,建立一个食饵—捕食者系统的数学模型如下：①x1(t),x2(t)分别是食饵、捕食者（鲨鱼）在t时刻的数量； ②r1,r2是食饵、捕食者的固有增长率；③λ1是捕食者掠取食饵的能力,λ2是食饵对捕食者的供养能力；给定一组具体数据求解.

有一个不认识的字符 # 出现在不应该出现的位置.

应该这样program exam real v open(2,file='6.txt',status='unknown') v=sqrt(4.0)write(2,*) vclose(2)call system('c:windowsnotepad.exe 6.txt')end既然V是realsqrt(4.0)中间应该是4.0

代码：program data_sum
implicit none
integer :: num, sum = 0

do
read (*,*) num
if (num .eq. -1) then
exit
end if
if (num .gt. 0) then
sum = sum + num
end if
end do
write (*,'(a, i6)') 'sum = ', sum

end program data_sum运行结果：1
2
3
4
5
6
-1
sum =21

我挑了其中一些有代表性的语句进行了注释,惊叹号(!)是Fortran里的注释.
很多语句可以举一反三.
integer row,i,j,k,n
!声明了五个整型变量
real s
!声明实型变量
real,dimension(:,:),allocatable::a,l,u
!声明实型二维数组
real,dimension(:),allocatable::b,x,y
!声明实型一维数组
read*,row
!从键盘读数据row
allocate(a(row,row))
!给二维数组a分配内存空间
read*,((a(i,j),j=1,row),i=1,row)
!从键盘读数组a
allocate(b(row))
read*,(b(i),i=1,row)
allocate(y(row))
allocate(x(row))
do k=1,row-1
!循环语句,k从1循环到row-1
do i=k+1,row
a(i,k)=a(i,k)/a(k,k)
do j=k+1,row
a(i,j)=a(i,j)-a(i,k)*a(k,j)
end do
!变量j的循环体结束
end do
end do
allocate(l(row,row))
do i=1,row
do j=1,row
if(i>j)then
!判断语句,i是否大于j
l(i,j)=a(i,j)
elseif(i==j)then
l(i,j)=1
else if(ij)then
u(i,j)=0
end if
end do
end do
print*,'a=',((a(i,j),j=1,row),i=1,row)
!在屏幕上输出数组a
print*,'l=',((l(i,j),j=1,row),i=1,row)
print*,'u=',((u(i,j),j=1,row),I=1,row)
n=10
do while(n>0)
!满足n>0时循环
y(1)=b(1)/l(1,1)
do i=2,row
sum=0
do j=1,i-1
sum=sum+l(i,j)*y(j)
end do
y(i)=(b(i)-sum)/l(i,i)
end do
print*,(y(i),i=1,row)
x(row)=y(row)/u(row,row)
do i=row-1,1,-1
!i从row-1循环到1,步长是-1,也就是倒着循环
sum1=o
do j=i+1,row
sum1=sum1+u(i,j)*x(j)
end do
x(i)=(y(i)-sum1)/u(i,i)
end do
print*,(x(i),i=1,row)
s=abs(x(1))
!abs是绝对值函数
do j=2,row
if(abs(x(j))>s)then
s=abs(x(j))
end if
end do
do i=1,row
b(i)=x(i)/s
end do
n=n-1
end do
end
!最后的end代表程序结束
!这里应该是开始了一个子程序
INTEGER ROW,COL,I,J,N
REAL,DIMENSION(:,:),ALLOCATABLE::B
REAL,DIMENSION(:),ALLOCATABLE::X,G
READ*,row,col
allocate(b(row,col))
read*,((b(i,j),j=1,col),i=1,row)
allocate(x(col),g(col))
read*,(x(i),i=1,col)
read*,(g(i),i=1,col)
read*,n
do while(n.ne.0)
!.ne.是not equal的意思,也就是判断n是否不等于0
do i=1,col
sum=0
do j=1,col
sum=sum+b(i,j)*x(j)
end do
x(i)=sum+g(i)
end do
n=n-1
print*,(x(i),i=1,col)
end do
end

正确程序如下：----------------------------
!主程序
program main
external f1,f2
double precision x,y,E,a,b,f1,f2,Ex,Ey
open(9,file="v(x,y).dat",status='unknown')
a=3.14/3.0
b=2.0*3.14/3.0
do i=1,n
x=0.0-2.0+i*0.05
if(x

如果 str 字符串的第 i 个字母等于 0,则退出（循环）
ichar 是取字符的 ASCII 码.为 0 的字符与 C 语言的 类似.

看不到 我小学生 看楼下怎么解释 我也学习一下

建议认真学习 Fortran 语法.
作业是你提高的很好方式,也是锻炼自己的重要途径.
如果你连作业题都求助,让别人给你现成的,你永远不会进步.

如果 block 等于 'gene',调用 rgene 函数

NINT是四舍五入取整
在python中用round直接完成.
算法上其实是加0.5后再截尾取整
3.4 +0.5得3.9 ...截尾得3
3.6 +0.5得4.1 ...截尾得4

10 20 30 是标号,分别是：
10代表 Y=X*X+1
20 代表Y=SIN(X+1)
30 代表Y=SQRT(X+2)
IF(X)10,20,30的意思是判断x的值三种情况,小于0则转到10代表的那一行,等于0则转到20,大于0则转到30那一行

你给出的 random 没有任何虚参,不存在类型不一致的情况.
1.随机数,Fortran 有内置的函数random_number可以实现.
2.如果你要找出你的问题,建议你给出完整代码.

这不是代码,而是代码的解释.[] 中的部分为可写,可不写.
意思是,可以写为：
Do s i = e1 ,e2
或
Do i = e1 ,e2
或
Do s i = e1 ,e2 ,e3
或
Do i = e1 ,e2 ,e3
这句是循环语句,意思是让 i 从 e1 循环到 e2,每次增加 e3（步长）.如果 e3 没有写,默认为 1
s 的意思是指定结束语句.如果没有指定结束语句,则以 End Do 结束

这是一个输出语句,输出目标是 6 号文件.自由格式.
在某些编译器下,6号文件默认是显示器.
输出内容为：Speed of sound = 和 rho 的值.
假设 rho 为 10,则输出 Speed of sound = 10
这些都是Fortran最基本的赋值和函数调用语句.
我相信你的代码远远不止这一句,而这一句,也不是你全部的障碍.其他语句你肯定也看得吃力.
如果你连这个还看不明白,那么你需要的是系统的学习 Fortran.
这样断章取义的抽出一两个语句来求助,对你没有太大的帮助.

比如 a%b 表示结构体变量 a 的 b 成员
类似于 C 语言的 a.b,结构体（struct）
inpKey%nKeys
inpKey 是一个结构体,你可以搜索一下 Type inpKey 找到它的定义
nKeys 是它的一个结构体成员.就这样

print*,mod(x,3),mod(x,5)

呵呵,智能天线啊.
不过这个没给参数别人没法给你编啊,你可以找论文作者要程序,前提是人家愿意给你.
波束成型其实就是些向量相乘,程序不会太难,关键是看懂算法,搞清楚物理意义啊.

举个最简单的例子,Jacobi迭代解x+y=2,2x+3y=5
program Jacobi
integer::N,NMAX=10000
real::X,Y,X0,Y0,ERROR=0.0001
X0=0.0
Y0=0.0
do N=1,NMAX
X=2-Y0
Y=(5.0-2*X0)/3.0
if (abs((X-X0)/X0)

我想你的问题,可能是符号计算的问题.
一般学Fortran刚开始都是做数值计算,而不是做符号计算.数值计算在计算机实现里更容易.
如果你想做符号计算,可以用matlab,Mathematica 这类的数学工具来进行.
当然,Fortran是一门编程工具,理论上matlab能实现的他都能,他也能实现符号计算,但是恐怕需要你自己去书写大量的代码来实现.而不由语法里直接提供

我从来不用format语句,如果要求不高的话可以
write(*,*) ((a(i,j),j=1,n),i=1,n)
或者
do i=1,6
write(*,*) a(i,:)
enddo
如果你对格式要求很高,那么最不济的办法可以用循环并分实部和虚部来输出,这样就全是实数了.

#include
void main()
{
int XY[102][102];
int step[4][2] = {{1,0},{0,1},{-1,0},{0,-1}};//走动的方向
int x,y;
for(y=0;y