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

Fig 4.11
========



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=3
00013 DL=0.
00014 DO 1 J=1,N
00015 X=X1+(X2-X1)/(N-1.)*(J-1.)
00016 1 Z0(J)=4.*X*(1.-X)
00017 CALL PTICR0(AK,U,X1,X2,Y1,Y2,N,M)
00018 CALL PTICR3(U,Z0,U0,N,M)
00019 DO 4 J=1,N
00020 Z(J)=0.
00021 4 CONTINUE
00022 CALL PTISR(AK,U0,X1,X2,Y1,Y2,N,M,Z,AN,
*ITER,DL,IMAX,IC,U,41*260,IERR)
00023 PRINT 5,Z0,Z,DL,AN,ITER,IERR
00024 STOP
00025 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)
00026 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:

.000000 .097500 .190000 .277500 .360000
.437500 .510000 .577500 .640000 .697500
.750000 .797500 .840000 .877500 .910000
.937500 .960000 .977500 .990000 .997500
1.000000 .997500 .990000 .977500 .960000
.937500 .910000 .877500 .840000 .797500
.750000 .697500 .640000 .577500 .510000
.437500 .360000 .277500 .190000 .097500
.000000

Approximate solution:

.000000 .097502 .189998 .277498 .360000
.437501 .510001 .577501 .639999 .697499
.749999 .797500 .840000 .877500 .910000
.937500 .960000 .977500 .990000 .997500
1.000000 .997500 .990000 .977500 .960000
.937500 .910000 .877500 .840000 .797500
.749999 .697499 .639999 .577501 .510001
.437501 .360000 .277498 .189998 .097502
.000000

Right side error square : .000000D+00
Discrepancy : .946349D-18
Number of iterations : 90
Return code (IERR) : 0