Контроль точности при повторении хаотических карт в Matlab

Question

Я хотел бы иметь возможность точно контролировать десятичную точность при повторении хаотической карты. В частности, я хотел бы убедиться, что Matlab использует не больше / не меньше значения p, которое я указываю в следующей функции:

function out = chaotic_logistic(seed,n,p,varargin)
%Returns a vector containing the seed and n iterations 
%of the logistic map: u*x*(1-x), which is chaotic on the
%unit interval when u is contained in (~3.56995,4)

%If a third argument is passed, this will be used as the u parameter.
%Otherwise u=4 by default.

%Finally, if 'last' is passed as the fourth argument, only the last
%iteration will be returned.

%NOTE: As currently designed, this function can be broken 
%by passing 'last' as the THIRD argument, etc. SO BE CAREFUL OR FIX IT.

digits(p) %sets desried accuracy to p decimal places??

if length(varargin)>=1 %default param 
    param = varargin{1};
else
    param = 4;
end

out = [seed zeros(1,n-1)]; %preallocate

apply_fcn = vpa(seed); %see vpa doc??

for i=1:n-1
    apply_fcn = vpa(param.*apply_fcn.*(1-apply_fcn)); %note vpa again
    out(i+1) = apply_fcn;
end

if length(varargin)==2 && strcmp(varargin{2},'last')
    out = out(end);
end

Если я правильно понимаю, digits() Я должен делать свое дело, но я не уверен, и очень важно, чтобы у меня была точность, указанная при вызове функции (как это будет знать любой, кто знаком с хаосом!).

Заранее спасибо!