P-file 1.0L%6Y4bigint_installansxNLDOS_NAMENEW_FILEBAD_FILEADD_BINAADD_TEXTEND_FILEPROGRESSmfiledosfilemfuncmtext>?\ns:3'?:PCWINfi<findstr @*@$:3'?:VAX_VMS<findstr @*@$:3'?4if findstr(computer, 'PCWIN'), mfile = dosfile; end;3'?[fp, msg] = fopen(mfile, 'w');B3'?if(fp < 0), error(msg), end;3'?fwrite(fp, mtext, 'uchar');3' ? fprintf(fp, mtext);3' ?  fclose(fp);3' ? disp([' Unbundled: ' mfile])3' ?  ##:3-? ' ## Generated by BUNDLE: bigint_install:3-? ## Created: 2000.04.07.08.42.48:3-?7 ## Copyright (C) 1994-9 Dr. Charles R. Denham, ZYDECO.:3-? ## All Rights Reserved.:3-? ## E-Mail: cdenham@usgs.gov. #:3-?- ## WWW: http://crusty.er.usgs.gov/~cdenham/.na:3-? ##:3-?:c: ?" ## BIGINT_INSTALL -- PLEASE NOTE: :3-?) ## This installer creates and populates #:3-?. ## a folder called "bigint", located in your :3-?, ## present directory. After the installer:3-?+ ## has run, be sure to adjust your Matlab:3-?1 ## path to reflect the presence of the "bigint" #:3-? ## folder.:3-? ## Press any key to continue..:3-?pause #;3-? ## mf:3-?  makedir.mdo3'  makedir.mdo3'makedir3'?!%;3-?"%;3-?#%;3-?$function makedir(varargin)%%#%% makedir -- Make a new directory.%B%% makedir('thePath') creates a new sub-directory named 'thePath'%%?%% (no embellishments), located within the current directory.%>%% A dialog is invoked if manual intervention is needed. No%%=%% action is taken if the desired directory already exists.%%%=%% The "present-working-directory" itself remains the same.%%%<%% makedir "p q r" assumes that the arguments represent one%5%% path-name with single-blanks between components.%%%1%% makedir p q r is the same as makedir "p q r".%%%%6%% Copyright (C) 1996-7 Dr. Charles R. Denham, ZYDECO.%%%% All Rights Reserved.%9%% Disclosure without written explicit consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 08-Jul-1997 08:36:09.%#%% Revised 04-Nov-1998 21:57:53.%#%% Revised 05-Nov-1998 15:24:50.%%+if nargin < 1, help(mfilename), return, end%% thePath = '';fo%for i = 1:length(varargin) %) if i > 1, thePath = [thePath ' ']; end %# thePath = [thePath varargin{i}];%end%%6if thePath(1) == '"' & thePath(length(thePath)) == '"' %* thePath = thePath(2:length(thePath)-1);i%end%3'?F% ;3-?G% itExists = 1;tr%trystr = 'cd(thePath); cd ..';c%catchstr = 'itExists = 0;';%eval(trystr, catchstr)i% if itExists%8 disp([' ## Directory already exists: "' thePath '"'])% returnla%end%% lasterr('')%*warn = 'warning(lasterr); lasterr('''');';t%% thePWD = pwd;if%$if thePWD(length(thePWD)) == filesep% thePWD(length(thePWD)) = '';%end%%5if exist('mkdir', 'builtin') | exist('mkdir', 'file') %+ theCommand = ['mkdir(''' thePath ''')'];%& elseif any(findstr(computer, 'MAC')) %/ theCommand = ['newfolder(''' thePath ''')'];%& elseif any(findstr(computer, 'VMS')) %C theCommand = ['!create/directory "' thePWD filesep thePath '"'];% else %8 theCommand = ['!mkdir "' thePWD filesep thePath '"'];%end%%disp([' ## ' theCommand])ev%eval(theCommand, warn)t%%,trystr = 'itExists = 1; cd(thePath); cd ..';%3'?i% ;3-?jcatchstr = 'itExists = 0;';%eval(trystr, catchstr)i%% if itExists%1 disp([' ## Directory created: "' thePath '"'])th% returnla%end%%,thePrompt = ['Make New ' thePath ' Folder'];%%theInstruction = 'Create, Then Save';wh%while ~itExists%. if any(uisetdir(thePrompt, theInstruction)) % cd ..% eval(trystr, catchstr)% else % break% end % thePrompt = [thePrompt '!'];%end%% if itExists%1 disp([' ## Directory created: "' thePath '"'])th% else %0 disp([' ## Unable to create: "' thePath '"'])%end%%-function theResult = newfolder(theFolderName)%%%%.%% newfolder -- Create a new Macintosh folder.%%6%% newfolder('theFolderName') creates a new Macintosh%%8%% folder of 'theFolderName' in the current directory.%5%% The current directory setting remains unchanged.%%%3'?% ;3-?6%% If the folder already exists, no action is taken.%%6%% The result is logical(1) if successful; otherwise%%%% it is logical(0).%% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 05-Sep-1997 17:01:27.%%4%% The applescript commands in file "newfolder.mac".%%%%%0%% set theTarget to thePath & ":" & theNewFolder%&%% if not (exists item theTarget) then%%'%% make new folder at folder thePath%L%% set name of folder "untitled folder" of folder thePath to theNewFolder% %% end ifif%% if nargin < 1 % help(mfilename)i% returnla%end%%!if ~any(findstr(computer, 'MAC')) %N disp([' ## No action taken: "' mfilename '" requires Macintosh computer.'])t% returnla%end%%thePath = pwd;w%)while thePath(length(thePath)) == filesep %! thePath(length(thePath)) = '';th%3'?% ;3-?end%thePath = ['"' thePath '"'];%%'theNewFolder = ['"' theFolderName '"'];%%2result = feval('applescript', 'newfolder.mac', ... %> 'thePath', thePath, 'theNewFolder', theNewFolder);i%%&if ~isempty(result), disp(result), endr%%"result = logical(isempty(result));i%%'if nargout > 0, theResult = result; end%%8function theStatus = uisetdir(thePrompt, theInstruction)%%6%% uisetdir -- Open the destination folder via dialog.%%D%% uisetdir('thePrompt', 'theInstruction') presents the "uiputfile"%D%% dialog with 'thePrompt' and 'theInstruction', for selecting the%C%% desired destination folder. The returned status is logical(1)%*%% if successful; otherwise, logical(0).%% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 03-Jul-1997 09:16:00.%%=if nargin < 1, thePrompt = 'Open The Destination Folder'; endif%3if nargin < 2, theInstruction = 'Save If Okay'; end%%3'?% ;3-?theFile = 0; thePath = 0;[t%:[theFile, thePath] = uiputfile(theInstruction, thePrompt);s%%status = 0; % if isstr(thePath) & any(thePath)% status = 1; %% eval('cd(thePath)', 'status = 0;')if%end%%1if nargout > 0, theStatus = any(any(status)); endne%3'?% ;3-?% ;3-?% ;3-? newfolder.mac453'  newfolder.mac453' newfolder.mac453'?%;3-?%;3-?%;3-?@F@F@*@F@F@@@P@\@Y@@X@@]@Y@@@@X@@@@[@Y@@]@@@Y@[@[@Y@Y@@\@G@*@F@F@*@F@F@@@P@[@\@^@@\@Z@@Y@Z@]@@@D@P@D@@@H@L@L@K@@@Q@\@G@@@P@Z@X@@\@[@Y@@\@@@T@G@@@Q@Y@@[@Z@X@@[@@F@@@V@V@@Q@Q@@P@S@G@*@F@F@@@@@P@@[@[@@@T@Z@@Y@Z@]@\@@@T@Y@@\@Y@@\@]@Y@@Y@G@*@F@F@@@@@@@Q@Z@@\@X@[@[@\@]@@\@Y@@@@]@Z@@]@Z@[@]@@]@@@Y@@^@\@[@Z@@X@Z@@]@@@]@\@Z@@]@]@Y@@[@@@X@[@[@\@Y@@[@]@@@Y@\@[@[@@@@]@Z@Y@@*@F@F@@@@@@@@@X@[@\@^@@\@Z@@Y@Z@]@@@[@]@[@Y@@\@@@Y@[@Y@@\@@@[@[@]@@@X@[@[@\@]@Z@@]@]@@]@Y@@@@\@]@@X@[@Z@@X@X@@]@Z@@[@[@G@*@*@F@F@@@U@Y@@\@\@Z@@[@[@@@[@Y@@@H@J@F@T@Y@@\@F@H@L@L@K@@@H@K@M@J@H@M@H@L@G@*@F@F@*@]@Y@@[@[@@@X@@\@\@[@Z@@X@X@@]@Z@@[@[@@@A@Q@Z@@[@Y@Y@@\@A@*@"@F@F@*@"@\@Y@@]@@@U@Q@@T@U@R@@S@Q@@@]@[@@@H@*@"@\@Y@@]@@@U@Q@@T@U@R@@S@Q@@@]@[@@@H@*@"@Z@@Y@@@U@Q@@T@U@R@@S@Q@@@Z@@\@@@H@@@]@Z@Y@@[@*@"@"@\@Y@@]@@@]@Z@Y@@T@X@@]@Z@@@]@[@@@A@\@\@Z@@X@@\@]@@\@M@P@@P@@P@@A@*@"@"@\@Y@@]@@@]@Z@Y@@S@Y@@]@Q@[@[@Y@Y@@\@@@]@[@@@A@[@@^@@S@Y@@]@Q@[@[@Y@Y@@\@A@*@"@Y@@[@Y@@@Z@@Y@*@"@F@F@*@"@\@Y@@]@@@]@Z@Y@@U@X@@\@Y@Y@@]@@@]@[@@@]@Z@Y@@T@X@@]@Z@@@C@@@A@M@A@@@C@@@]@Z@Y@@S@Y@@]@Q@[@[@Y@Y@@\@*@"@Z@@Y@@@[@[@]@@@D@Y@@^@Z@@\@]@\@@@Z@@]@Y@@[@@@@]@Z@Y@@U@X@@\@Y@Y@@]@D@@@]@Z@Y@@[@*@"@"@[@@X@@Z@Y@@@@[@Y@@]@@@Y@[@[@Y@Y@@\@@@X@@]@@@Y@[@[@Y@Y@@\@@@]@Z@Y@@T@X@@]@Z@*@"@"@\@Y@@]@@@[@X@@[@@Y@@@@[@Y@@@Y@[@[@Y@Y@@\@@@A@]@@[@]@Z@@]@[@Y@@Y@@@Y@[@[@Y@Y@@\@A@@@[@Y@@@Y@[@[@Y@Y@@\@@@]@Z@Y@@T@X@@]@Z@@@]@[@@@]@Z@Y@@S@Y@@]@Q@[@[@Y@Y@@\@*@"@Y@@[@Y@@@Z@@Y@*@"@F@F@*@Y@@[@Y@@@]@Y@@[@[@*3'?% ;3-?% ;3-?% ;3-?makedir('bigint')cd;3-? cd bigintcl;3-?  makedir.mdo3'  makedir.mdo3'makedir3'? %;3-? %;3-? %;3-?function makedir(varargin)%%#%% makedir -- Make a new directory.%B%% makedir('thePath') creates a new sub-directory named 'thePath'%%?%% (no embellishments), located within the current directory.%>%% A dialog is invoked if manual intervention is needed. No%%=%% action is taken if the desired directory already exists.%%%=%% The "present-working-directory" itself remains the same.%%%<%% makedir "p q r" assumes that the arguments represent one%5%% path-name with single-blanks between components.%%%1%% makedir p q r is the same as makedir "p q r".%%%%6%% Copyright (C) 1996-7 Dr. Charles R. Denham, ZYDECO.%%%% All Rights Reserved.%9%% Disclosure without written explicit consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 08-Jul-1997 08:36:09.%#%% Revised 04-Nov-1998 21:57:53.%#%% Revised 05-Nov-1998 15:24:50.%%+if nargin < 1, help(mfilename), return, end%% thePath = '';fo%for i = 1:length(varargin) %) if i > 1, thePath = [thePath ' ']; end %# thePath = [thePath varargin{i}];%end%%6if thePath(1) == '"' & thePath(length(thePath)) == '"' %* thePath = thePath(2:length(thePath)-1);i%end%3'?0% ;3-?1% itExists = 1;tr%trystr = 'cd(thePath); cd ..';c%catchstr = 'itExists = 0;';%eval(trystr, catchstr)i% if itExists%8 disp([' ## Directory already exists: "' thePath '"'])% returnla%end%% lasterr('')%*warn = 'warning(lasterr); lasterr('''');';t%% thePWD = pwd;if%$if thePWD(length(thePWD)) == filesep% thePWD(length(thePWD)) = '';%end%%5if exist('mkdir', 'builtin') | exist('mkdir', 'file') %+ theCommand = ['mkdir(''' thePath ''')'];%& elseif any(findstr(computer, 'MAC')) %/ theCommand = ['newfolder(''' thePath ''')'];%& elseif any(findstr(computer, 'VMS')) %C theCommand = ['!create/directory "' thePWD filesep thePath '"'];% else %8 theCommand = ['!mkdir "' thePWD filesep thePath '"'];%end%%disp([' ## ' theCommand])ev%eval(theCommand, warn)t%%,trystr = 'itExists = 1; cd(thePath); cd ..';%3'?S% ;3-?Tcatchstr = 'itExists = 0;';%eval(trystr, catchstr)i%% if itExists%1 disp([' ## Directory created: "' thePath '"'])th% returnla%end%%,thePrompt = ['Make New ' thePath ' Folder'];%%theInstruction = 'Create, Then Save';wh%while ~itExists%. if any(uisetdir(thePrompt, theInstruction)) % cd ..% eval(trystr, catchstr)% else % break% end % thePrompt = [thePrompt '!'];%end%% if itExists%1 disp([' ## Directory created: "' thePath '"'])th% else %0 disp([' ## Unable to create: "' thePath '"'])%end%%-function theResult = newfolder(theFolderName)%%%%.%% newfolder -- Create a new Macintosh folder.%%6%% newfolder('theFolderName') creates a new Macintosh%%8%% folder of 'theFolderName' in the current directory.%5%% The current directory setting remains unchanged.%%%3'?v% ;3-?w6%% If the folder already exists, no action is taken.%%6%% The result is logical(1) if successful; otherwise%%%% it is logical(0).%% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 05-Sep-1997 17:01:27.%%4%% The applescript commands in file "newfolder.mac".%%%%%0%% set theTarget to thePath & ":" & theNewFolder%&%% if not (exists item theTarget) then%%'%% make new folder at folder thePath%L%% set name of folder "untitled folder" of folder thePath to theNewFolder% %% end ifif%% if nargin < 1 % help(mfilename)i% returnla%end%%!if ~any(findstr(computer, 'MAC')) %N disp([' ## No action taken: "' mfilename '" requires Macintosh computer.'])t% returnla%end%%thePath = pwd;w%)while thePath(length(thePath)) == filesep %! thePath(length(thePath)) = '';th%3'?% ;3-?end%thePath = ['"' thePath '"'];%%'theNewFolder = ['"' theFolderName '"'];%%2result = feval('applescript', 'newfolder.mac', ... %> 'thePath', thePath, 'theNewFolder', theNewFolder);i%%&if ~isempty(result), disp(result), endr%%"result = logical(isempty(result));i%%'if nargout > 0, theResult = result; end%%8function theStatus = uisetdir(thePrompt, theInstruction)%%6%% uisetdir -- Open the destination folder via dialog.%%D%% uisetdir('thePrompt', 'theInstruction') presents the "uiputfile"%D%% dialog with 'thePrompt' and 'theInstruction', for selecting the%C%% desired destination folder. The returned status is logical(1)%*%% if successful; otherwise, logical(0).%% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 03-Jul-1997 09:16:00.%%=if nargin < 1, thePrompt = 'Open The Destination Folder'; endif%3if nargin < 2, theInstruction = 'Save If Okay'; end%%3'?% ;3-?theFile = 0; thePath = 0;[t%:[theFile, thePath] = uiputfile(theInstruction, thePrompt);s%%status = 0; % if isstr(thePath) & any(thePath)% status = 1; %% eval('cd(thePath)', 'status = 0;')if%end%%1if nargout > 0, theStatus = any(any(status)); endne%3'?% ;3-?% ;3-?% ;3-? newfolder.mac453'  newfolder.mac453' newfolder.mac453'?%;3-?%;3-?%;3-?@F@F@*@F@F@@@P@\@Y@@X@@]@Y@@@@X@@@@[@Y@@]@@@Y@[@[@Y@Y@@\@G@*@F@F@*@F@F@@@P@[@\@^@@\@Z@@Y@Z@]@@@D@P@D@@@H@L@L@K@@@Q@\@G@@@P@Z@X@@\@[@Y@@\@@@T@G@@@Q@Y@@[@Z@X@@[@@F@@@V@V@@Q@Q@@P@S@G@*@F@F@@@@@P@@[@[@@@T@Z@@Y@Z@]@\@@@T@Y@@\@Y@@\@]@Y@@Y@G@*@F@F@@@@@@@Q@Z@@\@X@[@[@\@]@@\@Y@@@@]@Z@@]@Z@[@]@@]@@@Y@@^@\@[@Z@@X@Z@@]@@@]@\@Z@@]@]@Y@@[@@@X@[@[@\@Y@@[@]@@@Y@\@[@[@@@@]@Z@Y@@*@F@F@@@@@@@@@X@[@\@^@@\@Z@@Y@Z@]@@@[@]@[@Y@@\@@@Y@[@Y@@\@@@[@[@]@@@X@[@[@\@]@Z@@]@]@@]@Y@@@@\@]@@X@[@Z@@X@X@@]@Z@@[@[@G@*@*@F@F@@@U@Y@@\@\@Z@@[@[@@@[@Y@@@H@J@F@T@Y@@\@F@H@L@L@K@@@H@K@M@J@H@M@H@L@G@*@F@F@*@]@Y@@[@[@@@X@@\@\@[@Z@@X@X@@]@Z@@[@[@@@A@Q@Z@@[@Y@Y@@\@A@*@"@F@F@*@"@\@Y@@]@@@U@Q@@T@U@R@@S@Q@@@]@[@@@H@*@"@\@Y@@]@@@U@Q@@T@U@R@@S@Q@@@]@[@@@H@*@"@Z@@Y@@@U@Q@@T@U@R@@S@Q@@@Z@@\@@@H@@@]@Z@Y@@[@*@"@"@\@Y@@]@@@]@Z@Y@@T@X@@]@Z@@@]@[@@@A@\@\@Z@@X@@\@]@@\@M@P@@P@@P@@A@*@"@"@\@Y@@]@@@]@Z@Y@@S@Y@@]@Q@[@[@Y@Y@@\@@@]@[@@@A@[@@^@@S@Y@@]@Q@[@[@Y@Y@@\@A@*@"@Y@@[@Y@@@Z@@Y@*@"@F@F@*@"@\@Y@@]@@@]@Z@Y@@U@X@@\@Y@Y@@]@@@]@[@@@]@Z@Y@@T@X@@]@Z@@@C@@@A@M@A@@@C@@@]@Z@Y@@S@Y@@]@Q@[@[@Y@Y@@\@*@"@Z@@Y@@@[@[@]@@@D@Y@@^@Z@@\@]@\@@@Z@@]@Y@@[@@@@]@Z@Y@@U@X@@\@Y@Y@@]@D@@@]@Z@Y@@[@*@"@"@[@@X@@Z@Y@@@@[@Y@@]@@@Y@[@[@Y@Y@@\@@@X@@]@@@Y@[@[@Y@Y@@\@@@]@Z@Y@@T@X@@]@Z@*@"@"@\@Y@@]@@@[@X@@[@@Y@@@@[@Y@@@Y@[@[@Y@Y@@\@@@A@]@@[@]@Z@@]@[@Y@@Y@@@Y@[@[@Y@Y@@\@A@@@[@Y@@@Y@[@[@Y@Y@@\@@@]@Z@Y@@T@X@@]@Z@@@]@[@@@]@Z@Y@@S@Y@@]@Q@[@[@Y@Y@@\@*@"@Y@@[@Y@@@Z@@Y@*@"@F@F@*@Y@@[@Y@@@]@Y@@[@[@*3'?% ;3-?% ;3-?% ;3-? makedir.mdo3'  makedir.mdo3'makedir3'?%;3-?%;3-?%;3-?function makedir(varargin)%%#%% makedir -- Make a new directory.%B%% makedir('thePath') creates a new sub-directory named 'thePath'%%?%% (no embellishments), located within the current directory.%>%% A dialog is invoked if manual intervention is needed. No%%=%% action is taken if the desired directory already exists.%%%=%% The "present-working-directory" itself remains the same.%%%<%% makedir "p q r" assumes that the arguments represent one%5%% path-name with single-blanks between components.%%%1%% makedir p q r is the same as makedir "p q r".%%%%6%% Copyright (C) 1996-7 Dr. Charles R. Denham, ZYDECO.%%%% All Rights Reserved.%9%% Disclosure without written explicit consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 08-Jul-1997 08:36:09.%#%% Revised 04-Nov-1998 21:57:53.%#%% Revised 05-Nov-1998 15:24:50.%%+if nargin < 1, help(mfilename), return, end%% thePath = '';fo%for i = 1:length(varargin) %) if i > 1, thePath = [thePath ' ']; end %# thePath = [thePath varargin{i}];%end%%6if thePath(1) == '"' & thePath(length(thePath)) == '"' %* thePath = thePath(2:length(thePath)-1);i%end%3'?% ;3-?% itExists = 1;tr%trystr = 'cd(thePath); cd ..';c%catchstr = 'itExists = 0;';%eval(trystr, catchstr)i% if itExists%8 disp([' ## Directory already exists: "' thePath '"'])% returnla%end%% lasterr('')%*warn = 'warning(lasterr); lasterr('''');';t%% thePWD = pwd;if%$if thePWD(length(thePWD)) == filesep% thePWD(length(thePWD)) = '';%end%%5if exist('mkdir', 'builtin') | exist('mkdir', 'file') %+ theCommand = ['mkdir(''' thePath ''')'];%& elseif any(findstr(computer, 'MAC')) %/ theCommand = ['newfolder(''' thePath ''')'];%& elseif any(findstr(computer, 'VMS')) %C theCommand = ['!create/directory "' thePWD filesep thePath '"'];% else %8 theCommand = ['!mkdir "' thePWD filesep thePath '"'];%end%%disp([' ## ' theCommand])ev%eval(theCommand, warn)t%%,trystr = 'itExists = 1; cd(thePath); cd ..';%3'?;% ;3-?<catchstr = 'itExists = 0;';%eval(trystr, catchstr)i%% if itExists%1 disp([' ## Directory created: "' thePath '"'])th% returnla%end%%,thePrompt = ['Make New ' thePath ' Folder'];%%theInstruction = 'Create, Then Save';wh%while ~itExists%. if any(uisetdir(thePrompt, theInstruction)) % cd ..% eval(trystr, catchstr)% else % break% end % thePrompt = [thePrompt '!'];%end%% if itExists%1 disp([' ## Directory created: "' thePath '"'])th% else %0 disp([' ## Unable to create: "' thePath '"'])%end%%-function theResult = newfolder(theFolderName)%%%%.%% newfolder -- Create a new Macintosh folder.%%6%% newfolder('theFolderName') creates a new Macintosh%%8%% folder of 'theFolderName' in the current directory.%5%% The current directory setting remains unchanged.%%%3'?^% ;3-?_6%% If the folder already exists, no action is taken.%%6%% The result is logical(1) if successful; otherwise%%%% it is logical(0).%% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 05-Sep-1997 17:01:27.%%4%% The applescript commands in file "newfolder.mac".%%%%%0%% set theTarget to thePath & ":" & theNewFolder%&%% if not (exists item theTarget) then%%'%% make new folder at folder thePath%L%% set name of folder "untitled folder" of folder thePath to theNewFolder% %% end ifif%% if nargin < 1 % help(mfilename)i% returnla%end%%!if ~any(findstr(computer, 'MAC')) %N disp([' ## No action taken: "' mfilename '" requires Macintosh computer.'])t% returnla%end%%thePath = pwd;w%)while thePath(length(thePath)) == filesep %! thePath(length(thePath)) = '';th%3'?% ;3-?end%thePath = ['"' thePath '"'];%%'theNewFolder = ['"' theFolderName '"'];%%2result = feval('applescript', 'newfolder.mac', ... %> 'thePath', thePath, 'theNewFolder', theNewFolder);i%%&if ~isempty(result), disp(result), endr%%"result = logical(isempty(result));i%%'if nargout > 0, theResult = result; end%%8function theStatus = uisetdir(thePrompt, theInstruction)%%6%% uisetdir -- Open the destination folder via dialog.%%D%% uisetdir('thePrompt', 'theInstruction') presents the "uiputfile"%D%% dialog with 'thePrompt' and 'theInstruction', for selecting the%C%% desired destination folder. The returned status is logical(1)%*%% if successful; otherwise, logical(0).%% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 03-Jul-1997 09:16:00.%%=if nargin < 1, thePrompt = 'Open The Destination Folder'; endif%3if nargin < 2, theInstruction = 'Save If Okay'; end%%3'?% ;3-?theFile = 0; thePath = 0;[t%:[theFile, thePath] = uiputfile(theInstruction, thePrompt);s%%status = 0; % if isstr(thePath) & any(thePath)% status = 1; %% eval('cd(thePath)', 'status = 0;')if%end%%1if nargout > 0, theStatus = any(any(status)); endne%3'?% ;3-?% ;3-?% ;3-? newfolder.mac453'  newfolder.mac453' newfolder.mac453'?%;3-?%;3-?%;3-?@F@F@*@F@F@@@P@\@Y@@X@@]@Y@@@@X@@@@[@Y@@]@@@Y@[@[@Y@Y@@\@G@*@F@F@*@F@F@@@P@[@\@^@@\@Z@@Y@Z@]@@@D@P@D@@@H@L@L@K@@@Q@\@G@@@P@Z@X@@\@[@Y@@\@@@T@G@@@Q@Y@@[@Z@X@@[@@F@@@V@V@@Q@Q@@P@S@G@*@F@F@@@@@P@@[@[@@@T@Z@@Y@Z@]@\@@@T@Y@@\@Y@@\@]@Y@@Y@G@*@F@F@@@@@@@Q@Z@@\@X@[@[@\@]@@\@Y@@@@]@Z@@]@Z@[@]@@]@@@Y@@^@\@[@Z@@X@Z@@]@@@]@\@Z@@]@]@Y@@[@@@X@[@[@\@Y@@[@]@@@Y@\@[@[@@@@]@Z@Y@@*@F@F@@@@@@@@@X@[@\@^@@\@Z@@Y@Z@]@@@[@]@[@Y@@\@@@Y@[@Y@@\@@@[@[@]@@@X@[@[@\@]@Z@@]@]@@]@Y@@@@\@]@@X@[@Z@@X@X@@]@Z@@[@[@G@*@*@F@F@@@U@Y@@\@\@Z@@[@[@@@[@Y@@@H@J@F@T@Y@@\@F@H@L@L@K@@@H@K@M@J@H@M@H@L@G@*@F@F@*@]@Y@@[@[@@@X@@\@\@[@Z@@X@X@@]@Z@@[@[@@@A@Q@Z@@[@Y@Y@@\@A@*@"@F@F@*@"@\@Y@@]@@@U@Q@@T@U@R@@S@Q@@@]@[@@@H@*@"@\@Y@@]@@@U@Q@@T@U@R@@S@Q@@@]@[@@@H@*@"@Z@@Y@@@U@Q@@T@U@R@@S@Q@@@Z@@\@@@H@@@]@Z@Y@@[@*@"@"@\@Y@@]@@@]@Z@Y@@T@X@@]@Z@@@]@[@@@A@\@\@Z@@X@@\@]@@\@M@P@@P@@P@@A@*@"@"@\@Y@@]@@@]@Z@Y@@S@Y@@]@Q@[@[@Y@Y@@\@@@]@[@@@A@[@@^@@S@Y@@]@Q@[@[@Y@Y@@\@A@*@"@Y@@[@Y@@@Z@@Y@*@"@F@F@*@"@\@Y@@]@@@]@Z@Y@@U@X@@\@Y@Y@@]@@@]@[@@@]@Z@Y@@T@X@@]@Z@@@C@@@A@M@A@@@C@@@]@Z@Y@@S@Y@@]@Q@[@[@Y@Y@@\@*@"@Z@@Y@@@[@[@]@@@D@Y@@^@Z@@\@]@\@@@Z@@]@Y@@[@@@@]@Z@Y@@U@X@@\@Y@Y@@]@D@@@]@Z@Y@@[@*@"@"@[@@X@@Z@Y@@@@[@Y@@]@@@Y@[@[@Y@Y@@\@@@X@@]@@@Y@[@[@Y@Y@@\@@@]@Z@Y@@T@X@@]@Z@*@"@"@\@Y@@]@@@[@X@@[@@Y@@@@[@Y@@@Y@[@[@Y@Y@@\@@@A@]@@[@]@Z@@]@[@Y@@Y@@@Y@[@[@Y@Y@@\@A@@@[@Y@@@Y@[@[@Y@Y@@\@@@]@Z@Y@@T@X@@]@Z@@@]@[@@@]@Z@Y@@S@Y@@]@Q@[@[@Y@Y@@\@*@"@Y@@[@Y@@@Z@@Y@*@"@F@F@*@Y@@[@Y@@@]@Y@@[@[@*3'?% ;3-?% ;3-?% ;3-? clear makedirma;3-?makedir @bigint;3-? cd @biginta;3-?abs.mab3' abs.mab3'abs3'?%;3-?%;3-?%;3-?function theResult = abs(self)%%%2%% bigint/abs -- Square-root of a "bigint" object.%%1%% abs(self) returns the absolute value of self,%%%%% a "bigint" object.% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 16:49:24.%%+if nargin < 1, help(mfilename), return, end%%result = sign(self, +1);%%if nargout > 0 % theResult = result;%else%% assignin('caller', 'ans', result) d% disp(result)an%end%3'?% ;3-?% ;3-?% ;3-?and.man3' and.man3'and3'?%;3-?%;3-?%;3-?%function theResult = and(self, other)%%%%5%% bigint/and -- Logical-AND of two "bigint" objects.%%%0%% and(self, other) returns (self & other), for%8%% the given objects, one of which must be a "bigint".%%%%%%% Also see: or, not.%%% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 13:16:19.%%+if nargin < 2, help(mfilename), return, end%%,result = (abs(self) ~= 0 & abs(other) ~= 0);%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'?% ;3-? % ;3-?!% ;3-?#base.mb3' base.mb3'base3'?$%;3-?%%;3-?&%;3-?'function theResult = base(self)%%1%% bigint/base -- Base of "bigint" (always 1000).%%%3%% base(self) returns the base of self, a "bigint"%&%% object. The base is always 1000.%% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 22-Feb-1998 21:53:06.%%+if nargin < 1, help(mfilename), return, end%%result = 1000; %%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)an%end%3'?A% ;3-?B% ;3-?C% ;3-?Ebigint.m3' bigint.m3'bigintf3'?F%;3-?G%;3-?H%;3-?I function self = bigint(varargin)%%3%% bigint/bigint -- Constructor for "bigint" class.%6%% bigint(theValue) constructs a 'bigint" object from%%9%% theValue, a vector of signed base-1000 coefficients,%%%7%% or a string of decimal digits, possibly signed and%5%% comma-delimited. A "bigint" is a large integer,%%%6%% suitable for representing exact integers that are%%3%% too large for double-precision representation.%-%% They respond to the following operators:%%%-%% + - * / ^ < <= == >= > ~= | & ~%%%8%% Other functions are available; see "methods bigint"%4%% for the complete list. Some functions, such as%5%% division and "sqrt", are computed via bisection.%%%3%% bigint(...) returns a "bigint" constructed from%6%% the left-to-right concatenation of the arguments.%%7%% For a negative-number, use non-negative arguments,%4%% then change the "bigint" sign with unary-minus.%+%% bigint(aBigInt) returns aBigInt intact.%!%% bigint([]) returns bigint(0).%%%'%% bigint('demo') demonstrates itself.%+%% bigint (no argument) returns bigint(0).%%%%%%% See also: methods('bigint').% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 12:37:20.%%3'?k% ;3-?lif nargin < 1 & nargout < 1% eval('setdef(mfilename)', '') % return%end%%"if nargin < 1, varargin = {0}; endi%%if isequal(varargin{1}, 'demo')% help(mfilename)%5 x = bigint(123,456); %% Same as bigint('123,456'). y%1 y = bigint(789); %% Same as bigint('789'). % sum = x + y;% difference = x - y;% product = x * y; q% quotient = x / y; % remainder = rem(x, y); %% greatest_common_divisor = gcd(x, y); l%# least_common_multiple = lcm(x, y);% y_squared = y^2;% y_cubed = y^3; % two_to_the_y_power = 2^y; s% square_root_of_y = sqrt(y);% cube_root_of_y = root(y, 3); % log_base_two_of_y = log2(y);% one_googal = bigint(10)^100; d% disp(display(x)) d% disp(display(y)) d% disp(display(sum))% disp(display(difference)) % disp(display(product))% disp(display(quotient))% disp(display(remainder)) d%3'?% ;3-?' disp(display(greatest_common_divisor))%% disp(display(least_common_multiple)) d% disp(display(y_squared)) d% disp(display(y_cubed))%" disp(display(two_to_the_y_power)) % disp(display(square_root_of_y))% disp(display(cube_root_of_y)) %! disp(display(log_base_two_of_y)) d% disp(display(one_googal))i% return%end%%0if isa(varargin{1}, 'bigint') %% Pass-through.% result = varargin{1};el%Aelseif nargin == 1 & isempty(varargin{1}) %% Same as bigint(0). % result = bigint(0);%$else %% Concatenate the arguments.%! if isa(varargin{1}, 'double') % theSign = +1; % theValue = []; %& for i = 1:prod(size(varargin)) % v = varargin{i};% v = v(:)'; %$ theValue = [theValue v];% end%# elseif isa(varargin{1}, 'char')% theValue = ''; %" for i = 1:length(varargin) %. theValue = [theValue varargin{i}]; % end%- theSign = 1 - 2*any(theValue == '-'); %- theValue = strrep(theValue, '-', ''); %3'?% ;3-?- theValue = strrep(theValue, ',', ''); %+ while any(rem(length(theValue), 3))%& theValue = ['0' theValue]; % end%) v = zeros(1, length(theValue)/3); % k = 0; %$ for i = 1:3:length(theValue)% k = k+1;%) v(k) = eval(theValue(i:i+2)); % end% theValue = v; % end%2%% Define an empty structure for the private data. % theStruct.itsMagnitude = [];% theStruct.itsSign = [];%#%% Allocate the object via "class".%( result = class(theStruct, 'bigint');%1%% Initialize the private data via class methods. %) result = magnitude(result, theValue);%# result = sign(result, theSign);% result = trim(result);%%end%% %% Return. %%if nargout > 0 % self = result;d%else%% assignin('caller', 'ans', result) d% disp(result)an%end%3'?% ;3-?% ;3-?% ;3-?digits.m3' digits.m3'digitsf3'?%;3-?%;3-?%;3-?!function theResult = digits(self)%%%%;%% bigint/digits -- Number of base-10 digits in a "bigint".%<%% digits(self) returns the number of digits in the base-10%/%% representation of self, a "bigint" object.% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 23-Feb-1998 11:32:17.%%+if nargin < 1, help(mfilename), return, end%%a = magnitude(trim(self));r%9result = bigint(length(int2str(a(1))) + 3*(length(a)-1));di%%if nargout > 0 % theResult = result;%else%% assignin('caller', 'ans', result) d% disp(result)%end%3'?% ;3-?% ;3-?% ;3-?disp.mf3' disp.mf3'disp3'?%;3-?%;3-?%;3-?function theResult = disp(self)%%,%% bigint/disp -- Display a "bigint" object.%4%% disp(self) displays or returns the value of self%)%% as text, for self a "bigint" object.%%%%%%%%% Also see: display.%%% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 13:10:01.%%+if nargin < 1, help(mfilename), return, end%%if prod(size(self)) > 1% for i = 1:prod(size(self)) % disp(self(i)) % end% returns%end%%self = trim(self);t%theMagnitude = magnitude(self);%theSign = sign(self);s %s = '';%%if theSign < 0 % s = '-';%else%3'?!% ;3-?" s = '+';%end%%for i = 1:length(theMagnitude) % if i > 1, s = [s, ',']; end%! t = num2str(theMagnitude(i)); % if i > 1%- while length(t) < 3, t = ['0' t]; end % end% s = [s t]; %end%%if nargout > 0 % theResult = s;%else% disp(' ') % disp(s)%end%3'?7% ;3-?8% ;3-?9% ;3-?; display.mdi3'  display.mdi3'display3'?<%;3-?=%;3-?>%;3-??"function theResult = display(self)%%%/%% bigint/display -- Display a "bigint" object.%7%% display(self) displays or returns the value of self%)%% as text, for self a "bigint" object.%%%%%%%%% Also see: display.%%% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 13:10:01.%%+if nargin < 1, help(mfilename), return, end%%if nargout > 0 %- theResult = [inputname(1) ' = ' disp(self)]; d%else% disp(' ') % disp([inputname(1) ' =']) % disp(self)%end%3'?Z% ;3-?[% ;3-?\% ;3-?^double.m3' double.m3'doublef3'?_%;3-?`%;3-?a%;3-?b!function theResult = double(self)%%%%4%% bigint/double -- Double-equivalent of a "bigint".%8%% double(self) returns the double-precision equivalent% %% of self.%%% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 18:41:34.%%+if nargin < 1, help(mfilename), return, end%%5result = sign(self) * polyval(magnitude(self), 1000);eq%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'?|% ;3-?}% ;3-?~% ;3-?eq.m3' eq.m3'eqf3'?%;3-?%;3-?%;3-?$function theResult = eq(self, other)%%9%% bigint/le -- Equal-to comparison for "bigint" objects.%%%:%% le(self, other) returns true if self <= other, for the%%7%% two items, one of which must be a "bigint" object.% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 17:04:31.%%+result = ~(self < other) & ~(self > other);%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'?% ;3-?% ;3-?% ;3-?gcd.mgc3' gcd.mgc3'gcd3'?%;3-?%;3-?%;3-?.function theResult = gcd(self, other, verbose)%%%=%% bigint/gcd -- Greatest-common-divisor of "bigint" objects.%%%8%% gcd(self, other) returns the greatest-common-divisor%=%% between self and other, one of which must be a "bigint".%%%%%%%%% Also see: lcm.%%% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 22-Feb-1998 23:52:25.%%+if nargin < 2, help(mfilename), return, end%if nargin < 3, verbose = 0; end%%a = abs(bigint(self));b%b = abs(bigint(other));%% result = 1;%%if a ~= 0 & b ~= 0 % while (1) % r = rem(a, b);% if r < 2, break, endif% if (verbose)%1 disp([' ## ' mfilename ' remainder = ' disp(r)]) %end% a = b;% b = r;%3'?% ;3-? end%end%% if r == 0 % result = b;%else% result = 1;%end%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)an%end%3'?% ;3-?% ;3-?% ;3-?ge.m3' ge.m3'gef3'?%;3-?%;3-?%;3-?$function theResult = ge(self, other)%%>%% bigint/ge -- Greater-than-or-equal-to for "bigint" objects.%%:%% ge(self, other) returns true if self >= other, for ther%7%% two items, one of which must be a "bigint" object.% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 17:04:31.%%result = ~(self < other);gt%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'?% ;3-?% ;3-?% ;3-?gt.m3' gt.m3'gtf3'?%;3-?%;3-?%;3-?$function theResult = gt(self, other)%%<%% bigint/gt -- Greater-than comparison of "bigint" objects.%9%% gt(self, other) returns true if self > other, for there%7%% two items, one of which must be a "bigint" object.% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 17:04:31.%%result = (other < self);%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'?% ;3-?% ;3-?% ;3-?half.mh3' half.mh3'half3'?%;3-?%;3-? %;3-?!function theResult = half(self)%%%%% bigint/half -- half of a "bigint".%%%"%% half(self) returns (self / 2).%% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 08-Jul-1998 17:23:07.%%+if nargin < 1, help(mfilename), return, end%%%% Use the magnitude.a %%a = magnitude(self);%len = length(a);%%%% Borrow the odd digits.r %%r = rem(a, 2);a% a = a - r;a%%a = a + base(self) .* [0 r(1:len-1)];%%%% %% Halve.re%%!result = magnitude(self, a ./ 2);lc%%if nargout > 0 % theResult = result;%else%3'?C% ;3-?D" assignin('caller', 'ans', result) % disp(result)an%end%3'?J% ;3-?K% ;3-?L% ;3-?Nlcm.mlc3' lcm.mlc3'lcm3'?O%;3-?P%;3-?Q%;3-?R%function theResult = lcm(self, other)%%%%;%% bigint/lcm -- Least-common-multiple of "bigint" objects.%6%% lcm(self, other) returns the least-common-multiple%%=%% between self and other, one of which must be a "bigint",%%%<%% and both of which must be positive numbers. The least-%4%% common-multiple is self*other/gcd(self, other).%%%%%%% Also see: gcd.if% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 22-Feb-1998 23:52:25.%%+if nargin < 2, help(mfilename), return, end%%if self < 1 | other < 1%. error(' ## Both arguments must be positive.')c%end%%c = gcd(self, other);re%%result = (self * other) / c;%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)an%3'?t% ;3-?uend%3'?y% ;3-?z% ;3-?{% ;3-?}le.m3' le.m3'lef3'?~%;3-?%;3-?%;3-?$function theResult = le(self, other)%%;%% bigint/le -- Less-than-or-equal-to for "bigint" objects.%:%% le(self, other) returns true if self <= other, for the%%7%% two items, one of which must be a "bigint" object.% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 17:04:31.%%result = ~(self > other);lo%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'?% ;3-?% ;3-?% ;3-?log10.m3' log10.m3'log10fu3'?%;3-?%;3-?%;3-? function theResult = log10(self)%%:%% bigint/log10 -- Base-10 logarithm of a "bigint" object.%%=%% log10(self) returns the base-10 logarithm of the absolute%%%/%% value of self, a non-zero "bigint" object.%%%%%%% Also see: log2, digits.% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 16:49:24.%%+if nargin < 1, help(mfilename), return, end%%result = digits(abs(self)) - 1;%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'?% ;3-?% ;3-?% ;3-?log2.ml3' log2.ml3'log23'?%;3-?%;3-?%;3-?function theResult = log2(self)%%8%% bigint/log2 -- Base-2 logarithm of a "bigint" object.%;%% log2(self) returns the base-2 logarithm of the absolute%:%% value of self, a "bigint" object. The algorithm uses%%%% bisection.%%%%%"%% Also see: root, mpower, log10.s% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 16:49:24.%%+if nargin < 1, help(mfilename), return, end%%self = abs(bigint(self));re%%result = bigint(1);%%two = bigint(2);%%while two^(2*result) <= self% result = 2*result;d%end%%delta = result;%%while delta > 1% delta = half(delta); t%3'?% ;3-? trial = result + delta;% if two^trial <= self % result = trial; e% end%end%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'?% ;3-?% ;3-?% ;3-?lt.m3' lt.m3'ltf3'?%;3-?%;3-?%;3-?$function theResult = lt(self, other)%%9%% bigint/lt -- Less-than comparison of "bigint" objects.%%%9%% lt(self, other) returns true if self < other, for the%%%7%% two items, one of which must be a "bigint" object.%4%% This routine is the basis for all other logical%%% comparisons.re% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 17:04:31.%% result = 0;%%self = bigint(self);%other = bigint(other);i%%if sign(self) < sign(other)% result = 1;%else% a = magnitude(self);% b = magnitude(other); %/ while length(a) < length(b), a = [0 a]; end%/ while length(b) < length(a), b = [0 b]; end% f = find(a < b);% g = find(a > b);% if any(f) & ~any(g)% result = 1;%( elseif any(f) & any(g) & f(1) < g(1)%3'?!% ;3-?" result = 1;% else% result = 0;% end%end%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'?1% ;3-?2% ;3-?3% ;3-?5 magnitude.m3'  magnitude.m3' magnitudefu3'?6%;3-?7%;3-?8%;3-?92function theResult = magnitude(self, theMagnitude)%%%6%% bigint/magnitude -- Magnitude of a "bigint" object.%%2%% magnitude(self) returns the magnitude of self,%%4%% expressed in base-1000 as a vector of "double".%4%% magnitude(self, theMagnitude) sets the magnitude%%% of self to theMagnitude.%%%%%%%%% Also see: sign.%% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 13:23:47.%% if nargin < 2 % result = self.itsMagnitude;%else%% self.itsMagnitude = theMagnitude; % result = self;m%end%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)an%end%3'?Z% ;3-?[% ;3-?\% ;3-?^minus.m3' minus.m3'minusfu3'?_%;3-?`%;3-?a%;3-?b'function theResult = minus(self, other)%%6%% bigint/minus -- Difference of two "bigint" objects.%%2%% minus(self, other) returns (self - other), for%%8%% the given objects, one of which must be a "bigint".%%%%%%% Also see: plus.[% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 13:16:19.%%+if nargin < 2, help(mfilename), return, end%%self = bigint(self);%other = bigint(other);i%%![self, other] = pad(self, other);a %%!a = magnitude(self) * sign(self);b %#b = magnitude(other) * sign(other);%%if self < other% result = -bigint(b - a);%else% result = bigint(a - b);%end%%if nargout > 0 %3'?% ;3-? theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'?% ;3-?% ;3-?% ;3-?mpower.m3' mpower.m3'mpowerf3'?%;3-?%;3-?%;3-?+function theResult = mpower(self, thePower)%%9%% bigint/mpower -- Raise a "bigint" to an integer power.%%%;%% mpower(self, thePower) returns self raised to thePower,%@%% for self, a "bigint", and thePower, a non-negative integer.%9%% This routine uses the successive-squaring algorithm.%%%%%%%%% Also see: sqrt, root.% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 22-Feb-1998 15:24:29.%%+if nargin < 2, help(mfilename), return, end%%self = bigint(self);%%thePower = abs(thePower);if%%if thePower < 1%$ [num, den] = rat(double(thePower));%$ result = root(self^num(1), den(1));%else% self = bigint(self); t% thePower = bigint(thePower); r% result = bigint(1);% factor = self;% while thePower > 0% remainder = rem(thePower, 2);%3'?% ;3-? if remainder % result = result * factor;% end %( thePower = half(thePower - remainder);% factor = factor * factor;% end%end%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'?% ;3-?% ;3-?% ;3-? mrdivide.mm3'  mrdivide.mm3'mrdivide3'?%;3-?%;3-?%;3-?:function [theResult, theRemainder] = mrdivide(self, other)%%%7%% bigint/mrdivide -- Quotient of two "bigint" objects.%5%% mrdivide(self, other) returns (self / other), for%%%8%% the given objects, one of which must be a "bigint".%"%% The algorithm uses bisection.%%5%% [theResult, theRemainder] = mrdivide(self, other)%%%6%% also returns theRemainder when called explicitly.%%,%% Otherwise, use "rem" for the remainder.%%%%%%% Also see: mtimes.% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 13:16:19.%%+if nargin < 2, help(mfilename), return, end%%self = bigint(self);%other = bigint(other);i%%!theSign = sign(self)*sign(other);se%%self = abs(self);ot%other = abs(other);%%result = bigint(0);%%if self >= other%3'?% ;3-? result = bigint(1);% while result*other*2 <= self % result = result*2;% end% delta = result;% while delta > 1% delta = half(delta); % trial = result + delta; % if trial*other <= self% result = trial; % end % end% result = sign(result, theSign);%end%%%if nargout > 0 % theResult = result;%8 if nargout > 1, theRemainder = self - other*result; end%else%" assignin('caller', 'ans', result) % disp(result)%end%3'?% ;3-?% ;3-?% ;3-?mtimes.m3' mtimes.m3'mtimesf3'?%;3-?%;3-?%;3-?(function theResult = mtimes(self, other)%%4%% bigint/mtimes -- Product of two "bigint" objects.%3%% mtimes(self, other) returns (self * other), for%8%% the given objects, one of which must be a "bigint".%%%%%%% Also see: mrdivide.c% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 13:16:19.%%+if nargin < 2, help(mfilename), return, end%%self = bigint(self);%other = bigint(other);i%%,c = conv(magnitude(self), magnitude(other));%%1result = sign(bigint(c), sign(self)*sign(other));ne%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'?7% ;3-?8% ;3-?9% ;3-?;ne.m3' ne.m3'nef3'?<%;3-?=%;3-?>%;3-??$function theResult = ne(self, other)%%=%% bigint/ne -- Not-equal-to comparison for "bigint" objects.%%%:%% ne(self, other) returns true if self ~= other, for ther%7%% two items, one of which must be a "bigint" object.% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 17:04:31.%%result = ~(self == other);n%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'?W% ;3-?X% ;3-?Y% ;3-?[not.mno3' not.mno3'not3'?\%;3-?]%;3-?^%;3-?_function theResult = not(self)%%%2%% bigint/not -- Logical-NOT of a "bigint" object.%%8%% not(self, other) returns ~self, for self a "bigint".%%%%%%% Also see: and, or.re% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 13:16:19.%%+if nargin < 1, help(mfilename), return, end%%result = (abs(self) == 0);o%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'?z% ;3-?{% ;3-?|% ;3-?~or.m3' or.m3'orf3'?%;3-?%;3-?%;3-?!function result = or(self, other)%%%%3%% bigint/or -- Logical-OR of two "bigint" objects.%/%% or(self, other) returns (self & other), for%8%% the given objects, one of which must be a "bigint".%%%%%%% Also see: and, not.r% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 13:16:19.%%+if nargin < 2, help(mfilename), return, end%%,result = (abs(self) ~= 0 | abs(other) ~= 0);%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'?% ;3-?% ;3-?% ;3-?pad.mpa3' pad.mpa3'pad3'?%;3-?%;3-?%;3-?1function [theResult, theOther] = pad(self, other)%%%%9%% bigint/pad -- Pad two "bigint" objects to same length.%%%5%% [theResult, theOther] = pad(self, other) pads the%%%5%% magnitudes of self and other to the same length.%%%6%% At least one of the arguments must be a "bigint".%%%%%%%% Also see: trim.%% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 22-Feb-1998 14:41:49.%%self = bigint(self);%other = bigint(other);i%%a = magnitude(self);%b = magnitude(other);wh%%+while length(a) < length(b), a = [0 a]; end%+while length(b) < length(a), b = [0 b]; end%%self = magnitude(self, a);o%other = magnitude(other, b);%%if nargout > 0 % theResult = self; % theOther = other; r%else%3'?% ;3-? result = {result, other};p%% assignin('caller', 'ans', result) d% disp(result)an%end%3'?% ;3-?% ;3-?% ;3-?plus.mp3' plus.mp3'plus3'?%;3-?%;3-?%;3-?&function theResult = plus(self, other)%%%.%% bigint/plus -- Sum of two "bigint" objects.%%1%% plus(self, other) returns (self + other), for%%%8%% the given objects, one of which must be a "bigint".%%%%%%% Also see: minus.% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 13:16:19.%%+if nargin < 2, help(mfilename), return, end%%self = bigint(self);%other = bigint(other);i%%![self, other] = pad(self, other);a %%!a = magnitude(self) * sign(self);b %#b = magnitude(other) * sign(other);%%result = bigint(a + b);%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%3'?% ;3-?end%3'?% ;3-?% ;3-?% ;3-?rem.mre3' rem.mre3'rem3'?%;3-?%;3-?%;3-?*function theResult = rem(self, theModulus)%%%/%% bigint/rem -- Remainder for "bigint" object.%4%% rem(self, theModulus) returns the remainder from%3%% integer division of self by theModulus, one of% %% which is a "bigint" object.% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 22-Feb-1998 17:26:42.%%+if nargin < 2, help(mfilename), return, end%%quotient = self / theModulus;te%%$temp = self - quotient * theModulus;% if temp < 0% temp = self + temp;%! quitient = temp / theModulus;re%end%%&result = self - quotient * theModulus;r%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'?*% ;3-?+% ;3-?,% ;3-?.root.mr3' root.mr3'root3'?/%;3-?0%;3-?1%;3-?2(function theResult = root(self, theRoot)%%4%% bigint/root -- Integer root of a "bigint" object.%9%% root(self, theRoot) returns the desired root of self,%%%8%% a "bigint" object, for theRoot, a positive-integer.%"%% The algorithm uses bisection.%%%%%%%% Also see: sqrt, mpower.t% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 22-Feb-1998 15:24:29.%%+if nargin < 2, help(mfilename), return, end%%self = abs(bigint(self));re%theRoot = abs(bigint(theRoot));%%result = bigint(1);% while (result*2)^theRoot <= self% result = result*2;%end%%delta = result;%%while delta > 1% delta = half(delta); t% trial = result + delta;% if trial^theRoot <= selfsi%3'?T% ;3-?U result = trial; e% end%end%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'?b% ;3-?c% ;3-?d% ;3-?fsign.ms3' sign.ms3'sign3'?g%;3-?h%;3-?i%;3-?j(function theResult = sign(self, theSign)%%,%% bigint/sign -- Sign of a "bigint" object.%5%% sign(self) returns -1 or +1 for the sign of self,%%%%% a "bigint" object.%2%% magnitude(self, theSign) sets the sign of self%%%% to theSign.%%%%%%% Also see: magnitude.% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 13:23:47.%% if nargin < 2 % if self.itsSign < 0% result = -1;% else% result = +1;% end%else% if theSign < 0 % theSign = -1; % else% theSign = +1; % end% self.itsSign = theSign;% result = self;m%end%3'?% ;3-?%if nargout > 0 % theResult = result;%else%% assignin('caller', 'ans', result) d% disp(result)an%end%3'?% ;3-?% ;3-?% ;3-?sqrt.ms3' sqrt.ms3'sqrt3'?%;3-?%;3-?%;3-?function theResult = sqrt(self)%%3%% bigint/sqrt -- Square-root of a "bigint" object.%0%% sqrt(self) returns the largest integer whose%7%% square does not exceed the absolute value of self,%6%% a "bigint" object. The algorithm uses bisection.%%%%%%%% Also see: root, mpower.% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 16:49:24.%%+if nargin < 1, help(mfilename), return, end%%self = abs(bigint(self));re%%result = bigint(1);%while result*result*4 <= self i% result = result*2;%end%%delta = result;%%while delta > 1% delta = half(delta); t% trial = result + delta;% if trial^2 <= self% result = trial; e%3'?% ;3-? end%end%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'?% ;3-?% ;3-?% ;3-?trim.mt3' trim.mt3'trim3'?%;3-?%;3-?%;3-?function theResult = trim(self)%%5%% bigint/trim -- Trim leading zeros from a "bigint".%%%2%% trim(self) adjusts self, a "bigint" object, to%%2%% mod(1000) form and removes its leading zeros.%%%%%%%% Also see: pad.%%% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 21-Feb-1998 13:39:08.%%+if nargin < 1, help(mfilename), return, end%%%% Trim.%%a = magnitude(self);%/while length(a) > 1 & a(1) == 0, a(1) = []; end%%%% Change sign.%%-if a(1) < 0, self = -magnitude(self, -a); end%%%%a = magnitude(self);%% %% Borrow.t%%theBase = base(self);wh%while any(a < 0)%3'?% ;3-?0 while length(a) > 1 & a(1) == 0, a(1) = []; end% for i = 2:length(a)% a(i-1) = a(i-1) - 1; % a(i) = a(i) + theBase;% end%end%% %% Carry.ca%% carry = 0;f%for i = length(a):-1:1 % a(i) = a(i) + carry;%# remainder = rem(a(i), theBase);%) carry = (a(i) - remainder) / theBase; % a(i) = remainder;%%%end%%%% Left-over carry.%%while carry > 0% a = [carry a]; %# remainder = rem(a(1), theBase);%) carry = (a(1) - remainder) / theBase; % a(1) = remainder;re%end%%%% Trim.%%/while length(a) > 1 & a(1) == 0, a(1) = []; end%%result = magnitude(self, a);%%3'? % ;3-? if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'? %% ;3-? &% ;3-? '% ;3-? )uminus.m3' uminus.m3'uminusf3'? *%;3-? +%;3-? ,%;3-? -!function theResult = uminus(self)%%%%5%% bigint/uminus -- Unary-minus of a "bigint" object.%%%5%% uminus(self) changes the sign of self, a "bigint"%%% %% object.%%%%%%% Also see: uplus.% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 22-Feb-1998 14:33:07.%%!result = sign(self, -sign(self));up%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'? G% ;3-? H% ;3-? I% ;3-? Kuplus.m3' uplus.m3'uplusfu3'? L%;3-? M%;3-? N%;3-? O function theResult = uplus(self)%%3%% bigint/uplus -- Unary-plus of a "bigint" object.%0%% uplus(self) returns self, a "bigint" object.%%%%%%% Also see: uminus.% %%%4%% Copyright (C) 1997 Dr. Charles R. Denham, ZYDECO.%%% All Rights Reserved.%9%% Disclosure without explicit written consent from the%%%6%% copyright owner does not constitute publication. % %%%#%% Version of 22-Feb-1998 14:33:07.%%result = self;c%%if nargout > 0 % theResult = result;%else%" assignin('caller', 'ans', result) % disp(result)%end%3'? h% ;3-? i% ;3-? j% ;3-? kcd .. #;3-? lcd .. #;3-? m %%:3-? n' ## REMEMBER to put the "bigint" folder:3-? o ## in your Matlab path. #:3-? p %%:3-? q/ ## NOTE: "bigint('demo')" demonstrates itself.:3-G