Function EstimateTideValues() As [Function]
Dim Pi As [Double] = Math.PI
Dim sum As [Double] = harmonics(0, 1)
For k As Int16 = 1 To harmonics.GetUpperBound(0)
Dim period As [Double] = 2 * Pi / harmonics(k, 0)
Dim coefCos As [Double] = harmonics(k, 1)
Dim coefSin As [Double] = harmonics(k, 2)
sum = sum + coefCos * Math.Cos(time * period)
sum = sum + coefSin * Math.Sin(time * period)
Next
Return sum
End Function
double EstimateTideValues()
{
double Pi = Math.PI;
double sum = harmonics[0, 1];
for (int k=1; k <= harmonics.GetUpperBound(0); k++)
{
double period = 2 * Pi / harmonics[k, 0];
double coefCos = harmonics[k, 1];
double coefSin = harmonics[k, 2];
sum += coefCos * Math.Cos(time * period);
sum += coefSin * Math.Sin(time * period);
}
return sum;
}