MATLAB 求和CN=C1*e^(-N)+e^(-N)*∑(e^(CN-1-n)其中,sigma的下标n=0,上标N-1;C1是已知的,求N=50时的CN值?（N是下标）.大家帮我看看怎么编写这个程序.

MATLAB 求和CN=C1*e^(-N)+e^(-N)*∑(e^(CN-1-n)其中,sigma的下标n=0,上标N-1;C1是已知的,求N=50时的CN值?（N是下标）.大家帮我看看怎么编写这个程序.
MATLAB 求和
CN=C1*e^(-N)+e^(-N)*∑(e^(CN-1-n)
其中,sigma的下标n=0,上标N-1;C1是已知的,求N=50时的CN值?（N是下标）.
大家帮我看看怎么编写这个程序.

MATLAB 求和CN=C1*e^(-N)+e^(-N)*∑(e^(CN-1-n)其中,sigma的下标n=0,上标N-1;C1是已知的,求N=50时的CN值?（N是下标）.大家帮我看看怎么编写这个程序.
%by dynamic
%2009.2.17
%
%这函数如果使用解析解法,估计Matlab无能为力
%下面我使用数值解法帮你求解
%
C1=10;
N=50;
fun=@(CN)C1*exp(-N)+exp(-N)*sum(exp(CN-1:CN-1-(N-1)))-CN;
[x,fval]=fsolve(fun,rand)
Optimization terminated:first-order optimality is less than options.TolFun.
x =
0
fval =
1.9287e-021