clear
x0=[1;0];
tspan=[0 0.5];
[t,x]=ode45('keadaan',tspan,x0);
plot(t,x(:,1)),grid on
pause
plot(t,x(:,2)),grid on
