設圓心坐標（a,b)
因為圓C與y軸相切所以r=a
根據幾何關系得d^2=r^2-(√7)^2
根據點到直線距離公式得d^2=[(a-b)^2]/2
所以r^2-(√7)^2=[(a-b)^2]/2①
又因為圓心C在直線l:x-3y=0上
所以a=3b②
聯立①②解得
a=±3 b=±1
所以圓C的方程為（x±3)^2+(y±1)^2=9