% ed.m: excess demand function

   function edv=ed(p);

     edv=5*(1000^(.2))*(p^(-.5)) -11-5*p;
