l22

clear all;
ti = [10 10];
for j = 1:length(ti)
    h = 0.1;
    t = 0:h:100;
    N = length(t);
    x = t;
    x(1) = ti(j);
    for i=2:N
        k1 = h*(abs(sin(t(i-1)))*x(i-1)-x(i-1)*x(i-1));
        k2 = h*(abs(sin(t(i-1)+1/2*h))*(x(i-1)+1/2*k1)-(x(i-1)+1/2*k1)*(x(i-1)+1/2*k1));
        k3 = h*(abs(sin(t(i-1)+1/2*h))*(x(i-1)+1/2*k2)-(x(i-1)+1/2*k2)*(x(i-1)+1/2*k2));
        k4 = h*(abs(sin(t(i-1)+h))*(x(i-1)+k3)-(x(i-1)+k3)*(x(i-1)+k3));
       x(i) = x(i-1)+1/6*(k1+2*k2+2*k3+k4);
    end
    plot(t,x,'g')
    hold on
end