Displaying 1 result from an estimated 1 matches for "tukeysolve".
2007 Sep 24
1
Root finding problem
...o. :-( Again, what I want is to find the root of the
above-mentioned function, for arbitrary values of x, a, and lambda,
and just specifying one value for u as a starting guess. Could anyone
please give me a suggestion of how to do that? Thanks in advance.
In MATLAB I do the following:
function u=tukeysolve(a,lambda,xgrid)
u=zeros(length(xgrid),1);
old=0.001; % just a guess
for i=1:length(xgrid)
x=xgrid(i);
u(i) = fsolve(@tuki,old,optimset('Display','Off','tolf',5e-8),x,a,lambda);
old=u(i); % use previous value as new start value
end
function f=tuki(u,x,a,lambda)
f = x -...