F0（R0 , G0 , B0）＝（R × I , G × I , B × I）

f0（r0 , g0 , b0）＝（R0 ×（1 － S）, G0 ×（1 － S）, B0 ×（1 － S））

基準点と参照点の距離を Dt＝（（Xt － Xs）²＋（Yt － Ys）²＋（Zt － Zs）²）^1/2　とすると
Vd（Xd , Yd , Zd）＝（Xt － Xs , Yt － Ys , Zt － Zs）／ Dt

D ＝（（X － Xp）²＋（Y － Yp）²＋（Z －Zp）²）^1/2

Sh ＝1 ／（kc ＋ kl × D ＋ kq × D²）

Vp（Xp , Yp , Zp）＝（X － Xs , Y － Ys , Z － Zs）／ D

θ＝arccos（Xd × Xp ＋ Yd × Yp ＋ Zd × Zp）

F（RF , GF , BF）＝ 0.0

（θ≧ｂ）

F（RF , GF , BF）＝F0 × Sh ×（ｂ－θ）／（ｂ－ａ）

（ｂ＞θ＞ａ）

F（RF , GF , BF）＝F0 × Sh

（ａ≧θ）

Cosφ＝　－（Xp × Xv ＋ Yp × Yv ＋ Zp × Zv）
（Cosφ ＜ 0 の場合、Cosφ ＝ 0）

E（RE , GE , BE）＝ F × Cosφ