Hi!

I'm currently using Maple 17 and I'm trying to programatically export .eps 3D plots using the following code:

plotsetup(eps, plotoutput = `square_of_x`, plotoptions = `color=rgb , portrait, noborder, height=4in,width=4in, shrinkby=0.1`)

The problem that I'm facing is:

1- The z-axis label is outside of my figure (see uploaded figure). This happens even if I try to save the figure in .png and it seems that the 'shrinkby' option its not working properly;

Could you help me solving this?

Thanks!

Hi

I have a set P which is filled with 2D elements e.g.

P:={[1,1],[1,2],[1,3]...}

Later I'm wanting to add a new element e.g. {[1,10]}
At the moment I have

P union {[10, 10]}

If I were to print this I get what I want, but I want to update P.
When I try and do

 P:= P union {[10,10]}

I get an error "Error, illegal use of an object as a name"

What am I doing wrong and how can I solve this? Much appreciated!

I have upload the attachment here. in this problem I want to find the comparison. Please check it and please tell me the correction

Ask.mw

AB44.mwAdded these tasks by highlighting object and then choosing
Edit > Create Task

Is there some way to delete a task from my Tasks palette

At some point I created "My second palette" is there a way
to delete individual items from this palette and how are tasks
added to a customized palette.

Any help appreciated, Thanks! Les

AB44.mw

Hello!

Hope everything fine with you. Please check the following program

restart; a := 0; b := 1; Rd := .5; Ra := .2; Nc := .6; e := .2; Nr := .3; H := .9; Digits := 20;
with(LinearAlgebra);
L[0] := 1; L[1] := a+1-x;
M := 4; K := 1;
for n while n <= M do
L[n+1] := expand(-((x-2*n-a-1)*L[n]+(n+a)*L[n-1])/(n+1))
end do; printlevel := 2;
for m from 0 while m < M do
for n to 2^(K-1) do
psi[n, m] := simplify(sqrt(m+1/2.)*2^((1/2)*K)*(eval(L[m], x = 2^K*x-2*n+1)))
end do
end do;
f1 := simplify(sum(sum(c[p, q]*psi[p, q], q = 0 .. M-1), p = 1 .. 2^(K-1)));
u := collect(f1, [x], recursive);

u := -3.5276684147527874541*c[1, 3]*x^3+(4.4721359549995793928*c[1, 2]+21.166010488516724724*c[1, 3])*x^2+(-3.4641016151377545871*c[1, 1]-13.416407864998738178*c[1, 2]-34.394767043839677677*c[1, 3])*x+c[1, 0]+3.4641016151377545871*c[1, 1]+7.8262379212492639373*c[1, 2]+14.992590762699346680*c[1, 3]

f := 1+sum(A[i]*(x^(i+1)-(i+1)*x), i = 1 .. M-2)

f:=1+A[1]*(x^2-2*x)+A[2]*(x^3-3*x)

I want to equate the like powers of "x" in f ( mention as red) and u ( mention as red) for any values of M make a system of equations as

Eq[1]:=-3.5276684147527874541*c[1, 3] = A[2];

Eq[2]:=4.4721359549995793928*c[1, 2]+21.166010488516724724*c[1, 3] = A[1];

Eq[3]:=-3.4641016151377545871*c[1, 1]-13.416407864998738178*c[1, 2]-34.394767043839677677*c[1, 3] = -2*A[1]-3*A[2];

Eq[4]:=c[1, 0]+3.4641016151377545871*c[1, 1]+7.8262379212492639373*c[1, 2]+14.992590762699346680*c[1, 3] = 1;

and solve them for any values of A[1] and A[2];

Also attached the program Equate.mw

With my best regards and sincerely.

Muhammad Usman

School of Mathematical Sciences 
Peking University, Beijing, China

For the PDE, I can't seem to plot T' as the y-axis, it gives me the error, "Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct". Anyone knows what seem to be the problem? I am open to all ideas and would appreciate any help:)

fyp3.mw

Hello, 

I have just started using Maple, and it seems very powerful. I am trying to solve trigonometric equations and get all the solutions in a range, but when I use fsolve I only get one solution. 

Is this by design of the function or is there another way to do this? 

Tom

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)

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