ainitial = 0;
arange=[0 2];
[a,A] = ode45(@rkfunc, arange, ainitial);

binitial = A(end); //This is where I know I need a better assignment
brange=[2 4];
[b,B] = ode45(@rkfunc, brange, binitial);