function singularity = q3_A(x,n,a) if(x>a) singularity = (x-a)^n; else singularity = 0; end