%%diferansiyel denklem çözelim%%%
clear
%%ilk değer%%
x1(1)=0.5;
x2(1)=0.5;
x3(1)=0.5;
x_v(:,1)=[x1(1); x2(1); x3(1)];
%%parametreler%%
alpha11=0;
alpha12=-5;
alpha13=0;
alpha21=1;
alpha22=-2;
alpha23=0;
alpha31=0;
alpha32=0;
alpha33=-1;
a_matris=[alpha11 alpha12 alpha13; 
          alpha21 alpha22 alpha23; 
          alpha31 alpha32 alpha33];
%%bazı büyüklükler%%
iterasyon=5000;
mu=0.01;
%%dif denklem%
for k=1:iterasyon
    x_v(:,k+1)=x_v(:,k)+mu*(a_matris*x_v(:,k));
end
subplot(4,1,1), plot(x_v(1,:))
subplot(4,1,2), plot(x_v(2,:))
subplot(4,1,3), plot(x_v(3,:))
subplot(4,1,4), plot3(x_v(1,:),x_v(2,:),x_v(3,:)),grid on,