function [ori oriP oriN]=orivectors(n,d);

    %to be memorized oris (circle representation, "continous" model)
    ori=linspace(0,360,n+1);
    ori=ori(1:end-1);            % ori vector in degrees
    ori=ori+unique(diff(ori))/2; % rotate awyy from cardinal
    %(I think was by half of the inter-ori distance?)
    
    % test oris (CCW and CW)
    oriP=ori+d; 
    oriN=ori-d; 

    % make everything radians
    ori=deg2rad(ori);        
    oriP=deg2rad(oriP);        
    oriN=deg2rad(oriN); 
        
end