Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- syms lamda;
- %
- % for i = 1:N
- % por(i) = lamda * exp( - lamda * data(i));
- % end
- % clear i;
- %
- % fun = prod(por);
- sumEl = sum(data);
- len = length(data);
- fun = lamda^len * exp(-lamda * sumEl);
- Lg = log(fun);
- dLg = diff(Lg, lamda)==0;
- figure;
- hold on;
- fplot(Lg, lamda);
- hold off;
- dLg = str2sym(char(dLg));
- sa = double(vpasolve(dLg, lamda));
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement