function [p]=price()
global tau N alpha A L K w r;
alpha_= alpha.^(-alpha).*(1-alpha).^(alpha-1);
supply = alpha_.*A.*(L.^alpha).*(K.^(1-alpha));
q=(w.^alpha).*(r.^(1-alpha))./A;
share=consshare(supply,q);
p=sum(share.*(tau.*(q*ones(1,N))))';
return