计算两个三角函数乘积的积分,写了如下一个子函数
function y=gam(a,m,n)
% cos(m*pi*x/a)*sin(n*pi*x/a)
if m==0
if n==m
y=0;
else
y=-a*(-1+(-1)^n)/(n*pi);
end
else
if n==m
y=0;
else
y=a*n*(-1+(-1)^(m+n))/(pi*(m^2-n^2));
end
end
想请问一下如何对其进行矢量化,m,n都是矩阵,然后当其中的元素相等的情况下,积分为0,不等的时候为 y=a*n*(-1+(-1)^(m+n))/(pi*(m^2-n^2));。对其进行矢量化后当m=n时,结果为无穷大。