%задание точок
a=4;
b=4;
u = (-0.5:0.1:0.5)';
v = 0:0.2*pi:2*pi;
X = a*ones(size(u))*cos(v);
Y = b*ones(size(u))*sin(v)+15;
Z = u*ones(size(v))+5;
% figure('Color','w')
% hS=mesh(X,Y,Z);
xlabel('ось X'),ylabel('ось Y'),zlabel('ось Z')
%обчислення інтегралів
for i=1:11
for j=1:11
par1=X(i,j);
par2=Y(i,j);
par3=Z(i,j);
XX(i,j) = quad (@myX,0,1, 1.0e-5, 0, par1, par2,par3);
YY(i,j) = quad (@myY,0,1, 1.0e-5, 0, par1, par2,par3);
ZZ(i,j) = quad (@myZ,0,1, 1.0e-5, 0, par1, par2,par3);
end
end
grid on
%випромінювач
x1=0;y1=15;z1=5;
m=0.1; n=0; p=0;
t=0:0.05:1; x=x1+m*t;y=y1+n*t;z=z1+p*t;
plot3(x,y,z)
xlabel('ось X'),ylabel('ось Y'),zlabel('ось Z')
hold on
quiver3 (X,Y,Z,XX,YY,ZZ)