The Math Forum

Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Math Forum » Discussions » Software » comp.soft-sys.matlab

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: Error : Inner matrix dimensions must agree When computing the three dimensional integral.
Replies: 7   Last Post: Jun 24, 2013 11:14 AM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]

Posts: 5
Registered: 6/19/13
Re: Error : Inner matrix dimensions must agree When computing the three dimensional integral.
Posted: Jun 20, 2013 12:47 PM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply


I do appreciate your reply!
You are right. I also checked the triplequad function detailedly yesterday. And as I defined, z should be a vector and then I can use the function triplequad.

So as you said, the two matrixs don't work. Then how should I define the f function?
The exponentional part of the f function is just the part of the three dimension
joint Gaussian distributation as following:

exp( ([(x-m1),(y-m2),(z-m3)]*(cov1)*[(x-m1);(y-m2);(z-m3)])./(-2) )./((sqrt(2*pi)).^3.*(sqrt(det(cov))))

It has three matrixs operation. So how should I define?
Thanks for your nice help!!!

"Steven_Lord" <> wrote in message <kpv4sv$otv$>...
> "The first input, fun, is a function handle. fun(x,y,z) must accept a vector
> x and scalars y and z, and return a vector of values of the integrand."
> When your function is called with a vector z and scalar x and y, the
> expressions (x-m1) and (y-m2) will be scalars but the expression (z-m3) will
> be a vector. One of the concatenation steps will not work:
> [(x-m1),(y-m2),(z-m3)] % Doesn't work if z is a column vector
> [(x-m1);(y-m2);(z-m3)] % Doesn't work if z is a row vector
> I believe z will be a row vector, so the latter step will not work.
> --
> Steve Lord
> To contact Technical Support use the Contact Us link on

Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© The Math Forum at NCTM 1994-2018. All Rights Reserved.