Here is an example. You can use univariate procedure for any quantile value.

data t1;
   call streaminit(123);
   do i=1 to 100;
      x=rand('chisquare',2);
          output;
        end;
run;

proc means data=t1  noprint p25 p50 p75;
var x;
output out=q p25=x_p25 p50=x_p50 p75=x_p75;
run;

data _null_;
  set q;
  call symputx( 'x_p25',x_p25);
  call symputx( 'x_p50',x_p50);
  call symputx( 'x_p75',x_p75);
run;

%put >>>>p25=&x_p25 p50=&x_p50 p75=&x_p75<<<<;