エカつきのブログ = Eka tsuki no blog

clc;
clear;
x0=1.5;
y0=3.5;
disp('Solusi Persamaan Nirlanjar f(x,y)=x^2+xy-10=0');
disp(' f(x,y)=y+3xy^2-57=0');
disp('Dengan Metode Newton-Raphson');
disp('iterasi akar1 akar2');
for iterasi=0:50;
u=x0^2+x0*y0-10;
v=y0+3*x0*y0^2-57;
x=x0-((u*(1+6*x0*y0)+ v*x0)/((2*x0+y0)*(1+6*x0*y0)-x0*3*y0^2));
y=y0+((u*(3*y0^2)- v*(2*x0+y0))/((2*x0+y0)*(1+6*x0*y0)-x0*3*y0^2));
fprintf(' %3g %10.6f %10.6f\n', iterasi, x, y);
if (abs(x-x0)<0.000001) || (abs(y-y0)<0.000001);
 break;
end;
x0=x;
y0=y;
end;
fprintf('Jumlah iterasi = %g\n',iterasi);
fprintf('Akar akarnya adalah %10.6f dan %10.6f\n',x, y);