function [g_Gaussian,g_vonMises] = MotorPrimitives(P)
    mu = linspace(P.theta/P.nj,P.theta, P.nj);
    k = 1/((pi*P.sigma)/180)^2;
    theta = 1:P.theta;
    g_Gaussian = 1/(P.sigma*sqrt(2*pi)) * exp(-0.5*(theta-mu'./P.sigma).^2);
    g_vonMises = exp(k*cos((pi*(theta-mu'))/180))/(2*pi*besseli(0,k));
    