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

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

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

a = [A / 36525]
b = A - 36525 * a
c = [(b + 1) / 365.25]
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
Y < 1のとき、紀元前となるので、B.C.は1-Y年となります。
解が12月32日の場合は、翌年の1月1日としてください。

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

ユリウス通日 :