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

clc;
clear;
x0=1.5;
y0=3.5;
epsilon=0.000001;
disp('Solusi Persamaan Nirlanjar f(x,y)=x^2+xy-10=0');
disp(' f(x,y)=y+3xy^2-57=0');
disp('Dengan Metode Titik Tetap');
disp('r x y |x(r+1)-x(r)| |y(r+1)-x\y(r)|');
for iterasi=1:50;
x=10/(x0+y0);
y=50/(1+3*x*y0);
fprintf('%3g %12.7f %12.7f %12.7f %12.7f\n', iterasi, x, y, abs(x-x0),abs(y-y0));
if (abs(x-x0) break;
end;
x0=x;
y0=y;
end;
akar1=x;
akar2=y;
fprintf('Akar 1 adalah %10.7f\n',x);
fprintf('Akar 2 adalah %10.7f\n',y);
fprintf('Jumlah iterasi = %3g\n', iterasi);