t = 0;
h = 0.2;
y = 0;

for k = 1: 5
k1 = tarea(t, y);
k2 = tarea(t+h/2, y+k1*h/2);
k3 = tarea(t+h/2, y+k2*h/2);
k4 = tarea(t+h, y+k3*h);
disp([k1,k2,k3,k4]);
t = t+h
y = y + (k1+2*k2+ 2*k3+ k4)*h/6
end