まずはユリウス通日を代入してAとBの値を求めます。

A = [ユリウス通日 + 68569.5]
B = ユリウス通日 + 0.5

AとBを用いて下記公式でユリウス通日からグレゴリオ暦を求めることができます。

a = [A / 36524.25]
b = A - [36524.25 * a + 0.75]
c = [(b + 1) / 365.25025]
d = b - [365.25 * c] + 31
e = [d / 30.59]
f = [e / 11]

u = 100 * (a - 49) + c + f
v = e - 12 * f + 2
w = d - [30.59 * e] + B - [B]
x = (w - [w]) * 24
y = (x - [x]) * 60
z = (y - [y]) * 60

グレゴリオ暦 (世界時) の年Y、月M、日D、時h、分m、秒sは、

Y = u
M = v
D = [w]
h = [x]
m = [y]
s = [z]

[X]は、実数Xに対してX以下の最大の整数です。
例1 : [135.79] = 135
例2 : [-135.79] = -136
解が12月32日の場合は、翌年の1月1日としてください。

ユリウス通日からグレゴリオ暦を求めるプログラム

ユリウス通日 :