Hallo,
im Zweifelsfall per for-Schleife:
L = ones(size(X1));
for I=1:numel(t)
L = L .* (X1./X2.*(t(I)./X2).^(X1-1).*exp(-(t(I)./X2).^X1));
end
Falls das zu lange rechnet, müsste man sich noch genauer ansehen, ob das beschleunigt werden kann.
Grüße,
Harald
↧