Resitev kvadratne enacbe
Kvadratno enacbo lahko rešimo s pomocjo kompleksnih števil

program kvadrat
implicit none
complex :: x1,x2
real :: a,b,c,d
write (*,*) "vnesi parametre kvadratne enacbe"
read (*,*) a, b, c
d=b**2-4e0*a*c
x1=cmplx(-b/2e0/a)+csqrt(cmplx(d))/2e0/a
x2=cmplx(-b/2e0/a)-csqrt(cmplx(d))/2e0/a
write (*,*)'x1=',real(x1),'+i(',imag(x1),')'
write (*,*)'x2=',real(x2),'+i(',imag(x2),')'
end program kvadrat

Prikaz



Nazaj na osnovno stran za numercne metode.