Решение систем гиперболических уравнений/Разностные схемы/Реализация на 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