178
>> dblquad('fof',0,1,1,2)
ans =
1.16777110966887
Натижа: 1.16777110966887
11 - масала.Символлик ҳисоблашлар ёрдамида
,
,
,
,
интегралларни ҳосил қилинг.
Бу ерда z=f(x,y)=xsin(y)+ysin(x) .
Дастур:
syms x y
z=sym('x*sin(y)+y*sin(x)');
i1=int(z,'x')
i2=int(z,'x',0,1)
i3=int(int(z,'x'),'y')
i4=int(int(z,'x',1,2),'y',0,1)
digits(14);
number4=vpa(i4)
i5=int(int(x+y,'y',x,1),'x',0,1)
i1 =
1/2*x^2*sin(y)-y*cos(x)
i2 =
1/2*sin(y)-y*cos(1)+y
i3 =
-1/2*x^2*cos(y)-1/2*y^2*cos(x)
i4 =
-1/2*cos(2)-cos(1)+3/2
number4 =
1.1677711124054
i5 =
1/2
Символлик
ҳисоблашлар
ҳисоблаш
методларининг
хатолигини
бермаганлиги ва ўзлари аниқроқлиги туфайли, d
blquad функция вергулдан кейин
7 рақам аниқликдаги натижани беришини кўриш мумкин.
в)
олий математикадан маълумки, аниқ ва икки каррали интегралларга
кўпгина бошқа
интеграл турлари, масалан, келтирилиши мумкин. Уни топишда
интеграл остида дифференциаллаш
фойдаланилганлиги туфайли, сонли
ҳисоблашларни фойдаланиш нокорректдир.
12 - масала.
1-тур сирт бўйича интегралниҳисобланг: бу ерда
S –
биринчи октантада ётган текистлик қисми (2 теорема бўича).
Дастур:
syms x y z f1 f2
f1=1-x-y;
179
f2=x*y*z;
fun=subs(f2,z,f1)
d=1+diff(f1,x)^2+diff(f1,y)^2
syms x1 x2 y1 y2
x1=sym('0');
x2=sym('1');
y1=sym('0');
y2=sym('1-x');
intpov1=int(int(fun*sqrt(d),'y',y1,y2),'x',x1,x2)
digits(10);
number=vpa(intpov1) fun =
x*y*(1-x-y)
d =
3
intpov1=
1/120*3^(1/2)
number =
1443375673e-1
Натижа: 1443375673e-1
