| 查看: 580 | 回复: 3 | ||
1312758691新虫 (小有名气)
|
[求助]
输出格式 已有1人参与
|
|
程序为： PROGRAM least_squares_fit IMPLICIT NONE ! 数据字典，声明常量 INTEGER,PARAMETER :: LU=18 ! I/O unit for disk I/O ! 数据字典，声明变量类型、定义和单位 CHARACTER(len=24) :: filename INTEGER :: ierror INTEGER :: n=0 REAL :: slope REAL :: sum_x=0. REAL :: sum_x2=0. REAL :: sum_xy=0. REAL :: sum_y=0. REAL :: sum_y2=0. REAL :: x REAL :: x_bar REAL :: y REAL :: y_bar REAL :: y_int ! 提示用户输入并获得输入文件名 WRITE (*,1000) 1000 FORMAT (1x,'This program performs a least-squares fit of an ',/,& 1x,'input data set to a straitght line. Enter the name',/& 1x,'of the file containing the input (x,y) pairs: ') READ (*,1010) filename 1010 FORMAT (A) ! 打开输入文件 OPEN (UNIT=LU,FILE=filename,STATUS='OLD',IOSTAT=ierror ) ! 查看OPEN操作是否失败 errorcheck: IF (ierror>0) THEN WRITE (*,1020) filename 1020 FORMAT(1x,'ERROR: File ',A,'does not exist! ') ELSE !如果文件成功打开，就从文件中读取数据对(x,y)。 DO READ (LU,*,IOSTAT=ierror) x,y !获得数据对 IF (ierror/=0) EXIT n=n+1 sum_x=sum_x+x sum_y=sum_y+y sum_x2=sum_x2+x**2 sum_xy=sum_xy+x*y END DO ! 现在计算斜率和截距。 x_bar=sum_x/real(n) y_bar=sum_y/real(n) slope=(sum_xy-sum_x*y_bar)/(sum_x2-sum_x*x_bar) y_int=y_bar-slope*x_bar ! 告知用户 WRITE (*,1030) slope,y_int,N 1030 FORMAT ('0','Regression coefficients for the least-squares line:',& /,1x,'slope(m)=',F12.3,& /,1x,'Intercept(b)=',F12.3,& /,1X,'No of points=',I12) ! 关闭输入文件并退出。 CLOSE (UNIT=LU) END IF errorcheck END PROGRAM least_squares_fit 运行结果为图片所示，为什么运行后图片箭头所指的“0”还在？这跟“1030 FORMAT ('0','Regression coefficients for the least-squares line:',& /,1x,'slope(m)=',F12.3,& /,1x,'Intercept(b)=',F12.3,& /,1X,'No of points=',I12)”这个语句的意思不符合啊！ 不懂，赐教！我用的CVF6.6的编译器！  1.jpg |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
baobiao007
木虫 (职业作家)
中国特色
- 应助: 201 (大学生)
- 金币: 6482.7
- 散金: 557
- 红花: 40
- 帖子: 3050
- 在线: 1009.9小时
- 虫号: 505962
- 注册: 2008-02-18
- 专业: 应用地球物理学
2楼2014-05-07 16:43:27
1312758691
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1334
- 散金: 12
- 帖子: 86
- 在线: 60.4小时
- 虫号: 2328063
- 注册: 2013-03-07
- 专业: 理论和计算化学
3楼2014-05-07 19:26:39
baobiao007
木虫 (职业作家)
中国特色
- 应助: 201 (大学生)
- 金币: 6482.7
- 散金: 557
- 红花: 40
- 帖子: 3050
- 在线: 1009.9小时
- 虫号: 505962
- 注册: 2008-02-18
- 专业: 应用地球物理学
4楼2014-05-07 19:32:39
20