%input of system of ordinary differential equations (as a vector)
 
funcODE=[ output_prm(k,t)-c-ddelta*k;...       %ODE: k
          ((1/ssigma)*(mpk_prm(k,t)-ddelta-rrho))*(c-p*k)-m*p*k+p*(output_prm(k,t)-c-ddelta*k); %ODE: c
          m*(m-(mpk_prm(k,t)-ddelta)+(1/ssigma)*(mpk_prm(k,t)-ddelta-rrho))]; %ODE: m

         %first state, then control