public static decimal A1(decimal t, decimal tay)
{
if (t != tay)
{
return(VectorHelper.Mult(
VectorHelper.Div(
Parametrization_ND.X0(tay),
Parametrization_ND.X0(t)),
Normal.GetNormal(t,
Parametrization_ND.X01d,
Parametrization_ND.X02d))
/ (decimal)Math.Pow(
(double)Normal.GetModul(
VectorHelper.Div(
Parametrization_ND.X0(t),
Parametrization_ND.X0(tay)))
, 2));
}
else
{
return(VectorHelper.Mult(
new[] { Parametrization_ND.X01d2(t), Parametrization_ND.X02d2(t) },
Normal.GetNormal(t,
Parametrization_ND.X01d,
Parametrization_ND.X02d))
/ (2M * (decimal)Math.Pow(
(double)Normal.GetModul(new[]
{
Parametrization_ND.X01d(t),
Parametrization_ND.X02d(t)
})
, 2)));
}
}
public static decimal H01(decimal t, decimal tay)
{
if (t != tay)
{
return((decimal)Math.Log(4 / Math.E *
Math.Pow(
Math.Sin(
(double)(t - tay) / 2
), 2)
/ Math.Pow((double)Normal.GetModul(
VectorHelper.Div(
Parametrization_ND.X0(t),
Parametrization_ND.X0(tay))), 2)) / 2);
}
else
{
return((decimal)Math.Log(1 /
(Math.E * Math.Pow((double)
Normal.GetModul(new[]
{
Parametrization_ND.X01d(t),
Parametrization_ND.X02d(t)
})
, 2))) / 2);
}
}
public static decimal H02(decimal t, decimal tay)
{
return((decimal)Math.Log(
(double)(1M / (
Normal.GetModul(
VectorHelper.Div(
Parametrization_ND.X0(t),
Parametrization_ND.X1(tay)))))));
}
public static decimal K11(decimal t, decimal tay)
{
return(VectorHelper.Mult(
VectorHelper.Div(
Parametrization_ND.X0(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.X0(tay)))
, 2));
}