Перевод из Visual Basic в С++

Sub ВолновоеУравнениеЯвнаяСхема()
Dim x() As Double, y() As Double, ynp1() As Double, ynm1() As Double
Dim h As Double, tau As Double, r As Double, a As Double
Dim hPr As Double, tauPr As Double, Tmax As Double, L As Double
Dim nn As Long, mm As Long, MPr As Long, a1 As Double, a2 As Double
Dim n As Long, m As Long, ngr As Integer
Dim t As Double, rr As Double, pi As Double
pi = Atn(1) * 4
’ Исходные данные для задачи
L = 1
a = 1 ’ Это скорость распространения волны
Tmax = 0.5 ’ Верхняя граница временной координаты
h = 0.01 ’ Это шаг по пространственной координате
r = 1 ’Параметр сетки
tau = r * h ’Это шаг по временной координате
hPr = 0.1 ’Шаг для вывода значений сеточной функции по координате x
tauPr = 0.1 ’ Это шаг вывода графиков по временной координате
’NPr = CLng(Tmax / tauPr) + 1 ’Число графиков
MPr= CLng(hPr/h) ’Через столько точек выводить результаты по оси Ox
mm = CLng(L/h) ’ Количество ячеек сетки по оси Ox
nn = CLng(Tmax/tau) ’Количество ячеек сетки по оси Ot
ReDim x(mm), y(mm), ynp1(mm), ynm1(mm) ’Размерность массивов
For m = 0 To mm
x(m) = h * m
’Подстановка начальных условий
’Нулевой временной слой
ynm1(m) = Sin(pi * x(m))
’Первый временной слой
y(m) = ynm1(m) + tau * 0.5 * Cos(0.5 * pi * x(m))
Next m
Dim NameSh As String ’Имя листа для вывода результатов и графиков
NameSh = "ВолновоеУравнение"
Sheets(NameSh).Select
Cells.Clear ’Очистить рабочий лист
Cells(1, 1) = "t\x"
Cells(2, 1) = "t=0"
ngr = 1 ’Номер текущего графика
For m = 0 To mm Step MPr
Cells(1, m \ MPr + 2) = CDbl(Format(x(m), "#0.00"))
Cells(ngr + 1, m \ MPr + 2) = CDbl(Format(ynm1(m), "#0.00"))
Cells(1, m \ MPr + 2).Interior.Color = RGB(220, 90, 220)
Next m
a1 = 2 * (1 - r ^ 2)
a2 = r ^ 2
For n = 2 To nn
t = tau * n
For m = 1 To mm - 1
ynp1(m) = a1 * y(m) + a2 * (y(m + 1) + y(m - 1)) - ynm1(m)
Next m
Call ГраничныеУсловияВолновое(ynp1, t, h, mm)
’Переставить активные временные слои для следующей итерации
For m = 0 To mm
ynm1(m) = y(m)
y(m) = ynp1(m)
Next m
’Записать решение на данном временном слое а
’Записываются с шагом tauPr
If Abs((ngr) * tauPr - t) < tau / 10 Then
ngr = ngr + 1
Cells(ngr + 1, 1) = "t=" + Format(t, "#0.00")
For m = 0 To mm Step MPr
Cells(ngr + 1, m \ MPr + 2) = CDbl(Format(y(m), "#0.00"))
Next m
End If
Next n
Range("a1:a" + Trim(Str(ngr + 1))).Interior.Color = _
RGB(220, 90, 220)
Dim NameDiagr As String ’Название диаграммы
NameDiagr = "Волновое уравнение"
Call ПостроитьГрафики(ngr, CInt(mm / MPr), NameSh, _
NameDiagr, 0, 1, 0, 1)
End Sub
Sub ГраничныеУсловияВолновое(ynp1() As Double, _
t As Double, h As Double, mm As Long)
Dim pi As Double
pi = Atn(1) * 4
ynp1(0) = 0 ’Левая граница
’Правая граница
ynp1(mm) = (ynp1(mm - 1) + h * _
0.1 * Sin(2 * pi * t)) / (1 + h)
End Sub