clear;
clf;
format long;
% Mencari solusi x(t), 0<t<ta, dengan diketahui:
% Persamaan differensial biasa (ODE) dx(t)/dt= f(x(t),t)
% yang ada pada file dfungsi.m dan keadaan awal tertentu x(0)
% Dari Keadaan awal pembagian 1 interval:
x0=1;
ta=10;
x(1)=x0;
alpha=1;
beta=alpha;
w2=1/(2*alpha);
w1=1-w2;
F1=ta*dfungsi(x(1),0);
F2=ta*dfungsi((x(1)+(beta*F1)),(alpha*ta));
x(2)=x(1)+(w1*F1)+(w2*F2);
xta(1)=x(2);
for n=2:500
   dt=ta/n;
   x(1)=x0;
   t(1)=0;
   for i=1:n
      t(i+1)=t(i)+dt;
      F1=dt*dfungsi(x(i),0);
		F2=dt*dfungsi((x(i)+(beta*F1)),(alpha*dt));
      x(i+1)=x(i)+(w1*F1)+(w2*F2);
   end
   clc
   k=n
   xta(n)=x(n+1);
   if (abs(xta(n)-xta(n-1)))<0.000001
    break
   end
end
plot(t,x), grid on