using System;
using System.Diagnostics;
using static System.Console;

public static class rndstart{

public static vector randomvec(vector a,vector b,Random rnd){
	var x=new vector(a.size);
	for(int i=0;i<a.size;i++)x[i]=a[i]+(b[i]-a[i])*rnd.NextDouble();
	return x;
	}

public static vector run
(Func<vector,double>f,vector a,vector b,int seconds=1,double acc=1e-3,int N=100)
{
Random rnd=new Random();
vector best=(a+b)/2;
double fbest=f(best);
var start_time=DateTime.Now;
do{
	vector x=randomvec(a,b,rnd);
	qnewton.min(f,ref x,acc:acc);	
	double fx=f(x);
	if(fx<fbest){
		best=x;
		fbest=fx;
		}
}while((DateTime.Now-start_time).TotalSeconds < seconds);
return best;
}//min

}//class