Предыдущая страница реферата | 1  2  3  4  5  6  7  8  9 | Следующая страница реферата

+ y(5,i)*Q(y(1,i)+1,y(2,i))+

+ y(6,i)*Q(y(1,i),y(2,i)+1)

c------------------------------------------------------------------ c вычисление погрешности на данной итерации c------------------------------------------------------------------ if (abs(NEXTQ-q(y(1,i),y(2,i))).gt.EI)

+ EI=abs(NEXTQ-q(y(1,i),y(2,i))) c print *,'x=',y(1,i),' y=',y(2,i) q(y(1,i),y(2,i))=NEXTQ enddo c print '(16h Итерация номер ,i5,13h погрешность=,E15.7)',S,EI if (EI.gt.E)goto 1

c do i=P,1,-1 c print '(21e10.3)',(q(j,i),j=1,M) c enddo end

c------------------------------------------------------------------ c ПОДПРОГРАММА АЛФАВИТНО-ЦИФРОВОГО,МОЗАИЧНОГО c ВЫВОДА РЕЗУЛЬТАТА c integer M-число узлов по оси X c integer P-число узлов по оси Y c real Q(M,P)-входной массив значений Y c c c------------------------------------------------------------------ subroutine outdata(M,P,q) character a(11)/'.','+','*','','','-','-','-','','-','-'/ integer M,P,I,J real q(M,P) do J=P,1,-1 print '(400A2)',(a(int(q(I,J)/21)+1),I=1,M),

+ (a(int(q(I,J)/21)+1),I=M-1,1,-1) enddo do I=1,10 print *,'''',a(I),'''','---> от ',20*(I-1),', до ',

+ 20*I,'(включительно)' enddo end c------------------------------------------------------------------ c ПОДПРОГРАММА ВЫЧИСЛЕНИЯ ОШИБКИ c real q-массив значений Y с шагом =2*h c real qq-массив значений Y с шагом =h c real E-значение погрешности c c------------------------------------------------------------------ subroutine mistake(M,P,q,qq,E) integer M,P,iq,jq,iqq,jqq real qq(M,P),q(int(M/2)+1,int(P/2)+1),max,E,other max=0 iq=0 do iqq=1,P,2 iq=iq+1 jq=0 do jqq=1,M,2 jq=jq+1 other=abs(q(jq,iq)-qq(jqq,iqq)) if (other.gt.max)max=other enddo enddo print *,M,' ',P,' ',max/3 if (max/3.lt.E) then call outdata(M,P,qq)

Stop endif end c------------------------------------------------------------------ c ОСНОВНАЯ ПРОГРАММА c c c------------------------------------------------------------------ integer N/90000/,M,P,flag/0/ real y(6,90000),q(300,300),H/.05/,K/.05/,E/.5/,qq(300,300) real EZ/.01/ c print *,'Введите шаг вдоль оси X ' c read (*,*)H c print *,'Введите шаг вдоль оси Y ' c read (*,*)K c print *,'Введите точность вычислений ' c read (*,*)E

M=.2/H+1

P=.4/K+1 call mkr(H,K,N,y,M,P,q) call zeidel(N,y,M,P,q,EZ)
111 H=H/2

K=K/2

M=.2/H+1

P=.4/K+1

N=90000 if (flag.eq.0)then flag=1 call mkr(H,K,N,y,M,P,qq) call zeidel(N,y,M,P,qq,EZ) call mistake(M,P,q,qq,E) else flag=0 call mkr(H,K,N,y,M,P,q) call zeidel(N,y,M,P,q,EZ) call mistake(M,P,qq,q,E) endif goto 111 end

Литература.

1. И.С.Березин,Н.П.Жидков ’Методы вычислений’,том 1,М.,1966,632 стр.
2.’ Численные методы решения задач на ЭВМ ’ , Учебное пособие ,
Г.Н.Рубальченко , К. , 1989 , 148 стр.
3.’Справочник языка ФОРТРАН’ , М.,1996 ,106 стр.
--------------------
Т=100 0 С

Т=200 0 С

Т=0 0 С

Т=0 0 С

0,1

0,4

Выбор начального шага

Составление системы ур-ний

№1

Рекомендуем скачать другие рефераты по теме: реферат техника, решебник 6 класс виленкин.



Предыдущая страница реферата | 1  2  3  4  5  6  7  8  9 | Следующая страница реферата

Поделитесь этой записью или добавьте в закладки