Cubic function problem

Help
1

while writing the program, I ran into an error that I can’t fix, the problem is on lines 35 and 63. I would ask for help.

subroutine cubic(a,b,c,d,x)
 
 real,INTENT(IN) :: a,b,c,d
 complex,INTENT(INOUT) :: x(*)
 real,parameter:: pi = 4.0d0*atan(1.0d0)
 real :: p,q,DD,phi,temp1,y1,y2,y3,temp2,u,v,y2r,y2i
 
 p = (3.0*c/a-(b/a)**2.0)/3.0
 q = (2.0*(b/a)**3.0-9.0*b*c/a/a + 27.0*d/a)/27.0
 
 DD = (p/3.0)**3.0 + (q/2.0)**2.0
 
 if(DD < 0.0)then

  phi=acos(-q/2.0/sqrt(abs(p)**3/27.0))
  temp1=2.0*sqrt(abs(p)/3.0)
  y1=temp1*cos(phi/3.0)
  y2=-temp1*cos((phi+pi)/3.0)
  y3=-temp1*cos((phi-pi)/3.0)
 else

  temp1=-q/2.0+sqrt(DD)
  temp2=-q/2.0-sqrt(DD)
  u=abs(temp1)**(1.0/3.0)
  v=abs(temp2)**(1.0/3.0)
  if(temp1 < 0.0) u=-u
  if(temp2 < 0.0) v=-v
  y1=u+v
  y2r=-(u+v)/2.0
  y2i=(u-v)*sqrt(3.0)/2.0
 end if

 temp1=b/a/3.0
 y1=y1-temp1
 y2=y2-temp1
 y3=y3-temp1
 y2r=y2r-temp1
 
 if(DD < 0.0) then
  x(1)=cmplx(y1,0.0)
  x(2)=cmplx(y2,0.0)
  x(3)=cmplx(y3,0.0)
 else if(DD > 0.0) then
  x(1)=cmplx(y1,0.0)
  x(2)=cmplx(y2r,y2i)
  x(3)=cmplx(y2r,-y2i)
 else
  x(1)=cmplx(y1,0.0)
  x(2)=cmplx(y2r,0.0)
  x(3)=cmplx(y2r,0.0)
 endif
end subroutine cubic

program cubicfunction
 implicit none

 real :: a,b,c,d
 complex :: x(1:3)

 print*,"Cubic function: ax^3 + bx^2 + cx + d = 0"
 read(*,*) a, b,c, d 

  call cubic(a,b,c,d,x)
  print *, "x1 =",x(1)
  print *, "x2 =",x(2)
  print *, "x3 =",x(3)

end program cubicfunction

when i try to activate the program i get this error
Zrzut ekranu 2022-05-07 202621

5

How should i set this in dd > 0?

7

i will try and thank you

8

If you calculate pi in double precision you should store it and all your other real and complex variables that way, e.g. with

integer,parameter:: dp = kind(1d0)

and then turn all the real and complex declarations into real(dp) and complex(dp).

1 Like