function [t,y]=javEE_zh22(diffegy,t0,y0,T,N) 
t=zeros(N+1,1); 
y=zeros(N+1,1);
t(1)=t0; 
y(1)=y0; 
h=(T-t0)/N;
for i=1:N
    t(i+1) = t(i) + h;
    k1 = diffegy(t(i),y(i));
    k2 = diffegy(t(i)+h,y(i)+h*k1);
    y(i+1) = y(i) + h * (k1+k2)/2; 
end