public static decimal K21(decimal t, decimal tay)
{
if (t != tay)
{
return(VectorHelper.Mult(
VectorHelper.Div(
Parametrization_ND.X1(tay),
Parametrization_ND.X1(t)),
Normal.GetNormal(t,
Parametrization_ND.X11d,
Parametrization_ND.X12d))
/ (decimal)Math.Pow(
(double)Normal.GetModul(
VectorHelper.Div(
Parametrization_ND.X1(t),
Parametrization_ND.X1(tay)))
, 2));
}
else
{
return(VectorHelper.Mult(
new[] { Parametrization_ND.X11d2(t), Parametrization_ND.X12d2(t) },
Normal.GetNormal(t,
Parametrization_ND.X11d,
Parametrization_ND.X12d))
/ (2M * (decimal)Math.Pow(
(double)Normal.GetModul(new[]
{
Parametrization_ND.X11d(t),
Parametrization_ND.X12d(t)
})
, 2)));
}
}