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

Fig 4.3
=======



00001 IMPLICIT REAL*8(A-H,O-Z)
00002 IMPLICIT INTEGER*4(I-N)
00003 DIMENSION U0(41),Z(41),U(5000),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 ALFA=0.000001
00014 HX=(X2-X1)/(N-1.)
00015 CALL PTICR0(AK,U,X1,X2,Y1,Y2,N,M)
00016 H=1.E-10
00017 DO 57 I=1,N
00018 X=X1+HX*(I-1.)
00019 Z0(I)=(DEXP(-(X-0.5)**2/0.06))
00020 57 CONTINUE
00021 PRINT 501,(Z0(II),II=1,N)
00022 501 FORMAT(1X,'Exact solution='/(5F11.7))
00023 CALL PTICR3(U,Z0,U0,N,M)
00024 DL=1.E-8
00025 PRINT 502,H,DL
00026 502 FORMAT(' .'/
*' Errors: in operator -',D16.9/
*' in right side -',D16.9)
00027 DO 34 I=1,N
00028 34 Z(I)=0.
00029 CALL PTIZR(AK,U0,X1,X2,Y1,Y2,N,M,Z,
*IC,AN2,DL,H,C1,ANGRD,IMAX,ALFA,U,
*10000,IERR)
00030 PRINT 503,(Z(II),II=1,N)
00031 PRINT 504,IERR,AN2,ALFA
00032 503 FORMAT(' .'/
*' Approximate solution:'/(5F11.7))
00033 504 FORMAT(' .'/
*' Return code :',I5/
*' Discrepancy :',D16.9/
*' Regularization parameter :',D16.9)
00034 STOP
00035 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 - .999999994D-08
.
Approximate solution:
.0211643 .0240784 .0328794 .0475972 .0684072
.0957261 .1302309 .1728098 .2243091 .2850577
.3545593 .4316402 .5146232 .6011371 .6878841
.7708297 .8457517 .9088284 .9568809 .9871938
.9975960 .9871938 .9568809 .9088284 .8457517
.7708297 .6878841 .6011371 .5146232 .4316402
.3545593 .2850577 .2243092 .1728098 .1302309
.0957261 .0684072 .0475972 .0328794 .0240785
.0211644
.
Return code : 0
Discrepancy : .151755068D-07
Regularization parameter : .929922913D-06