Per il primo pezzo messo così funziona, bisogna ovviamente vedere se il risultato è quello che voleva otternere il tuo amico:

codice:
clear 
n=73; 

X=sqrt(3)*(randn(n,1)+1); 
z=((X-ones(n,1)).^2)./(2*3); 
B=1/(sqrt(2*3.142)*sqrt(3))*exp(-z); 
%plot(X,B,'ro') 
d=ones(n,1); 
h=0.5; 
z=0; 

u=[];
I=ones(n,1); 

for d = 1:10 
    z=z+1; 
    u(:,z)=(X-z)./h;
    fepa(z)=(1/n*h)*sum(3/4*(ones(n,1)-u(z).^2).*I) 
    plot(X,fepa(z),'b+') 
    %Premi un tasto per andare al grafico successivo
    pause
end; 

for i=1:n; 
    if abs(u(i))>1; 
        I(i)=0; 
    end; 
end;
Nell'ultimo pezzo senza altre info non posso proprio dirti nulla, ciao.