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

Fig 4.2
=======


00001 IMPLICIT REAL*8(A-H,O-Z)
00002 IMPLICIT INTEGER*4(I-N)
00003 DIMENSION Z0(41),U0(41),U(5000),Z(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.0004
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.3)**2/0.03)+
+DEXP(-(X-0.7)**2/0.03))/0.955040800-
-0.0521309113
00020 57 CONTINUE
00021 CALL PTICR3(U,Z0,U0,N,M)
00022 PRINT 501,(Z0(II),II=1,N)
00023 501 FORMAT(1X,'Exact solution='/(5F11.7))
00024 DL=1.E-8
00025 PRINT 502,H,DL
00026 502 FORMAT(' .'/
*' Errors: in operator -',D16.9/
*' in right side -',D16.9)
00027 CALL PTIMR(AK,U0,X1,X2,Y1,Y2,N,M,Z,AN2,
*DL,H,C1,IMAX,ALFA,U,10000,IERR)
00028 PRINT 503,(Z(II),II=1,N)
00029 PRINT 504,IERR,AN2,ALFA
00030 503 FORMAT(' .'/
* ' Approximate solution:'/(5F11.7))
00031 504 FORMAT(' .'/
* ' Return code :',I5/
* ' Discrepancy :',D16.9/
* ' Regularization parameter :',D16.9)
00032 STOP
00033 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=
.0000000 .0320465 .0782459 .1415609 .2238816
.3251418 .4425163 .5699656 .6983832 .8164924
.9124512 .9758983 1.0000000 .9829996 .9288697
.8468924 .7502625 .6540357 .5728488 .5188144
.4998815 .5188144 .5728488 .6540357 .7502625
.8468924 .9288697 .9829996 1.0000000 .9758983
.9124512 .8164924 .6983832 .5699656 .4425163
.3251418 .2238816 .1415609 .0782459 .0320465
.0000000
.
Errors: in operator - .100000001D-09
in right side - .999999994D-08
.
Approximate solution:
.0173081 .0311578 .0715787 .1366338 .2238399
.3296361 .4487824 .5745464 .6991498 .8137768
.9084307 .9725681 .9974897 .9802096 .9256576
.8447169 .7506126 .6567675 .5762364 .5213060
.5017045 .5213060 .5762364 .6567675 .7506126
.8447169 .9256576 .9802096 .9974897 .9725681
.9084308 .8137768 .6991498 .5745464 .4487824
.3296361 .2238399 .1366338 .0715788 .0311579
.0173083
.
Return code : 0
Discrepancy : .178224579D-07
Regularization parameter : .244141280D-06