function y=pred_cor3(a,b,n,y0,f)

  h=(b-a)/n;
  x(1:n+1)=a+h*[0:n];
  y=RK44(a,b,n,y0,f);

for i=3:n

  aux1=h*feval(f,x(i),y(i));
  aux2=h*feval(f,x(i-1),y(i-1));
  aux3=h*feval(f,x(i-2),y(i-2));

  y_pr=y(i)+(1/12)*(23*aux1-16*aux2+5*aux3);
  y_cor=y(i)+(1/12)*(h*5*feval(f,x(i+1),y_pr)+8*aux1-aux2);
  y(i+1)=y_cor;

endfor


endfunction