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

Fig 4.6
=======



00001 IMPLICIT REAL*8(A-H,O-Z)
00002 IMPLICIT INTEGER*4(I-N)
00003 REAL*8 KERN
00004 DIMENSION U(6500),U0(32,32),Z0(32,32),
*Z(32,32),ALIM(12),KERN(32,32)
00005 EXTERNAL AK
00006 ALIM(5)=0.
00007 ALIM(6)=2.0
00008 ALIM(7)=0.
00009 ALIM(8)=2.0
00010 ALIM(9)=0.
00011 ALIM(10)=1.0
00012 ALIM(11)=0.0
00013 ALIM(12)=1.0
00014 ALPHA=1.000
00015 C1=1.000001D0
00016 N1=32
00017 N2=32
00018 IMAX=500
00019 R=0.
00020 ST1=(ALIM(6)-ALIM(5))/N1
00021 ST2=(ALIM(8)-ALIM(7))/N2
00022 DO 5 I=1,N1
00023 DO 5 J=1,N2
00024 S1=(I-0.5)*ST1-0.5
00025 S2=(J-0.5)*ST2-0.5
00026 Z0(I,J)=0.
00027 IF(S1.LT.0.0.OR.S1.GT.1.0.OR.S2.LT.0.0.OR.S2.GT.1.0)
* GOTO 5
00028 Z0(I,J)=((DEXP(-(S1-0.3)**2/0.03)+
+DEXP(-(S1-0.7)**2/0.03))/0.9550408-
-0.052130913)*DEXP(-(S2-0.5)**2/0.03)
00029 5 CONTINUE
00030 DO 554 K=1,N1
00031 DO 554 L=1,N2
00032 S1=(ALIM(9)+ALIM(10))/2.+ST1*(K-N1/2-1)
00033 S2=(ALIM(11)+ALIM(12))/2.+ST2*(L-N2/2-1)
00034 KERN(K,L)=0.
00035 IF(S1.LT.ALIM(9).OR.S1.GT.ALIM(10).
* OR.S2.LT.ALIM(11).OR.S2.GT.ALIM(12))
* GOTO 554
00036 KERN(K,L)=AK(S1,S2)
00037 554 CONTINUE
00038 HH=3.68D-11
00039 DO 771 I1=1,N1
00040 DO 771 I2=1,N2
00041 S=0.
00042 DO 770 J1=1,N1
00043 DO 770 J2=1,N2
00044 IND1=I1-J1+N1/2+1
00045 IND2=I2-J2+N2/2+1
00046 IF(IND1.LT.1.OR.IND1.GT.N1) GOTO 770
00047 IF(IND2.LT.1.OR.IND2.GT.N2) GOTO 770
00048 S=S+KERN(IND1,IND2)*Z0(J1,J2)
00049 770 CONTINUE
00050 U0(I1,I2)=S*ST1*ST2
00051 IF(R.LT.U0(I1,I2)) R=U0(I1,I2)
00052 771 CONTINUE
00053 DL=6.32D-07
00054 PRINT 553,DL,HH
00055 553 FORMAT(' Errors :'/
*' in right side: ',D16.9/
*' in operator: ',D16.9)
00056 CALL PTITR(AK,U0,ALIM,N1,N2,Z,DL,HH,
*C1,ALPHA,AN,U,6500,IMAX,IERR)
00057 PRINT 501,IERR,AN,ALPHA
00058 501 FORMAT(' .',/
*' Return code :',I5/
*' Discrepancy :',D16.9/
*' Regularization parameter :',D16.9)
00059 PRINT 503
00060 PRINT 502,(Z0(13,I),I=1,N2),
*(Z(13,I),I=1,N2)
00061 PRINT 505
00062 PRINT 502,(Z0(I,16),I=1,N1),
*(Z(I,16),I=1,N1)
00063 502 FORMAT(' Exact solution:'/8(4F8.4/))
00064 503 FORMAT('.'/' section X=0.28125 ')
00065 505 FORMAT('.'/' section Y=0.46875 ')
00066 1000 STOP
00067 END


00001 FUNCTION AK(X,Y)
00002 IMPLICIT REAL*8(A-H,O-Z)
00003 AK=DEXP(-20.*(X-0.5)**2-20.*(Y-0.5)**2)
00004 RETURN
00005 END







Errors :
in right side: .632000000D-06
in operator: .368000000D-10
.
Return code : 0
Discrepancy : .782989588D-06
Regularization parameter : .110838081D-07

section X=0.28125
Exact solution:
.0000 .0000 .0000 .0000
.0000 .0000 .0000 .0000
.0006 .0040 .0192 .0706
.2000 .4369 .7354 .9542
.9542 .7354 .4369 .2000
.0706 .0192 .0040 .0006
.0000 .0000 .0000 .0000
.0000 .0000 .0000 .0000

-.0019 -.0033 -.0036 -.0004
.0063 .0120 .0105 -.0009
-.0155 -.0162 .0194 .1083
.2493 .4178 .5706 .6616
.6616 .5706 .4178 .2493
.1083 .0194 -.0162 -.0155
-.0009 .0105 .0120 .0063
-.0004 -.0036 -.0033 -.0019


section Y=0.46875
Exact solution:
.0000 .0000 .0000 .0000
.0000 .0000 .0000 .0000
.0408 .1950 .4586 .7633
.9542 .9152 .7023 .5123
.5123 .7023 .9152 .9542
.7633 .4586 .1950 .0408
.0000 .0000 .0000 .0000
.0000 .0000 .0000 .0000

.0044 .0074 .0075 -.0018
-.0200 -.0353 -.0265 .0277
.1366 .2878 .4498 .5836
.6616 .6810 .6637 .6431
.6431 .6637 .6810 .6616
.5836 .4498 .2878 .1366
.0277 -.0265 -.0353 -.0200
-.0018 .0075 .0074 .0044