Неравенства для заштрихованной области
Для точек вне окружности радиуса R:
l1=x*x+y*y>=r*r
Для точек внутри квадрата:
l2=x>=-r.and.x<=r.and.y>=-r.and.y<=r
Для точек вне III четверти (т.е. для точек I, II, IV четвертей)
l3=.not.(x<=0.and.y<=0)
Пересечение полученных областей:
l=l1.and.l2.and.l3
Формула для подсчета B(X,Y):
b=exp(-(x*y)/(x-3.))*cos(2.*x-0.3*(10.**(-4.)))+((x*x)/(2.7+x*y))**(1./7.)
одз: x/=3.and.x*y/=-2.7
Текст программы
program lab5
implicit none
real alpha, x, y, b, r/1./
logical l, l1, l2, l3
open(1, file='circle.txt')
open(2, file='square.txt')
open(3, file='inside.txt')
open(4, file='result.txt')
do alpha=0,360,2
x=r*cosd(alpha)
y=r*sind(alpha)
write(1, *) x, y
enddo
do x=-r,r,2.*r/59.
write(2,*) x, r
write(2,*) x, -r
enddo
do y=-r,r,2.*r/59.
write(2,*) -r, y
write(2,*) r, y
enddo
do x=-2.*r,2.*r,4.*r/59
do y=-2.*r,2.*r,4.*r/59
l1=x*x+y*y>=r*r
l2=x>=-r.and.x<=r.and.y>=-r.and.y<=r
l3=.not.(x<=0.and.y<=0)
l=l1.and.l2.and.l3
if(l) then
write(3,*) x,y
if(x/=3.and.x*y/=-2.7) then
b=exp(-(x*y)/(x-3.))*cos(2.*x-0.3*(10.**(-4.)))+((x*x)/(2.7+x*y))**(1./7.)
write(4,*) 'b(',x,'; ',y,')=',b
endif
endif
enddo
enddo
end
Результаты
b( -9.830512E-01; 2.372877E-01)= 5.114928E-01
b( -9.830512E-01; 3.050843E-01)= 5.209492E-01
b( -9.830512E-01; 3.728809E-01)= 5.304161E-01
b( -9.830512E-01; 4.406776E-01)= 5.399021E-01
b( -9.830512E-01; 5.084742E-01)= 5.494164E-01
b( -9.830512E-01; 5.762708E-01)= 5.589689E-01
b( -9.830512E-01; 6.440673E-01)= 5.685703E-01
b( -9.830512E-01; 7.118639E-01)= 5.782319E-01
b( -9.830512E-01; 7.796605E-01)= 5.879667E-01
b( -9.830512E-01; 8.474571E-01)= 5.977883E-01
b( -9.830512E-01; 9.152537E-01)= 6.077119E-01
b( -9.830512E-01; 9.830503E-01)= 6.177543E-01
b( -9.152546E-01; 4.406776E-01)= 6.341263E-01
b( -9.152546E-01; 5.084742E-01)= 6.411638E-01
b( -9.152546E-01; 5.762708E-01)= 6.482534E-01
b( -9.152546E-01; 6.440673E-01)= 6.554029E-01
b( -9.152546E-01; 7.118639E-01)= 6.626207E-01
b( -9.152546E-01; 7.796605E-01)= 6.699157E-01
b( -9.152546E-01; 8.474571E-01)= 6.772979E-01
b( -9.152546E-01; 9.152537E-01)= 6.847780E-01
b( -9.152546E-01; 9.830503E-01)= 6.923679E-01
b( -8.474580E-01; 5.762708E-01)= 7.424949E-01
b( -8.474580E-01; 6.440673E-01)= 7.473198E-01
b( -8.474580E-01; 7.118639E-01)= 7.522199E-01
b( -8.474580E-01; 7.796605E-01)= 7.572017E-01
b( -8.474580E-01; 8.474571E-01)= 7.622719E-01
b( -8.474580E-01; 9.152537E-01)= 7.674379E-01
b( -8.474580E-01; 9.830503E-01)= 7.727079E-01
b( -7.796614E-01; 6.440673E-01)= 8.422801E-01
b( -7.796614E-01; 7.118639E-01)= 8.450403E-01
b( -7.796614E-01; 7.796605E-01)= 8.478854E-01
b( -7.796614E-01; 8.474571E-01)= 8.508196E-01
b( -7.796614E-01; 9.152537E-01)= 8.538482E-01
b( -7.796614E-01; 9.830503E-01)= 8.569762E-01
b( -7.118648E-01; 7.118639E-01)= 9.389609E-01
b( -7.118648E-01; 7.796605E-01)= 9.398923E-01
b( -7.118648E-01; 8.474571E-01)= 9.409123E-01
b( -7.118648E-01; 9.152537E-01)= 9.420240E-01
b( -7.118648E-01; 9.830503E-01)= 9.432307E-01
b( -6.440682E-01; 7.796605E-01)= 1.031030
b( -6.440682E-01; 8.474571E-01)= 1.030398
b( -6.440682E-01; 9.152537E-01)= 1.029854
b( -6.440682E-01; 9.830503E-01)= 1.029399
b( -5.762717E-01; 8.474571E-01)= 1.117013
b( -5.762717E-01; 9.152537E-01)= 1.115108
b( -5.762717E-01; 9.830503E-01)= 1.113285
b( -5.084751E-01; 9.152537E-01)= 1.195439
b( -5.084751E-01; 9.830503E-01)= 1.192570
b( -4.406784E-01; 9.152537E-01)= 1.268365
b( -4.406784E-01; 9.830503E-01)= 1.264790
b( -3.728818E-01; 9.830503E-01)= 1.327280
b( -3.050852E-01; 9.830503E-01)= 1.377039
b( -2.372886E-01; 9.830503E-01)= 1.410396
b( 2.372877E-01; -9.830512E-01)= 1.400240
b( 2.372877E-01; 9.830503E-01)= 1.536382
b( 3.050843E-01; -9.830512E-01)= 1.361850
b( 3.050843E-01; 9.830503E-01)= 1.524928
b( 3.728809E-01; -9.830512E-01)= 1.307299
b( 3.728809E-01; 9.830503E-01)= 1.487375
b( 4.406776E-01; -9.830512E-01)= 1.241032
b( 4.406776E-01; -9.152546E-01)= 1.246023
b( 4.406776E-01; 9.152537E-01)= 1.417778
b( 4.406776E-01; 9.830503E-01)= 1.425602
b( 5.084742E-01; -9.830512E-01)= 1.166846
b( 5.084742E-01; -9.152546E-01)= 1.171208
b( 5.084742E-01; 9.152537E-01)= 1.333189
b( 5.084742E-01; 9.830503E-01)= 1.340941
b( 5.762708E-01; -9.830512E-01)= 1.088178
b( 5.762708E-01; -9.152546E-01)= 1.091418
b( 5.762708E-01; -8.474580E-01)= 1.094783
b( 5.762708E-01; 8.474571E-01)= 1.220752
b( 5.762708E-01; 9.152537E-01)= 1.227568
b( 5.762708E-01; 9.830503E-01)= 1.234532
b( 6.440673E-01; -9.830512E-01)= 1.008201
b( 6.440673E-01; -9.152546E-01)= 1.009818
b( 6.440673E-01; -8.474580E-01)= 1.011566
b( 6.440673E-01; -7.796614E-01)= 1.013443
b( 6.440673E-01; 7.796605E-01)= 1.092005
b( 6.440673E-01; 8.474571E-01)= 1.097020
b( 6.440673E-01; 9.152537E-01)= 1.102177
b( 6.440673E-01; 9.830503E-01)= 1.107478
b( 7.118639E-01; -9.830512E-01)= 9.298452E-01
b( 7.118639E-01; -9.152546E-01)= 9.293516E-01
b( 7.118639E-01; -8.474580E-01)= 9.289814E-01
b( 7.118639E-01; -7.796614E-01)= 9.287319E-01
b( 7.118639E-01; -7.118648E-01)= 9.286008E-01
b( 7.118639E-01; 7.118639E-01)= 9.512073E-01
b( 7.118639E-01; 7.796605E-01)= 9.534684E-01
b( 7.118639E-01; 8.474571E-01)= 9.558400E-01
b( 7.118639E-01; 9.152537E-01)= 9.583232E-01
b( 7.118639E-01; 9.830503E-01)= 9.609190E-01
b( 7.796605E-01; -9.830512E-01)= 8.557795E-01
b( 7.796605E-01; -9.152546E-01)= 8.527161E-01
b( 7.796605E-01; -8.474580E-01)= 8.497552E-01
b( 7.796605E-01; -7.796614E-01)= 8.468915E-01
b( 7.796605E-01; -7.118648E-01)= 8.441200E-01
b( 7.796605E-01; -6.440682E-01)= 8.414365E-01
b( 7.796605E-01; 6.440673E-01)= 8.031211E-01
b( 7.796605E-01; 7.118639E-01)= 8.016261E-01
b( 7.796605E-01; 7.796605E-01)= 8.001733E-01
b( 7.796605E-01; 8.474571E-01)= 7.987619E-01
b( 7.796605E-01; 9.152537E-01)= 7.973909E-01
b( 7.796605E-01; 9.830503E-01)= 7.960594E-01
b( 8.474571E-01; -9.830512E-01)= 7.883811E-01
b( 8.474571E-01; -9.152546E-01)= 7.823378E-01
b( 8.474571E-01; -8.474580E-01)= 7.763596E-01
b( 8.474571E-01; -7.796614E-01)= 7.704370E-01
b( 8.474571E-01; -7.118648E-01)= 7.645612E-01
b( 8.474571E-01; -6.440682E-01)= 7.587240E-01
b( 8.474571E-01; -5.762717E-01)= 7.529176E-01
b( 8.474571E-01; 5.762708E-01)= 6.529185E-01
b( 8.474571E-01; 6.440673E-01)= 6.466583E-01
b( 8.474571E-01; 7.118639E-01)= 6.403257E-01
b( 8.474571E-01; 7.796605E-01)= 6.339161E-01
b( 8.474571E-01; 8.474571E-01)= 6.274250E-01
b( 8.474571E-01; 9.152537E-01)= 6.208476E-01
b( 8.474571E-01; 9.830503E-01)= 6.141795E-01
b( 9.152537E-01; -9.830512E-01)= 7.297008E-01
b( 9.152537E-01; -9.152546E-01)= 7.203331E-01
b( 9.152537E-01; -8.474580E-01)= 7.109750E-01
b( 9.152537E-01; -7.796614E-01)= 7.016109E-01
b( 9.152537E-01; -7.118648E-01)= 6.922261E-01
b( 9.152537E-01; -6.440682E-01)= 6.828068E-01
b( 9.152537E-01; -5.762717E-01)= 6.733398E-01
b( 9.152537E-01; -5.084751E-01)= 6.638125E-01
b( 9.152537E-01; -4.406784E-01)= 6.542126E-01
b( 9.152537E-01; 4.406776E-01)= 5.177940E-01
b( 9.152537E-01; 5.084742E-01)= 5.060382E-01
b( 9.152537E-01; 5.762708E-01)= 4.940500E-01
b( 9.152537E-01; 6.440673E-01)= 4.818185E-01
b( 9.152537E-01; 7.118639E-01)= 4.693331E-01
b( 9.152537E-01; 7.796605E-01)= 4.565825E-01
b( 9.152537E-01; 8.474571E-01)= 4.435558E-01
b( 9.152537E-01; 9.152537E-01)= 4.302415E-01
b( 9.152537E-01; 9.830503E-01)= 4.166281E-01
b( 9.830503E-01; -9.830512E-01)= 6.814286E-01
b( 9.830503E-01; -9.152546E-01)= 6.684728E-01
b( 9.830503E-01; -8.474580E-01)= 6.554526E-01
b( 9.830503E-01; -7.796614E-01)= 6.423438E-01
b( 9.830503E-01; -7.118648E-01)= 6.291232E-01
b( 9.830503E-01; -6.440682E-01)= 6.157688E-01
b( 9.830503E-01; -5.762717E-01)= 6.022597E-01
b( 9.830503E-01; -5.084751E-01)= 5.885754E-01
b( 9.830503E-01; -4.406784E-01)= 5.746962E-01
b( 9.830503E-01; -3.728818E-01)= 5.606028E-01
b( 9.830503E-01; -3.050852E-01)= 5.462762E-01
b( 9.830503E-01; -2.372886E-01)= 5.316979E-01
b( 9.830503E-01; 2.372877E-01)= 4.210525E-01
b( 9.830503E-01; 3.050843E-01)= 4.037956E-01
b( 9.830503E-01; 3.728809E-01)= 3.861194E-01
b( 9.830503E-01; 4.406776E-01)= 3.680043E-01
b( 9.830503E-01; 5.084742E-01)= 3.494304E-01
b( 9.830503E-01; 5.762708E-01)= 3.303775E-01
b( 9.830503E-01; 6.440673E-01)= 3.108249E-01
b( 9.830503E-01; 7.118639E-01)= 2.907515E-01
b( 9.830503E-01; 7.796605E-01)= 2.701356E-01
b( 9.830503E-01; 8.474571E-01)= 2.489551E-01
b( 9.830503E-01; 9.152537E-01)= 2.271875E-01
b( 9.830503E-01; 9.830503E-01)= 2.048093E-01
График
Блок-схема
Санкт-Петербургский Государственный Университет
Информационных технологий,
Механики и оптики
Отчет
К экзаменационной работе по информатике
Студент:
Преподаватель:
2004 г.