figure
t=0:pi/50:pi;
t=0:pi/50:pi;
m= 5
linecolor = ['r';'b';'g';'k';'y';'r'];
a = [1/2,1,2,3,3.43954,5]
for ii=1:length(a)
y=a(ii)*m^(-a(ii)) *t.^(a(ii)-1).*exp(-t/m);
type = linecolor(ii);
plot(t,y,type);
hold on
end
legend('a=0.5','a=1','a=2','a=3','a=3.43954','a=5')