Advertisement
Guest User

fdm

a guest
Apr 29th, 2018
68
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. module fdm
  2.     implicit none
  3.  
  4. contains
  5.     subroutine generate(A, X, N, beginc, endc, h)
  6.         implicit none
  7.         integer (kind = 4), intent(in) :: N, beginc, endc
  8.         real (kind = 8), intent(in) :: h
  9.         real (kind = 8), allocatable, dimension(:, :), intent(inout) :: A
  10.         real (kind = 8), allocatable, dimension(:), intent(inout) :: X
  11.         integer :: i
  12.  
  13.         A = 0
  14.         X = 0
  15.        
  16.         do i=2, N
  17.             A(i, i-1) = 1/(h**2)
  18.         end do
  19.  
  20.         do i=1, N
  21.             A(i, i) = -2/(h**2)
  22.         end do
  23.  
  24.         do i=1, N-1
  25.             A(i, i+1) = 1/(h**2)
  26.         end do
  27.  
  28.         X(1) = beginc
  29.         X(N) = endc
  30.  
  31.     end subroutine
  32. end module fdm
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement