xp[x_, y_] = 2  x - (x^2) - x*y;
yp[x_, y_] = y - (y^2) - 2*x*y;
x0 = 1;
y0 = 1;
t0 = 0;
tn = 1;
h =  .2;

f[t_, x_, y_] = xp[x, y];
g[t_, x_, y_] = yp[x, y];;
\(\(\(xi = x0\)\(;\)\(\ \)\(yi = y0\)\(;\)\(
\)\)\)
\(\(For[ti = t0, \ ti ≤ \((tn - h)\), \ ti += h, \ 
	
        k1 = h*f[ti, xi, yi]; 
	m1 = h*g[ti, xi, yi]; 
	
        k2 = h*f[ti + h\/2, xi + k1\/2, yi + m1\/2]; 
	
        m2 = h*g[ti + h\/2, xi + k1\/2, yi + m1\/2]; 
	
        k3 = h*f[ti + h\/2, xi + k2\/2, yi + m2\/2]; 
	
        m3 = h*g[ti + h\/2, xi + k2\/2, yi + m2\/2]; 
	
        k4 = h*f[ti + h, xi + k3, yi + m3]; 
	
        m4 = h*g[ti + h, xi + k3, yi + m3]; 
	
        K = \(k1 + 2  k2 + 2  k3 + k4\)\/6; 
	
        M = \(m1 + 2  m2 + 2  m3 + m4\)\/6; 
	xi += K; 
	
        yi += M;];\)\)
Print[xi];
Print[yi];