fortran中的立方根怎么算？

冬季极端降水

          急急急！大家快来可怜可怜我呀，我算x的立方根，按理说，我写的是x**(1.0/3.0)，可是结果就是出不来呀{:soso_e135:}，
快把我急坏了，出现的错误就是这样 ，求高手呀！！！

cnphoenix

用 ( 1.0/3.0 ) 表示1/3

言深深

楼上有理1/3是整型运算，直接变成0，1.0/3就可以了。
不过楼主的问题，彷佛不在这儿，你的问题在于开根号的里面出现不大于0的数。

mofangbao

搂住能输出这句话每个变量的数值出来看下么，一次的即可

残月

编程警示：永远不要试图对一个负数进行实数幂运算。

楼主的问题我个人认为：先加一个判断，如果底数大于零，就用 exp（1.0/3*log（real(。。。)））来代替原来的语句;如果底数小于零，则在前面加负号，即 -exp（1.0/3*log（abs（real(。。。)）））
’。。。‘为你括号里的东西

你自己试试，希望起作用

mofangbao

残月 发表于 2012-5-9 09:28

                               
登录/注册后可看大图

编程警示：永远不要试图对一个负数进行实数幂运算。

楼主的问题我个人认为：先加一个判断，如果底数大于 ...

抛开他的环境不说，判断加负号还不如直接用abs求绝对值呢

残月

mofangbao 发表于 2012-5-9 09:37

                               
登录/注册后可看大图

抛开他的环境不说，判断加负号还不如直接用abs求绝对值呢

直接加绝对值会影响结果吧，负数开立方是负的啊？这个怎么解决？请教啦

mofangbao

残月 发表于 2012-5-9 10:22

                               
登录/注册后可看大图

直接加绝对值会影响结果吧，负数开立方是负的啊？这个怎么解决？请教啦

额。。。我错了，哈哈

虫儿飞

两种方法：一、可以试试将（1.0/3.0）,直接取为0.333，你可以自己调整保留位数，只要精度在可以接受的范围内。
二、麻烦一点，需要专门编个子程序，但是也比较简单。假设一个数x，x的立方是要开根的数，x即为要求的数，设定x的步长为0.001，或者其他，根据你对结果精度的要求，立方之后，求差，差在一个你设定的范围内，此时x即为结果。

言深深

虫儿飞 发表于 2012-5-9 17:34

                               
登录/注册后可看大图

两种方法：一、可以试试将（1.0/3.0）,直接取为0.333，你可以自己调整保留位数，只要精度在可以接受的范围内 ...

哈哈，你的方法是数值方法求解方程的吧？这样的话如果保持高精度，计算机过程可能比较长