Matlab - Leyenda al representar una matriz de vectores

close all, clear all, clc

N=[5 10 15 20 25];

gamma_a=pi*8/10;

t = 2*(gamma_a^2)/(1+2*N(1));

[Mn1, rho] = MnRho(N(1),t);

M = Mn1;

for n=2 : length (N)

 t = 2*(gamma_a^2)/(1+2*N(n));

[Mn, rho] = MnRho(N(n),t);

M= [M;Mn];

end

figure, plot (rho, M), title ('Mn para N=[5 10 15 20 25]')

function [ Mn, rho ] = MnRho( N , t)

%function [ Mn, rho ] = MnRho( N , t)

    rho = 0:0.001:pi;

    f1=exp(-(rho.^2)/t);

    f2=zeros(size(rho));

    for k=0:N

        f2=f2+(((rho.^2)/t).^k)/factorial(k);

    end

    Mn=f1.*f2;

end