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;
                        }