Skip to content

Instantly share code, notes, and snippets.

View EzAccount's full-sized avatar
💭
(non) Productive mood

Mikhail Tikhonov EzAccount

💭
(non) Productive mood
View GitHub Profile
tgA[i] = 1 / n_curved(x_temp, i*dz + dz/2) * (n_x_curved(x_temp) * dz - )
tg[j, i] = tgA[i] + tg[j, i - 1]
x_temp = x[i - 1] + dz / 2. * tg[j, i]
x[i] = x[i - 1] + dz * tg[j, i]