function [v,u]=Coulomb() v=1; u=0; end