为什么我一直输0？

Question

我正在生成一个计算勒让德多项式的模块，但是一旦达到'P4（x）'，我开始得到0作为输出。

例如'P4（x）= 1/8 *（35x ^ 4 - 30x ^ 2 +3）'

所以我决定制作以下程序以检查任何可能的错误。

real (kind = 8) :: x,y
integer :: i

y = 0.0, x = -1.0

do i = 1, 30

y = ((1/8)*((35*(x**4))-(30*(x**2))+3)

write (6,*) y

cont = cont + 0.066

end do

此程序的输出始终相同0.000000

令人惊讶的是，它显示了符号的变化，但仍显示-0.000000。