I've been using the following syntax to set boundary condition which is a derivative, when passing it to pdsolve. Say we want to set u(r,theta,t) to have insulated boundary conditions at r=1. So the BC will be

For example, to set derivative of u w.r.t. "r" to zero when r=1

       eval(  diff(u(r,theta,t),r), r=1) = 0;  #(1)

or using this syntax

       D[1](u)(1, theta, t) = 0;  #(2)

But now I find, on one example below, that the above no longer works.  I have to use this syntax (which I did not know about) for it to work

        D[1]*u(1, theta, t) = 0;  #(3)

Has something changed? why when using (3) pdsolve now gives result, but when using (2) or (1) it returns unevaluated? are the three semantically equivalent? when to use which syntax?

I am using Physics updates 265, Latest Maple 2018.2 

Here is an example showing the (1,2)  syntax no longer works, but the (3) syntax works

#articolo example 6.9.2
restart;

#using (1) syntax
pde := diff(u(r, theta, t), t) = (diff(u(r, theta, t), r)+r*(diff(u(r, theta, t), r, r))+(diff(u(r, theta, t), theta, theta))/r)/(25*r);
bc_on_r := eval(diff(u(r,theta,t),r), r=1) = 0;
bc_on_theta:= u(r,0,t)=0, u(r,Pi,t)=0;
ic := u(r,theta,0)=(r-1/3*r^3)*sin(theta);
pdsolve([pde, bc_on_r,bc_on_theta,ic], u(r, theta, t), HINT = boundedseries(r = [0]))

does not solve it.

restart;

#using (2) syntax
pde := diff(u(r, theta, t), t) = (diff(u(r, theta, t), r)+r*(diff(u(r, theta, t), r, r))+(diff(u(r, theta, t), theta, theta))/r)/(25*r);
bc_on_r := D[1](u)(1, theta, t) = 0; 
bc_on_theta:= u(r,0,t)=0, u(r,Pi,t)=0;
ic := u(r,theta,0)=(r-1/3*r^3)*sin(theta);
pdsolve([pde, bc_on_r,bc_on_theta,ic], u(r, theta, t), HINT = boundedseries(r = [0]))

does not solve it.

restart;

#using(3) syntax
pde := diff(u(r, theta, t), t) = (diff(u(r, theta, t), r)+r*(diff(u(r, theta, t), r, r))+(diff(u(r, theta, t), theta, theta))/r)/(25*r);
bc_on_r := D[1]*u(1,theta,t)=0;
bc_on_theta:= u(r,0,t)=0, u(r,Pi,t)=0;
ic := u(r,theta,0)=(r-1/3*r^3)*sin(theta);
pdsolve([pde, bc_on_r,bc_on_theta,ic], u(r, theta, t), HINT = boundedseries(r = [0]))

I've used syntax (1) before many times and it works. Here is an example where all three syntax work

pde:=diff(u(x,t),t)=k*diff(u(x,t),x$2);
bc:=eval(diff(u(x,t),x),x=0)=0,u(L,t)=0;
ic:=u(x,0)=f(x);
sol:=pdsolve([pde,bc,ic],u(x,t));

pdsolve gives

 Using syntax (2)

pde:=diff(u(x,t),t)=k*diff(u(x,t),x$2);
bc:=D[1](u)(0,t)=0,u(L,t)=0;
ic:=u(x,0)=f(x);
sol:=pdsolve([pde,bc,ic],u(x,t));

gives same answer as (1) and using syntax (3)

pde:=diff(u(x,t),t)=k*diff(u(x,t),x$2);
bc:=D[1]*u(0,t)=0,u(L,t);
ic:=u(x,0)=f(x);
sol:=pdsolve([pde,bc,ic],u(x,t));

The answer also looks like different and simpler, but I assume they are equivalent for now without looking too much into it.

Which syntax should one use as now I am really confused.

It looks like (3) is the one that should be used? Why the others did not work on first example? i.e. pdsolve did not give an answer at all?  And if (1,2,3) syntax are supposed to be equivalent, whysecond example gives slightly different looking answer when using one syntax vs. the other?

And finally to make things more confusing, here is an example where syntax (3) does not work, but syntax 1 and 2 work:

#articolo example 8.4.3
restart;
pde := diff(u(x, t), t) = (1/20)*(diff(u(x, t), x, x))+t;
bc := u(0, t) = 5, (u(1, t)+ D[1](u)(1, t)) = 10;
ic:= u(x, 0) = -40*x^2*(1/3)+45*x*(1/2)+5;
pdsolve([pde, bc,ic], u(x, t))

gives answer.

restart;
pde := diff(u(x, t), t) = (1/20)*(diff(u(x, t), x, x))+t;
bc := u(0, t) = 5, (u(1, t)+eval(diff(u(x,t),x),x=1))  = 10;
ic:= u(x, 0) = -40*x^2*(1/3)+45*x*(1/2)+5;
pdsolve([pde, bc,ic], u(x, t))

gives same answer. But

restart;
pde := diff(u(x, t), t) = (1/20)*(diff(u(x, t), x, x))+t;
bc := u(0, t) = 5, (u(1, t)+ D[1]*u(1, t)) = 10;
ic:= u(x, 0) = -40*x^2*(1/3)+45*x*(1/2)+5;
pdsolve([pde, bc,ic], u(x, t))

does not work.

Clearly there is something I do not understand between these 3 syntaxes and when to use which.

Using 265 version

Physics:-Version()
 "C:\Maple_updates\Physics+Updates.maple", 2018, December 22, 

    10:41 hours, version in the MapleCloud: 265, version 

    installed in this computer: not installed

downloaded today.

When I try to install (as an example) physics update using Maple 2018.2.1 on windows 10, it keeps hanging in the middle as shown above.

I closed Maple, starting again and tried again, same thing happens. I waited for more than 10 minutes and nothing happens.  I click on the install button in the clouds windows from Maple itself to install it as I always did.

Do others have problem installing this?  I'll try again in few hrs, it might be the Maple server is having some issues but thought to ask.

I try to solve ODE with conditions, but it give answer only without conditions:

SOT.mw

Thank you.

Hello,

Suppose we have a set of quadratic equations of the form:

U_11 * a * q + U_12 * b * q + U_13 * c * q + U_14 * d * q + U_15 * a * w + U_16 * b * w + U_17 * c * w + U_18 * d * w + .. = C_i

.

.

.

Where terms in uppercase means constant while lowercase letters correspond to unknowns.

Now, I want to make a change of variable, so instead of having non-linear term `a * q` we would have `s_ij`, meaning i-th equation and j-th unknown after the change was done.

I'm trying to do so, because in my case I will be left with a linear system with a small number of unknowns(relative to N) and N equations, and so I could then solve it easily.

Any help is appreciated, thank you!

I need Maple code of Bisection method.

Download maple_prime.mw

Download int-too-many-levels-of-recursion.mw

Hello again, I posted a thread here earlier and received some great responses. I've made some progress in Maple since then but once again I've ended up on a question where I am completely stuck. I am only meant to solve this one in maple as the answer is written below the question. The question is:

https://imgur.com/rhFNh2l

I realise this is alot to ask for but I'm studying from afar and I don't have as many options for help as other students at the moment. I have managed to solve a), and I know how to solve b) the traditional way (pen and paper), but I have no idea how to do it in maple. My solution for a) is attached below. c) and d) strictly rely on the answer from b) so I'd greatly appreciate if I could have some help with it so I atleast could attempt the others on my own. This is the final question I have so once I'm done with this I'll pretty much be a master Maple... or not.. :P

I have commented the maple document so it is easier to understand what I've done and what I want help with. Also, I very much apologise if something that I write don't make sense, English is not my native language.

Hello, I just started working in Maple and I am struggling with a question that I found in my book. The question goes as follows:

Consider the set of vectors T = {v₁, v₂, v_{3​​​​,} v₄} in ℝ⁴ with

v₁ = (2,5,-2,0), v₂ = (-2,-4,b,4), v₃ = (-1,2,-2,-5), v₄ = (b,2,5,3)

Find all values of b, such that T is a linearly independent set

I am supposed to first solve it by hand and then in maple. I can solve it by hand but I have no clue how to do it in maple.The only thing I've done so far is create a matrix that contains the vectors in the set S as column entries.

I know that the columns of A are linearly independent if  det A ≠ 0. This is where my near non-existent knowledge in Maple stops me...

The answer is for all b {2}

I appreciate any tips...

Hi maple users

I'm about to do an assignment on the Joukowski transformation, where I have modelled and airfoil resembling the NACA-64015 airfoil. The parameters of the circle in z-plane are as follows: r=1.1241, x=-0.1241, y=0.0, where r is the radius of the circle, and x and y is the parameters used to offset the center of the circle. The problem is that I can't seem to plot the airfoil with contuours showing the fluid flow. I would therefore like to know if any of you guys has solved this before? Any help is much appreciated.

Yesterday, I accidentally discovered a nasty bug in a fairly simple example:

restart;
Expr:=a*sin(x)+b*cos(x);
maximize(Expr, x=0..2*Pi);
minimize(Expr, x=0..2*Pi);
                                    

I am sure the correct answers are  sqrt(a^2+b^2)  and  -sqrt(a^2+b^2)  for any real values  a  and  b .  It is easy to prove in many ways. The simplest method does not require any calculations and can be done in the mind. We will consider  Expr  as the scalar product (or the dot product) of two vectors  <a, b>  and  <sin(x), cos(x)>, one of which is a unit vector. Then it is obvious that the maximum of this scalar product is reached if the vectors are codirectional and equals to the length of the first vector, that is, sqrt(a^2+b^2).

Bugs in these commands were noted by users and earlier (see search by keywords bug, maximize, minimize) but unfortunately are still not fixed. 

How I can extract data about figure as matrix format.

I want to save data from figure as matrix for each curve in figure as i.e.,  x and y .

for example as this format.

fig10.mw

sol := Vector[column](4, [` 1001 x 4 `*Matrix, `Data Type: `*anything, `Storage: `*rectangular, `Order: `*Fortran_order])

Is it possible to export data from Maple to Excel where data are generated by Maple using Optimization:-LSSolve or DirectSearch:-SolveEquations function?

I have attached a file, where n is the primary variable. For n=1, 2, 3, .......... 20, I have got solution individually. But it is more time consuming event. Is it possible to use array or a vector for the variable n and solve equations at a time for all value of n and stored a excel file?
Any advice for me? @Preben Alsholm , @tomleslie or anybody
system_equations_(10.12.2008).mw