function [t,x]=bi4(f,x0,ti,tf,h)
//metodo de BRK4 con iteracion de Newton
//exec('jacobiant.sci')
//exec('rk4.sci')
t=[ti:h:tf];
x=zeros(length(x0),length(t));
x(:,1)=x0;
k=0;
I=eye(length(x0),length(x0))//matriz identidad
for tk=ti:h:tf-h
  k=k+1;
  xk12=rk4_step(f,x(:,k),t(k),h/2);
  x(:,k+1)=brk4_step(f,xk12,t(k)+h/2,h/2);
end
endfunction


