wolffd@0: function [errorcode,out1,out2,out3,out4] = distchck(nparms,arg1,arg2,arg3,arg4) wolffd@0: %DISTCHCK Checks the argument list for the probability functions. wolffd@0: wolffd@0: % B.A. Jones 1-22-93 wolffd@0: % Copyright (c) 1993-98 by The MathWorks, Inc. wolffd@0: % $Revision: 1.1.1.1 $ $Date: 2005/04/26 02:29:18 $ wolffd@0: wolffd@0: errorcode = 0; wolffd@0: wolffd@0: if nparms == 1 wolffd@0: out1 = arg1; wolffd@0: return; wolffd@0: end wolffd@0: wolffd@0: if nparms == 2 wolffd@0: [r1 c1] = size(arg1); wolffd@0: [r2 c2] = size(arg2); wolffd@0: scalararg1 = (prod(size(arg1)) == 1); wolffd@0: scalararg2 = (prod(size(arg2)) == 1); wolffd@0: if ~scalararg1 & ~scalararg2 wolffd@0: if r1 ~= r2 | c1 ~= c2 wolffd@0: errorcode = 1; wolffd@0: return; wolffd@0: end wolffd@0: end wolffd@0: if scalararg1 wolffd@0: out1 = arg1(ones(r2,1),ones(c2,1)); wolffd@0: else wolffd@0: out1 = arg1; wolffd@0: end wolffd@0: if scalararg2 wolffd@0: out2 = arg2(ones(r1,1),ones(c1,1)); wolffd@0: else wolffd@0: out2 = arg2; wolffd@0: end wolffd@0: end wolffd@0: wolffd@0: if nparms == 3 wolffd@0: [r1 c1] = size(arg1); wolffd@0: [r2 c2] = size(arg2); wolffd@0: [r3 c3] = size(arg3); wolffd@0: scalararg1 = (prod(size(arg1)) == 1); wolffd@0: scalararg2 = (prod(size(arg2)) == 1); wolffd@0: scalararg3 = (prod(size(arg3)) == 1); wolffd@0: wolffd@0: if ~scalararg1 & ~scalararg2 wolffd@0: if r1 ~= r2 | c1 ~= c2 wolffd@0: errorcode = 1; wolffd@0: return; wolffd@0: end wolffd@0: end wolffd@0: wolffd@0: if ~scalararg1 & ~scalararg3 wolffd@0: if r1 ~= r3 | c1 ~= c3 wolffd@0: errorcode = 1; wolffd@0: return; wolffd@0: end wolffd@0: end wolffd@0: wolffd@0: if ~scalararg3 & ~scalararg2 wolffd@0: if r3 ~= r2 | c3 ~= c2 wolffd@0: errorcode = 1; wolffd@0: return; wolffd@0: end wolffd@0: end wolffd@0: wolffd@0: if ~scalararg1 wolffd@0: out1 = arg1; wolffd@0: end wolffd@0: if ~scalararg2 wolffd@0: out2 = arg2; wolffd@0: end wolffd@0: if ~scalararg3 wolffd@0: out3 = arg3; wolffd@0: end wolffd@0: rows = max([r1 r2 r3]); wolffd@0: columns = max([c1 c2 c3]); wolffd@0: wolffd@0: if scalararg1 wolffd@0: out1 = arg1(ones(rows,1),ones(columns,1)); wolffd@0: end wolffd@0: if scalararg2 wolffd@0: out2 = arg2(ones(rows,1),ones(columns,1)); wolffd@0: end wolffd@0: if scalararg3 wolffd@0: out3 = arg3(ones(rows,1),ones(columns,1)); wolffd@0: end wolffd@0: out4 =[]; wolffd@0: wolffd@0: end wolffd@0: wolffd@0: if nparms == 4 wolffd@0: [r1 c1] = size(arg1); wolffd@0: [r2 c2] = size(arg2); wolffd@0: [r3 c3] = size(arg3); wolffd@0: [r4 c4] = size(arg4); wolffd@0: scalararg1 = (prod(size(arg1)) == 1); wolffd@0: scalararg2 = (prod(size(arg2)) == 1); wolffd@0: scalararg3 = (prod(size(arg3)) == 1); wolffd@0: scalararg4 = (prod(size(arg4)) == 1); wolffd@0: wolffd@0: if ~scalararg1 & ~scalararg2 wolffd@0: if r1 ~= r2 | c1 ~= c2 wolffd@0: errorcode = 1; wolffd@0: return; wolffd@0: end wolffd@0: end wolffd@0: wolffd@0: if ~scalararg1 & ~scalararg3 wolffd@0: if r1 ~= r3 | c1 ~= c3 wolffd@0: errorcode = 1; wolffd@0: return; wolffd@0: end wolffd@0: end wolffd@0: wolffd@0: if ~scalararg1 & ~scalararg4 wolffd@0: if r1 ~= r4 | c1 ~= c4 wolffd@0: errorcode = 1; wolffd@0: return; wolffd@0: end wolffd@0: end wolffd@0: wolffd@0: if ~scalararg3 & ~scalararg2 wolffd@0: if r3 ~= r2 | c3 ~= c2 wolffd@0: errorcode = 1; wolffd@0: return; wolffd@0: end wolffd@0: end wolffd@0: wolffd@0: if ~scalararg4 & ~scalararg2 wolffd@0: if r4 ~= r2 | c4 ~= c2 wolffd@0: errorcode = 1; wolffd@0: return; wolffd@0: end wolffd@0: end wolffd@0: wolffd@0: if ~scalararg3 & ~scalararg4 wolffd@0: if r3 ~= r4 | c3 ~= c4 wolffd@0: errorcode = 1; wolffd@0: return; wolffd@0: end wolffd@0: end wolffd@0: wolffd@0: wolffd@0: if ~scalararg1 wolffd@0: out1 = arg1; wolffd@0: end wolffd@0: if ~scalararg2 wolffd@0: out2 = arg2; wolffd@0: end wolffd@0: if ~scalararg3 wolffd@0: out3 = arg3; wolffd@0: end wolffd@0: if ~scalararg4 wolffd@0: out4 = arg4; wolffd@0: end wolffd@0: wolffd@0: rows = max([r1 r2 r3 r4]); wolffd@0: columns = max([c1 c2 c3 c4]); wolffd@0: if scalararg1 wolffd@0: out1 = arg1(ones(rows,1),ones(columns,1)); wolffd@0: end wolffd@0: if scalararg2 wolffd@0: out2 = arg2(ones(rows,1),ones(columns,1)); wolffd@0: end wolffd@0: if scalararg3 wolffd@0: out3 = arg3(ones(rows,1),ones(columns,1)); wolffd@0: end wolffd@0: if scalararg4 wolffd@0: out4 = arg4(ones(rows,1),ones(columns,1)); wolffd@0: end wolffd@0: end wolffd@0: