mtimes.m

matlab 多参数积分工具箱 multivariable calculus toolbox

function C=mtimes(A,B)
%*  Multiply.
%   V*W is the elementwise product of V and W,
%   where V and W are a scalar and a vector function.
%   If both V and W are vector functions, then V*W
%   will be the same as cross(V,W).

% Copyright (c) 2001-04-17, B. Rasmus Anthin.

As=isscalar(A);Bs=isscalar(B);
Av=isvector(A);Bv=isvector(B);
Al=length(A);Bl=length(B);
A=vector(A);
B=vector(B);
if isempty(str2num(A.fx)) &...
      isempty(str2num(A.fy)) &...
      isempty(str2num(A.fz)) &...
      ~strcmp(type(A),type(B))
   B=eval([type(B) '2' type(A) '(B)']);
elseif ~strcmp(type(A),type(B))
   A=eval([type(A) '2' type(B) '(A)']);
end
[x y z]=vars(A);
if (Av | Al==3) & (Bv | Bl==3)
   C=cross(A,B,inputname(1),inputname(2));
else
   C.x=[min(A.x(1),B.x(1)) max(A.x(2),B.x(2)) max(A.x(3),B.x(3))];
   C.y=[min(A.y(1),B.y(1)) max(A.y(2),B.y(2)) max(A.y(3),B.y(3))];
   C.z=[min(A.z(1),B.z(1)) max(A.z(2),B.z(2)) max(A.z(3),B.z(3))];
   C.xval=[];C.yval=[];C.zval=[];
   if ~isempty(A.xval) & ~isempty(B.xval)
      C.xval='*';
   end
   if ~isempty(A.yval) & ~isempty(B.yval)
      C.yval='*';
   end
   if ~isempty(A.zval) & ~isempty(B.zval)
      C.zval='*';
   end
   if isempty(inputname(1))
      Anamex=A.fx;Anamey=A.fy;Anamez=A.fz;
   elseif As
      Anamex=inputname(1);Anamey=inputname(1);Anamez=inputname(1);
   else
      Anamex=[inputname(1) '.' x];Anamey=[inputname(1) '.' y];Anamez=[inputname(1) '.' z];
   end
   if isempty(inputname(2))
      Bnamex=B.fx;Bnamey=B.fy;Bnamez=B.fz;
   elseif Bs
      Bnamex=inputname(2);Bnamey=inputname(2);Bnamez=inputname(2);
   else
      Bnamex=[inputname(2) '.' x];Bnamey=[inputname(2) '.' y];Bnamez=[inputname(2) '.' z];
   end
   findAx=any(Anamex=='+') | any(Anamex=='-');
   findAy=any(Anamey=='+') | any(Anamey=='-');
   findAz=any(Anamez=='+') | any(Anamez=='-');
   findBx=any(Bnamex=='+') | any(Bnamex=='-');
   findBy=any(Bnamey=='+') | any(Bnamey=='-');
   findBz=any(Bnamez=='+') | any(Bnamez=='-');
   if findAx & findBx
      C.fx=['(' Anamex ')*(' Bnamex ')'];
   elseif findAx & ~findBx
      C.fx=['(' Anamex ')*' Bnamex];
   elseif ~findAx & findBx
      C.fx=[Anamex '*(' Bnamex ')'];
   else
      C.fx=[Anamex '*' Bnamex];
   end
   if findAy & findBy
      C.fy=['(' Anamey ')*(' Bnamey ')'];
   elseif findAy & ~findBy
      C.fy=['(' Anamey ')*' Bnamey];
   elseif ~findAy & findBy
      C.fy=[Anamey '*(' Bnamey ')'];
   else
      C.fy=[Anamey '*' Bnamey];
   end
   if findAz & findBz
      C.fz=['(' Anamez ')*(' Bnamez ')'];
   elseif findAz & ~findBz
      C.fz=['(' Anamez ')*' Bnamez];
   elseif ~findAz & findBz
      C.fz=[Anamez '*(' Bnamez ')'];
   else
      C.fz=[Anamez '*' Bnamez];
   end
   C.Fx=['(' A.Fx ').*(' B.Fx ')'];
   C.Fy=['(' A.Fy ').*(' B.Fy ')'];
   C.Fz=['(' A.Fz ').*(' B.Fz ')'];
   C.coords=A.coords;
   C=vector(C);
end