dharr

Dr. David Harrington

9294 Reputation

22 Badges

21 years, 324 days
University of Victoria
Professor or university staff
Victoria, British Columbia, Canada

Social Networks and Content at Maplesoft.com

Maple Application Center
I am a retired professor of chemistry at the University of Victoria, BC, Canada. My research areas are electrochemistry and surface science. I have been a user of Maple since about 1990.

MaplePrimes Activity


These are replies submitted by dharr

@bashar27 Your explanation was in the title and cut off; we just got "mainly i need to collect independent term after...". So please upload your worksheet and put your explanation in the body of the post.

I think the point about a "time" variable is that the the associated conditions are "initial" conditions, so if it isn't actually time it can be treated as time if the conditions only specify one end of its domain, in contrast to a spatial variable where boundary conditions are specified at both ends of the domain. See the description of the "time" option. "These options are typically needed only for first order non-periodic problems where only one of the end points appears in the condition." So it looks to see if there are condition at one or both ends of the domain.

If you specify your rod problem in more detail it might be easier to answer, but I'm guessing you have conditions at both ends of the rod, and at both ends of the radius (r=0 and r=rod_radius), so no.

If collect(equ,tanh) doesn't do what you want, then you will need to explain exactly what you want to do. Your tanh's are inside JacobiSN etc, so probably you can only simplify the arguments of JacobiSN etc in some way. But what way? You also have sech, sinh, cosh; do you want to convert them all to tanh? Does convert(equ,tanh) help?

As far as I know, Maple can't directly solve 2-d recurrences, but it does have various tools to help you find a solution. If p was a known integer, then this probably could be done by first making generating functions in one of the variables, solving the 1-d recurrence in the generating functions, then finding the coefficients by Laplace inversion. I'm less certain if p is unspecified. So what are the restrictions on p?

But in any case, what is P(0,0)? You have P(0,m) = 1 but P(n,0)=0. Mariusz's solution gives P(0,0)=0 - is this correct?

@Rouben Rostamian  Yes, this also works in Windows (relative to currentdir(), which is usually the same as worksheetdir)

@C_R On Windows 11 and Maple 2026.1, print with Adobe or with "Microsoft print to pdf" give visually identical results - no strange characters. There are 3 t's on the second page before "new page". If there was supposed to be a page break in the .mw file I do not see it.

@janhardo In a document, ctrl-J inserts an execution group, but if you then type commands they are 2D if 2D is your global setting, or 1D if that is your global setting (i.e. the setting under File - options - display - input display.). ctrl-m/ctrl-r change to 1D/2D mode whether you are in an execution group or a document block.

Maple's help page ?DEtools,singularities uses "ordinary point" rather than "regular point" in the nonsingular case as does the DLMF, whiich I tend to use as an authoritative source. So I agree that ordinary point is probably better than the regular point used in the YouTube video, especially to avoid confusion with a regular singular point.

Those seem like sensible rules so that a quick reading in which some spaces aren't noticed doesn't lead to a misinterpretation.

a 2 c could be a variable name a2c but 2 a c can't be. 2 2 c could be 22 c.

See https://en.wikipedia.org/wiki/It%C3%B4%27s_lemma#Motivation

where the term with sigma is not involved in the expectation. But then I'm not sure why sigma would remain in the phase factor. 

@salim-barzani Notice that there is a u on the rhs side of the pde, so I expect to be able to cancel the exp parts, and this does work out.

You forgot to differentiate zeta(t) in the pde. And the paper had missing parentheses - if the exp bit is supposed to be the phase factor, then it must be exp(I*(X+Y)) not exp(I*X+Y).

But you still have to set sigma to zero, so I think this has to be part of the way to deal with stochastic pdes. Sigma still appears in the phase factor, so the solution has randomness still.

F-ode-Fpaper.mw

@salim-barzani I was doing this here. But there is still not enough information and the two papers are doing different things.

restart

pde := diff(Omega(x, t), t)-tau*(diff(B(t), t))

diff(Omega(x, t), t)-tau*(diff(B(t), t))

eq := Omega(x, t) = M(x, t)*exp(tau*B(t)-(1/2)*tau^2*t)

Omega(x, t) = M(x, t)*exp(tau*B(t)-(1/2)*tau^2*t)

The exp terms do not disappear unless the last term has tau=0 or diff(B(t),t)=0 or we make tau*B(t)=tau^2*t/2

ode := collect(eval(pde, eq), exp)

(diff(M(x, t), t)+M(x, t)*(tau*(diff(B(t), t))-(1/2)*tau^2))*exp(tau*B(t)-(1/2)*tau^2*t)-tau*(diff(B(t), t))

eval(ode, B(t) = (1/2)*tau*t)

diff(M(x, t), t)-(1/2)*tau^2

If we did this before the pde, we just get a constant term of -tau^2/2 in the pde - This is why you got cancellation of B(t) in pdetest.mw.

eval(pde, B(t) = (1/2)*tau*t)

diff(Omega(x, t), t)-(1/2)*tau^2

But if we do the following, which is what Claude implies (step 5 has -epsilon*M*V to get the cancellation (M should be diff(B(t),t))

pde2 := diff(Omega(x, t), t)-tau*Omega(x, t)*(diff(B(t), t))

diff(Omega(x, t), t)-tau*Omega(x, t)*(diff(B(t), t))

Then the exp terms can cancel

eval(pde2, eq); simplify(%/exp(tau*B(t)-(1/2)*tau^2*t))

(diff(M(x, t), t))*exp(tau*B(t)-(1/2)*tau^2*t)+M(x, t)*(tau*(diff(B(t), t))-(1/2)*tau^2)*exp(tau*B(t)-(1/2)*tau^2*t)-tau*M(x, t)*exp(tau*B(t)-(1/2)*tau^2*t)*(diff(B(t), t))

-(1/2)*M(x, t)*tau^2+diff(M(x, t), t)

But if this was correct, then we wouldn't need to set exp(tau*B(t))=exp(tau^2*t/2), which is one step in the first paper.
On the other hand for the second paper they do set up Eq 1.1 so this cancellation (for some terms) will occur. So the two papers are very different.

Conclusion is that there is insufficient information to know what to do. I suspect there is some hidden assumption in the way the Wiener process works. For example in the second paper, Eq 3.4 has a large epsilon (expectation?) - where did that come from. They do not even mention what it is.

NULL

NULL

Download idea.mw

@salim-barzani I agree about the tau=0. The real problem is earlier, where in both cases the exp() did not cancel. I'll work through the first paper again based on the claude.ai explanation. Some comments on the attached.

pde1.mw

pdetest.mw

@salim-barzani 

In the first paper, eq 4.1 is confusing as I said, but even more so since 4.15 is different yet also gets substituted into 4.2 and afterwards eq 4.20 still has the exp part from eq 4.1. So I can't figure this out.

@salim-barzani Please provide a link to this paper.

1 2 3 4 5 6 7 Last Page 1 of 105