sum of series with vectors

Hi,

I'm trying to solve power evolution using this equation:

I have all the values but can't get the summations to work. Thanks in advance!

Rhoo = 300e-5;
Caplambda = 1e-4;
Beta = [24 123 117 262 108 45];
lambda = [0.0127 0.031 0.116 0.3106 1.4006 3.876];
t = [0:1:30];
omega = [-6789970.6255578384124380944468818, -3.6511008223873248946050418708074, -1.1912661084218293186697922029422, -0.19487209040234811114948715764004, -0.068975792970862205151653629504031, -0.015127732106203404230259546907564, 0.00000038470100602825311618155332751786];

Answer 1

[Chunru]

Rhoo = 300e-5;
Caplambda = 1e-4;
Beta = [24 123 117 262 108 45];
lambda = [0.0127 0.031 0.116 0.3106 1.4006 3.876];
t = [0:1:30];
omega = [-6789970.6255578384124380944468818, -3.6511008223873248946050418708074, -1.1912661084218293186697922029422, -0.19487209040234811114948715764004, -0.068975792970862205151653629504031, -0.015127732106203404230259546907564, 0.00000038470100602825311618155332751786];
nn0 = zeros(size(t));
for it = 1:numel(t)
    for j=1:6       % check your formular which has j=0
        s1 = sum(Beta.*lambda./(omega(j)+lambda).^2);
        nn0(it) = nn0(it) + exp(omega(j)*t(it))/(omega(j)*(Caplambda+s1));
    end
end
nn0 = nn0/Rhoo;