vvt 发表于 2018-1-23 22:59
tan 的定义域是以弧度制为单位，而不是值域。
所以，tan(lat) 要求 lat 是以弧度制为单位。其结果（值域 ...

好的，明白了

向着光亮 发表于 2018-1-23 21:18
tan计算的不就是一个数值吗？为何是弧度制？

tan 的定义域是以弧度制为单位，而不是值域。
所以，tan(lat) 要求 lat 是以弧度制为单位。其结果（值域），是个比值，无单位。

向着光亮 发表于 2018-1-23 21:18
tan计算的不就是一个数值吗？为何是弧度制？

请先回顾计算方法课程的相关内容：三角函数值
传送门：
http://blog.csdn.net/eric88/article/details/17031311

fcode 发表于 2018-1-23 12:58
tan 就是弧度制的，不是你想什么制，就什么制。

tan计算的不就是一个数值吗？为何是弧度制？

向着光亮 发表于 2018-1-23 12:40
我刚才又算了一遍，还是不行，大体思路我明白了，主要是角度和弧度的互换，太容易出错了，你刚才说的纬度 ...

坚持角度制的话，请查找手册，譬如：
https://gcc.gnu.org/onlinedocs/g ... an/ASIND.html#ASIND

9.20 ASIND — Arcsine function, degrees
Description:
ASIND(X) computes the arcsine of its X in degrees (inverse of SIND(X)).

This function is for compatibility only and should be avoided in favor of standard constructs wherever possible.

Standard:
GNU Extension, enabled with -fdec-math.

Class:
Elemental function

Syntax:
RESULT = ASIND(X)

Arguments:
X        The type shall be either REAL and a magnitude that is less than or equal to one - or be COMPLEX.
Return value:
The return value is of the same type and kind as X. The real part of the result is in degrees and lies in the range -90 \leq \Re \asin(x) \leq 90.

Example:
program test_asind
  real(8) :: x = 0.866_8
  x = asind(x)
end program test_asind
Specific names:
Name        Argument        Return type        Standard
ASIND(X)        REAL(4) X        REAL(4)        GNU Extension
DASIND(X)        REAL(8) X        REAL(8)        GNU Extension
See also:
Inverse function: SIND Radians function: ASIN

tan 就是弧度制的，不是你想什么制，就什么制。

pasuka 发表于 2018-1-23 11:21
依旧是问题代码：
1、明明打开同一个文件，open语句为何放在do循环内？
2、有open无close；

你好，听了你的建议我觉得确实应该严谨的 写才能避免出错，我还想问问你这个计算的时候所有的数字都得写的浮点型的吗，定义变量也要这样吗

pasuka 发表于 2018-1-23 11:21
依旧是问题代码：
1、明明打开同一个文件，open语句为何放在do循环内？
2、有open无close；

你说的对，这样也能算出来，把open放出来和在里面都是可以的，但是习惯很重要，所以我决定拿出来了，，至于这个单双精度浮点问题刚开始接受，需要慢慢适应