function [r0,e0,k,mass]=CO()
  r0=au(1.13,"angstrom");
  e0=au(1.076e6,"joule/mol");
  k=au(1870,"newton/m");
  mass=au(1.0,"carbon")+au(1.0,"oxygen");
endfunction
