clear;
clf;
%Mencari x(1) dengan diketahui
%Persamaan differensial dx(t)/dt= f(x(t),t)
% yang ada pada file dfungsi.m
% dan keadaan awal tertentu
k(1)=1;
% Dari Keadaan awal pembagian 1 interval:
x(1)=1+1*dfungsi(1,0);
for n=2:1000
   dt=1/n;
   xd(1)=1;
   td(1)=0;
   for i=1:n
      td(i+1)=td(i)+dt;
      xd(i+1)=xd(i)+dt*dfungsi(xd(i),td(i));
   end
   k(n)=n;
   x(n)=xd(n+1);
   if (abs(x(n)-x(n-1)))<0.000001
    break
   end
end
N=k
X=x
plot(k,x), grid on
   