Высокопроизводительные вычисления на кластерах. Беликов Д.А - 139 стр.

UptoLike

ВУЗ: 

ТГУ | Томск

Составители: 

Рубрика: 

139
common/parallel/ rank, size, m
external rp0
! процедура вычисления правых частей уравнений системы
! ОДУ
c Занесение начальных условий в массив yy
do i=1,n
yy(0,i)=y(i)
end do
c
k=1
!k – счетчик номера шага метода для занесения решения в
!массив yy
do while (x.lt.3*h)
call rp0(x, y, k1)
do i=1,n
y1(i)=y(i)+k1(i)*h/2
end do
call rp0(x+h/2, y1, k2)
do i=1,n
y1(i)=y(i)+k2(i)*h/2
end do
call rp0(x+h/2, y1, k3)
do i=1,n
y1(i)=y(i)+k3(i)*h
end do
call rp0(x+h, y1, k4)
do i=1,n
y(i)=y(i)+h*(k1(i)+2*k2(i)+2*k3(i)+k4(i))/6
end do
x=x+h
do i=1,n
yy(k,i)=y(i)
end do
k=k+1
end do
return
end
В программе решается задача

Страницы