งานเก่า MATLAB

MATLAB  23/09/55
MATLAB::fNewton
MAIN FUNCTION::

function [xn,fxn] = fnewton (fx,dfx,x0)

k = 1;
xn = x0 - fx(x0)/dfx(x0);
disp([k x0 xn fx(xn)]);

ep = 0.001;

maxiter = 100;
while abs(xn-x0) > ep & k <= maxiter
%while k <= maxiter
x0 = xn;
k = k+1;
xn = x0 - fx(x0)/dfx(x0);
disp([k x0 xn fx(xn)]);
end
fxn = fx(xn);

RUN IN NEXT THIS::
clear
clc
fx = @(x) x^6 - x -1;
dfx = @(x) 6*x^5 - 1;

x0 = 1.5;
[x,y] = fNewton(fx,dfx,x0);

ต่อ

การบ้าน clear
clc
K = 0.016;
Ca0 = 42;
Cb0 = 28;
Cc0 = 4;
fx = @(x) (Cc0+x)-K*((Ca0-2*x)^2)*(Cb0-x);
dfx = @(x) (2*(8*x - 168)*(x - 28))/125 + (2*(2*x - 42)^2)/125 + 1;

x0 = 5;
[x,y] = fnewton(fx,dfx,x0);


ต่อ

**MATLAB**

clear
clc
fx = @(x) x^6 - x -1;
dfx = @(x) 6*x^5 - 1;

x0 = 1.5;
k = 1;
xn = x0 - fx(x0)/dfx(x0);

disp([k x0 xn fx(xn)]);

ep = 0.001;
maxiter = 100;
while abs(xn-x0) > ep & k <= maxiter
%while k <= maxiter
x0 = xn;
k = k+1;
xn = x0 - fx(x0)/dfx(x0);
disp([k x0 xn fx(xn)]);
end

แปะของเก่า ^^