[f,args]=Ms1v(z,m,q,ntrms)

function [f,args]=Ms1v(z,m,q,ntrms)
% [f,args]=Ms1v(z,m,q,ntrms)

% 3D array form of the modified Mathieu
% function Ms1(z,m,q) defined by 
% 20.6.9 & 20.6.10 of
% 'Handbook of Mathematical Functions'
% by Milton Abramowitz and Irene Stegun
% z     - vector of argument values
% m     - integer function order
% q     - scalar parameter value
% ntrms - number of series terms used
% f       array in which f(i,j,k)
%         contains Ms1(z,m,q) eval-
%         uated at z(i),m(j),q(k)
% args    a structure returning the
%         input data in elements
%         args.z, args.m, args.q

%          HBW 10/10/04

if nargin<4, ntrms=50; end
z=z(:); nz=length(z);
m=m(:)'; M=length(m); nq=length(q);
f=zeros(nz,M,nq);
for j=1:M, for k=1:nq
  f(:,j,k)=Ms1(z,m(j),q(k),ntrms);  
end, end
if nargout==2
  args=struct('z',{z},'m',{m},'q',{q});
end