Double integral where the limit of y is a function of x
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
I am trying to work out double integral, where the limit of y is a function of x, and the limits itself are vectors, like
n=1:10;
t=2;
a=(n-1).*t;
b=n.*t;
fun=@(x,y)lognpdf(b-x,2,1).*lognpdf(y,3,2);
ymax=@(x) b-x;
Qab=@(a,b)integral2(fun,a,b,0,ymax);
Q=arrayfun(Qab,a,b);
but this returns error
Error using -
Matrix dimensions must agree.
Error in @(x)b-x
Is this because b is a vector?
1 Kommentar
Antworten (1)
Andrei Bobrov
am 10 Dez. 2013
Bearbeitet: Andrei Bobrov
am 11 Dez. 2013
EDIT
n=1:10;
t=2;
a=(n-1).*t;
b=n.*t;
[ai,y] = ndgrid(a,b);
bi = y.';
f = @(x,y,y2)logpdf(y2-x,2,1).*logpdf(y,3,2);
out = arrayfun(@(x1,y1,y2)integral2(@(x,y)f(x,y,y2),x1,y1,0,@(x)y2-x),ai,bi,y);
Siehe auch
Kategorien
Mehr zu Circuits and Systems finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!