matlab - ガンマ分布の中央値を計算するにはどうすればよいですか?

Question

以下の分布があります。これの中央値を計算します。

x=0:0.01:10;
x=[x' x' x' x' x'];
a=ones(1,1001)';
a=[a*2 a*4 a*6 a*8 a*10];
b=2;
f = gampdf(x,a,b);
plot(x,f)
grid on