| 查看: 661 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有12人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
1312758691
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1334
- 散金: 12
- 帖子: 86
- 在线: 60.4小时
- 虫号: 2328063
- 注册: 2013-03-07
- 专业: 理论和计算化学
3楼2014-05-07 19:26:39
baobiao007
木虫 (职业作家)
中国特色
- 应助: 201 (大学生)
- 金币: 6484.2
- 散金: 557
- 红花: 40
- 帖子: 3050
- 在线: 1009.9小时
- 虫号: 505962
- 注册: 2008-02-18
- 专业: 应用地球物理学
2楼2014-05-07 16:43:27
baobiao007
木虫 (职业作家)
中国特色
- 应助: 201 (大学生)
- 金币: 6484.2
- 散金: 557
- 红花: 40
- 帖子: 3050
- 在线: 1009.9小时
- 虫号: 505962
- 注册: 2008-02-18
- 专业: 应用地球物理学
4楼2014-05-07 19:32:39
