Документ взят из кэша поисковой машины. Адрес оригинального документа : http://foroff.phys.msu.ru/illposed/programs/Fig4_7.txt Дата изменения: Mon Jul 7 00:04:43 2008 Дата индексирования: Mon Oct 1 22:04:57 2012 Кодировка:

Fig 4.7
=======

00001 IMPLICIT REAL*8(A-H,O-Z)
00002 IMPLICIT INTEGER*4(I-N)
00003 DIMENSION U0(41),Z(41),U(41,260),Z0(41)
00004 EXTERNAL AK
00005 X1=0.
00006 X2=1.
00007 Y1=-2.
00008 Y2=2.
00009 N=41
00010 M=41
00011 IMAX=400
00012 IC=2
00013 DL=0.
00014 C1=0.
00015 C2=1.
00016 DO 1 J=1,N
00017 X=X1+(X2-X1)/(N-1.)*(J-1.)
00018 1 Z0(J)=1.-X*X
00019 CALL PTICR0(AK,U,X1,X2,Y1,Y2,N,M)
00020 CALL PTICR3(U,Z0,U0,N,M)
00021 DO 4 J=1,N
00022 Z(J)=0.
00023 4 CONTINUE
00024 CALL PTIGR(AK,U0,X1,X2,Y1,Y2,N,M,Z,AN,
*ITER,DL,IMAX,C1,C2,IC,U,41*260,IERR)
00025 PRINT 5,Z0,Z,DL,AN,ITER,IERR
00026 STOP
00027 5 FORMAT(15X,'Exact solution:'/'.'/
* 8(5F11.6/),F11.6/'.'/
*15X,'Approximate solution:'/'.'/
* 8(5F11.6/),F11.6/'.'/
*10X,'Right side error square :',
* D14.6/
*10X,'Discrepancy :',
* D14.6/
*10X,'Number of iterations :',
* I14/
*10X,'Return code (IERR) :',
* I14)
00028 END


00001 FUNCTION AK(X,Y)
00002 IMPLICIT REAL*8(A-H,O-Z)
00003 AK=1./(1.+100.*(X-Y)**2)
00004 RETURN
00005 END





Exact solution:

1.000000 .999375 .997500 .994375 .990000
.984375 .977500 .969375 .960000 .949375
.937500 .924375 .910000 .894375 .877500
.859375 .840000 .819375 .797500 .774375
.750000 .724375 .697500 .669375 .640000
.609375 .577500 .544375 .510000 .474375
.437500 .399375 .360000 .319375 .277500
.234375 .190000 .144375 .097500 .049375
.000000

Approximate solution:

.988767 .988767 .988576 .987790 .986477
.984712 .981902 .974724 .963316 .951637
.939958 .928279 .916600 .896694 .876788
.856882 .836975 .814705 .791996 .769288
.746579 .723870 .701114 .676169 .647221
.616290 .581252 .545292 .508441 .469801
.430246 .390691 .351135 .311580 .272024
.232469 .192914 .150570 .108227 .065883
.023540

Right side error square : .000000D+00
Discrepancy : .207510D-06
Number of iterations : 400
Return code (IERR) : 0