public double chi(double t, PhasePoint x)
{
double norm = Math.Sqrt(x
.Select(a => a * a)
.Sum());
if (norm > r)
{
return((norm / r) - 1);
}
else
{
return(0);
}
}
public double sigma(PhasePoint x)
{
return(Math.Sqrt(x
.Select(a => a * a)
.Sum()));
}