MATLAB: Goel-Okumoto software reliability

software reliability

Hi everyone,

i have this code for the Goel-Okumoto software reliability model.

function [a,b]=go(N) k=length(N); l = 0; u = 3; n = sum(N); while ((u-l)>0.001) b=(l+u)/2; summa = 0; for i=2:k namn = exp(-b*(i-1))-exp(-b*i); summa = summa+N(i)*(i*exp(-b*i)-(i-1)*exp(-b*(i-1)))/namn; end; summa = summa+N(1)*(1*exp(-b))/(1-exp(-b*1)); summa = summa-n*k*exp(-b*k)/(1-exp(-b*k)); if (summa < 0) u=b; else l=b; end; end; a=n/(1-exp(-b*k));

I would like to know what are the values " l " and " u "????

Thanks in advance.

Best Answer

"u" appears to be "upper bound", and "l" appears to be "lower bound".

Related Solutions

MATLAB: I have to find the sum of the elements in the upper dimension matrix and this is the code i came up with. When I submitted this code it’s showing error.

Why not just:

[a,b]=size(A);
if a==b
    summa=sum(triu(A),'all');
end

MATLAB: Please help me in completing this program

3. There are no comments after the initialization comments. No comments about what the loops and ode15s are supposed to do.

4. No initialization of du early in the function so that an error in the program will throw an additional, secondary error about du not being assigned.

5. This code:

tau= 1:20;
for tau=1:20
    u(tau)=1;
end

can be replaced by this:

u = ones(1, 20);

6. Since u is 1 for any tau, this code probably does not do what you expect:

if(u(tau)>b/l)
        f(u(tau))=u(tau)-b/l;
    end
    if(u(tau)> -b/l && u(tau)<b/l)
        f(u(tau))=0;
    end
    if(u(tau)<-b/l)
        f(u(tau))=u(tau)+b/l;
    end

You're basically setting f(1) to something.

On the plus side though, you did explain the initial variables and eve gave their units. (Don't know what "pairin" means though.)

Best Answer

Related Solutions

MATLAB: I have to find the sum of the elements in the upper dimension matrix and this is the code i came up with. When I submitted this code it’s showing error.

MATLAB: Please help me in completing this program

Related Question