[matlab] programmer recurrence

bonjour,
je dois programmer une relation de recurrence avec matlab afin de trouver alpha:
for k=2:6
alpha(k+1) = alpha(k)-(alpha(k)-alpha(k-1))*h(alpha(k))/(h(alpha(k))-h(al­pha(k-1))) ;
end

ou h est la fonction definie par:
h=1+psi(1/alpha)/alpha-sum(V(i))^(alpha)*log(sum(V(i))/sum(V­(i))^alpha+log(alpha/N*sum(V(i))^(alpha)/alpha))

mais je ne sais pas comment rentrer cette fonction, vu qu'on me dit qu'il ne connait pas alpha...ce que justement je veux calculer...

pourriez vous m'aider?
merci

Configuration: Windows XP
Internet Explorer 7.0