Im attempting some very basic calculations where numbers are stored in variables with correct units. However i wish to take a variable to the power 2 without altering it's unit. This is a really easy but strange calculation to make which is why i guess Maple is having trouble with it. 

e.g. I have a variable A with a number 4 stored inside with unit in Meters. I want to take A² = 16 meters and and not A² = 16 meters². 

Is there a function or way to select the base number without having to write an entire line about it? 

convert(A,'unit_free')^(2)*Unit('m')

can do the trick but i believe there has to be an easier, less complicated, way to achieve this? 

i need to find the graph of exp(v(r)) from 0 to 2e6.

i obtained the graph of v(r) by taking log but how to obtained in exp form

graph_of_exp_form_of_u(r).mw 

restart; with(LinearAlgebra); with(Student[LinearAlgebra]); 
Z := Matrix([[3, 1, 5, 4, 5], [0, 0, 0, 0, 0], [28, 6, 4, 5, 9], [98, 5, 82, 2, 4], [24, 55, 23, 22, 90]]); 
B := Matrix([[0], [0], [0], [0], [0]]); 
A := <Z|B>:
if Row(A, 2) = ZeroVector[row](6) then A; end if;

I used the code above but it didn't work...any help appreciated

Hello,

I would like to know if it is possible to create a dll function from a maple function (a procedure) which would be directly usable in a excel sheet.

In other words, similary to the VBA function that it is possible to build in Excel, i would like to create a dll function that i could use in Excel.

If it is possible, can you help me to use a good process?

Thank you for your help and feedback.

how to find the integration of z(x) form 0 to x with the given condition...
 

Download integration.mw

Hi,

I'm just doing a small project/HW for school where I am trying to create a histogram of the appearance of the digits of Pi. My problem is that I lose connection with the Maple Kernel whenever I go over 1000 digits. Why? Is this because my function scales badly? Any help would be appreciated.

Thanks in advance.

d:=1500;
Digits:=d;
                              1500
                              1500
valpi:=evalf(Pi):

dg:=proc(n,k)
floor(n/10^(ceil(log10(n)-k))) mod 10;
end:

with(Statistics):
L:=[seq(dg(valpi,i),i=1..d)]:
Histogram(L)

I apologise if this information is available on the support pages or on these forums already.

I have found the general solution to a differential equation and now need to plot multiple solution curves that correspond to a different value of a constant.

For example if I had y(x) = x + k, how would I go about plotting this for all the different values of k.

Many Thanks. 

There are many conventions for the Euler angles or other angles used to define o rotation of a 3d plot.
In Maple these angles are in the plot option orientation, but I think that the help page is not correct about them.
The same info appears in a worksheet (see ?rotateplot), so I am even more intrigued. [Note that many authors also switch phi and theta in spherical coordinates].

The help file says:

orientation=[theta, phi, psi]
This orientation specified by these angles is obtained by rotating the plot
1. psi about the x-axis,
2. then phi about the (transformed) z-axis, and
3. then theta about the (transformed) y-axis.
  These angles, given in degrees, are the Euler angles for the transformation matrix, using the axes specified. The angle psi is optional and is assumed to be 0 if not given. If the orientation option is not specified, the default value used is [55, 75, 0].

After some tests it seems that y and z should be switched, i.e. keeping the names (and order) for the angles ==>
1. psi about the x-axis,
2. then phi about the (transformed) y-axis, and
3. then theta about the (transformed) z-axis.
 
Am I right?

Download y(x).mw

I want to have a display() function that plots a certain expression that I already have defined, but then I also want to include a variable number of pointplots into the same graph. I basically want the user to define a set of x and y coordinates into two different arrays in the start of the Maple Worksheet, then later down the worksheet, I want to plot these x,y coordinates as a pointplot, for each point.

But for example, if there are 5 points, I would define the pointplots as p1,p2,p3,p4,p5, which I already have a loop to do automatically depending on however many points that has been defined in the previous arrays. But when I want to plot these in a display() function, I have to write display(p1,p2,p3,p4,p5) individually, and I can't therefore seem to find a way to make it like display(p1,...,pn). I want the display to add more pointplots depending on how many points are defined in the arrays.

Can anyone help me? Sorry if it was difficult to understand, I can explain further if you didn't understand.

I have four matrix equations

P1, P2, P3 are known 4x4 matrix.

A1 A2 A3 A4 are known 1x4 matrix.

x1 x2 x3 x4 are 1x1 known matrix.

U is 4x4 unknown matrix.

These equations are 

(A1^T*U*P1*A1) +( (P2*A1)^T*U*P1*A1) + ( (P3*A1)^T*U*A1) + ( ( P3*A 2)^T*U*P1*A1) + x1 =0;

(A2^T*U*P1*A2) +( (P2*A2)^T*U*P1*A2) + ( (P3*A2)^T*U*A2) + ( ( P3*A2 )^T*U*P1*A2) + x2 =0;

(A3^T*U*P1*A3) +( (P2*A3)^T*U*P1*A3) + ( (P3*A3)^T*U*A3) + ( ( P3*A3 )^T*U*P1*A3) + x3 =0;

(A4^T*U*P1*A4) +( (P2*A4)^T*U*P1*A4) + ( (P3*A4)^T*U*A4) + ( ( P3*A4 )^T*U*P1*A4) + x4 =0;

How can i find 4x4 matrix U by using these above four equations??

Thank you

Is it possible to display an internet webpage in a maple?

Hello everybody,

In the linearised theory of gravity, I want to do some symbolic calculations.  First, I need to set that:

Then I want to see how the Christoffel symbols will change by putting the above in this:

Any hint someone?  I really appreciate the help for learning the Physics package.  Thank you in advance.

Mario

I am running MAPLE 2016.1

I am able to set the color option to a function in fieldplot, e.g.

plots[fieldplot]([y, -sin(x)-(1/10)*y], x = -10 .. 10, y = -10 .. 10, arrows = SLIM, color = sin(x)*sin(y))

works as expected.  

However, for plots[fieldplot3d] I cannot set color = f(x,y,z) . E.g.

plots[fieldplot3d]([2*x, 2*y, 2*z], x = -1 .. 1, y = -1 .. 1, z = -1 .. 1, grid = [5, 5, 5], arrows = SLIM, color = sin(x)*sin(y)*sin(z));

Error, (in plot/color) invalid color specification: sin(x)*sin(y)*sin(z)

Please can anyone help?

MRB  

First Question: How to define nx1 matrix Y:=(y1,y2,...,yn) ? (n is a Natural number while it is ungiven)

Second Question: How to derivative of matrix Y with respect to nx1 matrix X:=(x1,x2,...,xn) ?

My efforts for the first question:  (I know to define 6x1 matrix etc. , but I dont know to define nx1 matrix

restart: Matrix(1..6,1,symbol=y) 

My efforts for the second question:

restart; with(VectorCalculus);
Matrix([x^2, x*y, x*z]);
Jacobian([x^2, x*y, x*z], [x, y, z]);

Can you help me?