y0=[1;1;0.5;1;1];
dt=0.0001;t=0;

y=y0;
for n=1:10000
t=t+dt;
f=[1 - t*y(1), exp(-t) - 3*y(2), t^2-y(3), -y(4)*abs(y(4)), sin(t)-sqrt(abs(y(5)))];
k1=dt*f;
yd=y+k1';
td=t+dt;
k2=dt*[1 - td*y(1), exp(-td) - 3*yd(2), t^2-yd(3), -yd(4)*abs(yd(4)), sin(td)-sqrt(abs(yd(5)))];
y=y+dt*f'; 
if mod(n,1000)==0 [t,y'],end
end



ans =

  Columns 1 through 5 

    0.1000    1.0947    0.8228    0.4527    0.9091    0.9074
    0.2000    1.1775    0.6837    0.4119    0.8333    0.8292
    0.3000    1.2471    0.5736    0.3788    0.7692    0.7648
    0.4000    1.3024    0.4857    0.3545    0.7143    0.7131
    0.5000    1.3428    0.4148    0.3402    0.6666    0.6734
    0.6000    1.3682    0.3570    0.3368    0.6250    0.6446
    0.7000    1.3788    0.3095    0.3451    0.5882    0.6254
    0.8000    1.3754    0.2700    0.3660    0.5555   0.6148
    0.9000    1.3591    0.2369    0.4002    0.5263   0.6117
    1.0000    1.3313    0.2088    0.4482    0.5000   0.6147