#include "real.h"
#include <math.h>

void rkstep(int n,real x,real y0[],real y1[]
	,real (*f)(real,real*,real*)
	,real h,real* err,real* w1,real* w2){

int i; real dy;

f(x    ,y0,w1); for(i=1,i<n,i++) y1[i]=y0[i]+.5*h*w1[i];
f(x+h/2,y1,w2); for(i=1,i<n,i++) y1[i]=y0[i]+   h*w2[i];
*err=0; for(i=1,i<n,i++){
	dy=fabs(h*(w1[i]-w2[i])/3);
	if(dy>*err) *err=dy;}
}