clear all;
syms a;
syms x ;
syms y ;
a=0:0.1:pi/2;
x=0:0.1:pi/2;
y=0:0.1:pi/2;
P=cos(x).*cos(y)+sin(x).*sin(y).*cos(a);
Q=(cos(x)-cos(y)).^2+(sin(x).*cos(a)-sin(y)).^2+(sin(x).*sin(a)).^2;
H=P./sqrt(Q);
M=int(int(H,y),x);
plot(M,a)
这是我的源代码,函数,变量什么的均有定义,不知道哪点出了错。