수치해석 소스

본문내용

c
% xl, xu = lower and upper guesses
% es = (optional) stopping criterion (%) (default = 0.001)
% maxit = (optional) maximum allowable iterations (default = 50)
if func(xl)*func(xu)>0 %if guesses do not bracket a sign change
error(\'no bracketssss\') %display an error message and terminate
end
% default values
if nargin<5, maxit=50; end
if nargin<4, es=0.001; end
% false position
iter = 0;
xr = xl;
while (1)
xrold = xr;
xr = xu - func(xu)*(xl - xu)/(func(xl) - func(xu));
iter = iter + 1;
if xr ~= 0, ea = abs((xr - xrold)/xr) * 100; end
test = func(xl)*func(xr);
if test < 0
xu = xr;
elseif test > 0
xl = xr;
else
ea = 0;
end
if ea <= es | iter >= maxit, break, end
end
root = xr;
falsepos(inline(\'(1.5.*x)./((1+x.^2)).^2-0.65.*atan(1./x)+(0.65.*x)./(1+x.^2)\'),0.48,0.5 )
ans =
0.48086448538221

키워드

수치해석 소스,   수치,   수치해석