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

Fig 4.4
=======



00001 IMPLICIT REAL*8(A-H,O-Z)
00002 IMPLICIT INTEGER*4(I-N)
00003 DIMENSION U0(41),Z(41),U(10000),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=1000
00012 C1=1.001
00013 IC=0
00014 ALFA=0.000001
00015 HX=(X2-X1)/(N-1.)
00016 CALL PTICR0(AK,U,X1,X2,Y1,Y2,N,M)
00017 H=1.E-10
00018 DO 57 I=1,N
00019 X=X1+HX*(I-1.)
00020 Z0(I)=(DEXP(-(X-0.5)**2/0.06))
00021 57 CONTINUE
00022 PRINT 501,(Z0(II),II=1,N)
00023 501 FORMAT(1X,'Exact solution='/(5F11.7))
00024 CALL PTICR3(U,Z0,U0,N,M)
00025 DL=1.D-8
00026 PRINT 502,H,DL
00027 502 FORMAT(' .'/
*' Errors: in operator -',D16.9/
*' in right side -',D16.9)
00028 DO 34 I=1,N
00029 34 Z(I)=0.5
00030 CALL PTIPR(AK,U0,X1,X2,Y1,Y2,N,M,Z,IC,
*AN2,DL,H,C1,ANGRD,IMAX,ALFA,U,
*10000,IERR)
00031 PRINT 503,(Z(II),II=1,N)
00032 PRINT 504,IERR,AN2,ALFA
00033 503 FORMAT(' .'/
*' Approximate solution:'/(5F11.7))
00034 504 FORMAT(' .'/
*' Return code :',I5/
*' Discrepance :',D16.9/
*' Regularization parameter :',D16.9)
00035 STOP
00036 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=
.0155039 .0232740 .0342181 .0492711 .0694834
.0959671 .1298122 .1719732 .2231302 .2835359
.3528661 .4300946 .5134171 .6002454 .6872893
.7707304 .8464817 .9105104 .9591895 .9896374
1.0000000 .9896374 .9591895 .9105104 .8464817
.7707304 .6872893 .6002454 .5134171 .4300946
.3528661 .2835359 .2231302 .1719732 .1298122
.0959671 .0694834 .0492711 .0342181 .0232740
.0155039
.
Errors: in operator - .100000001D-09
in right side - .100000000D-07
.
Approximate solution:
.0210491 .0239928 .0328555 .0476350 .0684788
.0957898 .1302494 .1727642 .2242018 .2849067
.3543894 .4314756 .5144849 .6010413 .6878428
.7708505 .8458373 .9089764 .9570818 .9874308
.9978459 .9874308 .9570818 .9089764 .8458374
.7708505 .6878428 .6010413 .5144849 .4314756
.3543894 .2849068 .2242018 .1727642 .1302494
.0957897 .0684788 .0476350 .0328555 .0239929
.0210492
.
Return code : 67
Discrepancy : .124056379D-07
Regularization parameter : .838278424D-06