leasqr question

[Rainer]

Dear Colleagues

Below is a light change to the leasqrdemo to fit
my problem. But it doesnt work.
I get always the error:
"error: invalid column index = 3"
Has anybody an idea? Thanks

function leasqrtest
  % generate test data
  x = [1:0.1:100]';
  p = [1; 0.1; 150];
  data = p(1).*(p(3) - x).**(p(2));
  
%printf ("data = %f\n", data);

  F = @leasqrfunc;
  F
  
  wt1=0.1;
  
  dFdp = @leasqrdfdp; % exact derivative
  % dFdp = @dfdp;     % estimated derivative
  dp = [0.01; 0.01; 0.01];
  pin = [9; 1; 120]; 
  stol=0.001; niter=500;
  minstep = [0.0001; 0.0001; 0.001];
  maxstep = [0.5; 0.5; 0.5];
  options = [minstep, maxstep];

  global verbose;
  verbose=1;
  [f1, p1, kvg1, iter1, corp1, covp1, covr1,
stdresid1, Z1, r21] = ...
    leasqr (x, data, pin, F, stol, niter, wt1, dp,
dFdp, options);

function y = leasqrfunc(x,p)
    %y=p(1)*exp(-p(2)*x);
    p
    y=p(1).*(p(3) - x).**(p(2));
    
  endfunction

function y = leasqrdfdp(x,f,p,dp,func)
y=[(120-x).**(p(2)),p(1).*log(120-x).*(120-x).**(p(2))];
endfunction



      ________ 
Yahoo! Clever: Stellen Sie Fragen und finden Sie Antworten. Teilen Sie Ihr 
Wissen. www.yahoo.de/clever