I have a region x^2 + y^2 <= 1 and y>=0. It's temperature function is f(x,y) x^2 - 2y^2 + x + y. How do I find the max and min temperatures on the lower boundary y=0?

I took the derivatives with respect to x and with respect to y such that:

fx:=diff(f(x,y),x);

fy:=diff(f(x,y),x);

Then I used fsolve({fx=0, fy=0},{x,y}) which game me (-0.5, 0.25)

Is there really only one critical point on that lower bound...

Hi, I'm trying to have my half circle divided into equal gaps of 33 , but am not sure what to use. I tried using Array, Sequence, Matrix... But don't think i'm doing it right. Here is what I have :

gap :=pi/(16*2)

I want a sequence of [0,0.0982,0.1963,0.2945,0.3927etc] in a vector form. May I know what is the proper way to do it please?

Many thanks.

Please write for me a code this problem.

Write the equation of the plane passing through the two points A(1,-2,4), B(3,5,-1) and make a least angle with the plane x + y + z + 1 = 0.

Thank you very much.

Why in next code i must apply operator twice?

zzz.mw

-------

inert form eval procedure

Let A(1, 2, 3) be a point; d1: (x-2)/2 = (y+2)/(-1) = (z - 3) / 1

d2: (x--)/(-1) = (y - 1)/2 = (z  + 1) / 1.

Write the equation of the line passing through the point A, perpendicular to d1 and cut d2.

This is my code

restart;with(LinearAlgebra): with(geom3d):

A:=<1,2,3>:

B:=[2*t+2,-t-2,t+3]:

M:=<-m+1, 2*m+1,m-1>:

line(d1,B,t):

a:=convert(ParallelVector(d1),Vector):

v:=M-A:

Is there programs or algorithm to obtain resultant of multivariant polynomials for more than 2 polynomials

Dear Sirs

How to change font there?

==========

Debugger window

restart;
a:='a';
                               a
assume(a>0);
H:=a;
                               a
subs(a=b, H);
       ...

Let M(1; -1; 0) be a point, Delta: (x-2)/2 =(y+1)/(-1) = (z - 1)/1 be a line, (P): x + y + z - 2 = 0 be a plane. Find the coordinates of the point A lies on (P), knowing that the line AM perpendicular to the line Delta and distance from the point A to the line Delta equal to sqrt(33/2).

This is my code

 with(geom3d):

point(M,1,-1,3):

line(Delta,[2*t+2,-t-1,t+1],t):

plane(P,x+y+z+3=0,[x,y,z]):

a:=ParallelVector(Delta):

Is there a bug in the CodeTools:-EncodeName / DecodeName pair?

Consider the following minimal example:

with(CodeTools):
M := module() option package; export f; f := proc() end proc; end module;
EncodeName(M:-f);
DecodeName(%);

The EncodeName gives answer:

_Inert_ASSIGNEDLOCALNAME("f", "PROC", 0, _Inert_ATTRIBUTE(_Inert_NAME("protected", _Inert_ATTRIBUTE(_Inert_NAME("protected")))))

But applying...

Sometimes it is appealing to have a package export a procedure which does not show up when calling with().

For example, the procedure might be used elsewhere, but be otherwise so very technically obscure that nobody else would be interested. (The counter-argument is that what is good for the goose is good for the gander! If something else in Maple can make good use of it and need it, then you might too.)

Now, Maple's modules don't have the concept of "friends",...

how to find concavity of f(t)=(800-t)t/(100+t)

Hi:

Consider this trivial (very contrived) example:

restart:

a:=lambda+3;
b:=Vector([lambda-2, lambda+23]);

subs(lambda=1, a); #works as expected

subs(lambda=1, b); #works as expected

 
assume(lambda, real);  #seems like this causes a problem
subs(lambda=1, a);  #works as expected

subs(lambda=1, b);  #This substitution does NOT work - lambda stays as 'lambda'

I am using maple 15 on...

I converted an ode using the built-in "convert" tool to check some calculations I had done by hand. To my surprise, there was an inconsistency. I converted the ode using PDEtools[dchange], reproducing the steps I had followed manually, and they checked out. So my question is: is there a sign error in convert? (and therefore a bug) or are both conversions correct, and if so are there any lessons to be learned? (is it related to the equation's symmetries?)

Thanks for your comments.