برنامه اصلی رو که ران میکنم خطای زیر رو میده :
??? Input argument "q" is undefined.



Error in ==> f at 3

    ret=-sum(alpha.*q.^2+k.*(D-sum(Q)).*q);