急求卡方分布

function chichi
b=0:0.01:50;
c=chi2pdf(b,8);%自由度为8的卡方分布的密度函数值
d=chi2cdf(b,8);%卡房分布函数值
h=gampdf(b,9,0.5);%gamma分布的密度函数
g=gamcdf(b,9,0.5);%gamma分布函数
figure(1);clf
subplot(221);
plot(b,c,'linewidth',1.5)
xlim([0,30]);
ylim([0,max(c)+0.01]);
legend('卡方分布密度函数')
subplot(222);
plot(b,d,'r','linewidth',1.5)
xlim([0,30]);
ylim([0,1]);
legend('卡方分布的分布函数')
subplot(223);
plot(b,h,'g','linewidth',1.5)
xlim([0,10]);
ylim([0,0.3]);
legend('伽马分布密度函数')
subplot(224);
plot(b,g,'c','linewidth',1.5)
xlim([0,10]);
ylim([0,1]);
legend('伽马分布的分布函数')
end