Path: news.mathworks.com!not-for-mail
From: "David Doria" <daviddoria@gmail.com>
Newsgroups: comp.soft-sys.matlab
Subject: problem with subs
Date: Mon, 25 Feb 2008 22:08:05 +0000 (UTC)
Organization: The MathWorks, Inc.
Lines: 39
Message-ID: <fpve85$i5d$1@fred.mathworks.com>
Reply-To: "David Doria" <daviddoria@gmail.com>
NNTP-Posting-Host: webapp-03-blr.mathworks.com
Content-Type: text/plain; charset="ISO-8859-1"
Content-Transfer-Encoding: 8bit
X-Trace: fred.mathworks.com 1203977285 18605 172.30.248.38 (25 Feb 2008 22:08:05 GMT)
X-Complaints-To: news@mathworks.com
NNTP-Posting-Date: Mon, 25 Feb 2008 22:08:05 +0000 (UTC)
X-Newsreader: MATLAB Central Newsreader 1105197
Xref: news.mathworks.com comp.soft-sys.matlab:453584



I make a big 4x4 matrix of second partials:

------

syms A N sigma omega phi x t positive
lnf = -N/2 * log(2*pi*sigma) + -1/(2*sigma)*(x -
A*cos(omega*t - phi));
F=[diff(diff(lnf,A),A) diff(diff(lnf,A),omega)
diff(diff(lnf,A),phi) diff(diff(lnf,A),sigma);
	diff(diff(lnf,omega),A) diff(diff(lnf,omega),omega)
diff(diff(lnf,omega),phi) diff(diff(lnf,omega),sigma);
	diff(diff(lnf,phi),A) diff(diff(lnf,phi),omega)
diff(diff(lnf,phi),phi) diff(diff(lnf,phi),sigma);
	diff(diff(lnf,sigma),A) diff(diff(lnf,sigma),omega)
diff(diff(lnf,sigma),phi) diff(diff(lnf,sigma),sigma)];

----------

Then I want to substitute some values:
subs(F, {A,omega, phi, sigma}, {1, 1, 1, 1})

it doesn't work.

I tried
subs(F(1,2), {A,omega, phi, sigma}, {1, 1, 1, 1})

still no success.

However, if I do F(1,2) then manually copy and paste:
subs(1/2/sigma*sin(-omega*t+phi)*t, {A,omega, phi, sigma},
{1, 1, 1, 1})

it works properly.

Is there a way to do this automatically?

Thanks,

Dave