using System;
using static System.Math;
using static System.Console;
public class network{

public vector p;
public Func<double,double> f;
public readonly int n;

public network(int n, Func<double,double> f){
	this.n=n;
	this.f=f;
	this.p=new vector(3*n);
	}

public double feed(double x){
	double sum=0;
	for(int i=0;i<n;i++){
		double a=p[3*i+0];
		double b=p[3*i+1];
		double w=p[3*i+2];
		sum+=w*f((x-a)/b);
		}
	return sum;
	}

public void train(double[] xs,double[] ys){
	int ncalls=0;
	Func<vector,double> mismatch=(u)=>{
		ncalls++;
		p=u;
		double sum=0;
		for(int k=0;k<xs.Length;k++)
			sum+=Pow(feed(xs[k])-ys[k],2);
Error.Write($"sum={sum}\n");
		return sum/xs.Length;
		};
	vector v=p;
	//v=qnewton.sr1(mismatch,v,1e-3);
	v=simplex.downhill(mismatch,v,step:0.5,dx:1e-5);
	p=v;
Error.Write($"ncalls={ncalls}\n");
	}

}//network