Перейти к содержанию

Решение систем гиперболических уравнений/Разностные схемы/Реализация на Fortran-e

Материал из Викиучебника — открытых книг для открытого мира
     subroutine makestep_1d(U,U1,lambda,tau,h,MK)
!
!      Подпрограмма вычисляет значение U на n+1-м слое по времени и 
!      сохраняет значение в массиве U1
!      Использована простейшая схема "левый уголок".
!
       implicit none
       integer :: MK
       real*8,dimension(MK) :: U,U1
       real*8 :: lambda,tau,h
       
       integer :: m
       do m=2,MK
         U1(m)=U(m)-(tau*lambda/h)*(U(m)-U(m-1))
       end do
     end subroutine