         family = <span class="string">'qgauss'</span>;
0295         inc = 2;
0296         <span class="keyword">if</span> (any(size(ydata) ~= [n 1])); error(<span class="string">'y must be n*1 column vector'</span>); <span class="keyword">end</span>;
0297     <span class="keyword">end</span>
0298     <span class="keyword">if</span> (strcmp(varargin{na},<span class="string">'weights'</span>))
0299         wdata = double(varargin{na+1});
0300         inc = 2;
0301         <span class="keyword">if</spa