Система линейных уравнений в Фортране с использованием DGESV

Question

Я новичок в Fortran и пытаюсь решить ряд линейных уравнений с использованием DGESV. Я пытаюсь сделать простейшую программу возможной только для запуска метода, но получаю это сообщение об ошибке:

Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

Backtrace for this error:
#0  0x106cef092
#1  0x106cee3b0
#2  0x7fff8dda9f19
Segmentation fault: 11

Вот код:

   program linearEqSolver

   implicit none
   EXTERNAL         DGESV

   integer, parameter :: N = 3, NRHS = 1, LDA = N,LDB = N
   integer :: INFO
   real, dimension (N) :: IPIV, b
   real, dimension (N,N) :: A


   b = (/ 4, -6, 7/)
   A(:,1) = (/ 1, 1, -1/)
   A(:,2) = (/ 1, -2, 3/)
   A(:,3) = (/ 2, 3, 1/)

   CALL DGESV( N, NRHS, A, LDA, IPIV, B, LDB, INFO )

   end program