double[] v2 = new double[nvar-n2];
for(int i=0;i<v2.length;i++)
v2[i] = v[n2+i];
h1=(x2-x1)/100.0;
y=buildVector(load1(x1,v));
Output out = new Output();
StepperDopr853 s1 = new StepperDopr853();
Odeint integ1 = new Odeint(y,x1,xf,atol,rtol,h1,hmin,out,d,s1);
integ1.integrate();
f1=buildVector(score(xf,y));
y=buildVector(load2(x2,v2));