Non mi occupo di Matlab ma mi è arrivata una richiesta d'aiuto da qualcuno che ci deve studiare problemi di economia e .....cmq sotto posto il listato in cui viene segnalato errore alla riga in cui all'elemento dell'array u nel for verrebbe attribuito il valore. Premetto che nn conosco la sintassi di Matlab insomma fino a pochi minuti fa non l'avevo mai nemmeno visto :tongue: e se qualcuno è disposto ad aiutarmi anticipatamente lo ringrazio. Gli auguri di Buon Natale e Buon Anno li faccio indipendentemente da questo però
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;
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,'b+')
end;
I=ones(n,1);
for i=1:n;
if abs(u(i))>1;
I(i)=0;
end;
end;
ftri=(1/n*h)*sum(35/32*((ones(n,1)-u.^2).^3).*I)
fepa2=(1/n*1)*sum(3/4*(ones(n,1)-u.^2).*I)
fepa3=(1/n*2)*sum(3/4*(ones(n,1)-u.^2).*I)
ftri2=(1/n*1)*sum(35/32*((ones(n,1)-u.^2).^3).*I)
frti3=(1/n*2)*sum(35/32*((ones(n,1)-u.^2).^3).*I)
tvalue=1/(sqrt(2*3.142)*sqrt(3))*exp(-((0-1)^2)/(2*3))
sigma=std(X);
hopt=sigma/(n^0.2)