function da = DA( a, disa, eqeffort ) %Agents's disutility function des %eqeffort is one da=((a*eqeffort)^disa)*0.5; end