echo off;
kuvano=kuvano+1;
figure(kuvano);
m  = 2*m ,[T,Y] = rk4('f',a,b,ya,m);points = [T;Y];

%pause  % Press any key to see the list of solution points.

%clc;, clg;
axis([a b c d]);...
plot(T,Y,'-g');... 
hold on;... 
plot([a b],[0 0],'b',[0 0],[c d],'b');...

if m<=32,
  plot(T,Y,'or');
end;...
xlabel('t');...
ylabel('y');...
title('Runge-Kutta solution to y` = f(t,y)');...
grid;...
axis;...
hold off;...
shg; pause  % Press any key to continue.

Mx1 = 'Runge-Kutta solution to y` = f(t,y).';
Mx2 = '     t(k)               y(k)';
clc,echo off,diary output,...
%disp(''),disp(Mx1),...
%disp(''),disp(Mx2),disp(points'),diary off,echo on
