function [B,V,MSE,Y]=TPTest2025(m,s,mc)
    n=[10 50 100 250 500 1000 2000];
    N=length (n);
    for i=1:N
        b=[];
        v=[];
        for j=1:mc
            X=random('norm',m,s,n(i),mc);
            m1=mean(X);
            m2=sum(X)/(n(i)+1);
            m3=(max(X)+min(X))/2;
            b(j,:)=[mean(m1) mean(m2) mean(m3)]-m;
            v(j,:)=[var(m1) var(m2) var(m3)];
        end
        B(i,:)=mean(b);
        V(i,:)=mean(v);
    end
   MSE =B.^2+V;
   
   plot(n,B)
   title('Variation du biais en fonction de la taille de l''echantillon')
   legend('\hat{\mu}_1','\hat{\mu}_2','\hat{\mu}_3')
   figure
   plot(n,V)
   title('Variation de la variance en fonction de la taille de l''echantillon')
   legend('\hat{\mu}_1','\hat{\mu}_2','\hat{\mu}_3')
   
   figure
   plot(n,MSE)
   title('Variation du MSE en fonction de la taille de l''echantillon')
   legend('\hat{\mu}_1','\hat{\mu}_2','\hat{\mu}_3')
   for j=1:mc
            X=random('norm',m,2,1000,mc);
            m1=mean(X);
            m2=sum(X)/(n(i)+1);
            m3=(max(X)+min(X))/2;
            b(j,:)=[mean(m1) mean(m2) mean(m3)]-m;
            v(j,:)=[var(m1) var(m2) var(m3)];
   end
   Y=b./sqrt(v);
   figure
   histfit(Y(:,1))
   title('distribution du \mu_1')
    figure
   histfit(Y(:,2))
   title('distribution du \mu_2')
    figure
   histfit(Y(:,3))
   title('distribution du \mu_3')
end