Scilab

[code]function y = M_W(t,D,b,n);
y = sum((1-cos(b.^n*t))./(b.^((2-D)*n)));
endfunction

D = 1.5;
n=[-100:100];
t = [0:0.1:500];
c = [];

for b = [1.4,1.5,1.6,1.7,1.8,1.9]
c = [c, feval(t,list(M_W,D,b,n))'];
end

plot2d(t',c,style = [1...

http://www.invest74.ru/index.php?topic=10122.msg28874#msg28874