function P_v = ComputePrior(Vx,Vy,Prior_sigma)

P_v = exp(-0.5*(Vx.^2 +Vy.^2)/Prior_sigma^2);