How is vertical text created in a plot?  It would be similiar to using the labeldirection option set to "vertical", but my text in question will be an annotation within the field of the plot.  A more general question would be how to place rotated text within a plot.    

Hi everyone.

I wrote simple program with several functions. I would like to obtain parallel computing of my functions. I can't understand how to use nodes.

My programm.

Simple_example.mw

When I input

cos(k*Pi) assuming k::odd

The result is -1

but

(-1)^k assuming k::odd

The result is (-1)^k

Why doesn't the second expression work?

AOA... How are you all. I need the answer of the following question.

input in Maple: expand(exp(a+b)+exp(c+b))

output:  exp(a)*exp(b)+exp(c)*exp(b)

and

input in Maple: expand(exp(2a+b)+exp(3c+b))

output:  (exp(a))^2*exp(b)+(exp(c))^3*exp(b)

but i need exp(2a)*exp(b)+exp(3c)*exp(b)

PhD (Scholar)
Department of Mathematics

Say for a complicated function f(x), there are 10 solutions. I want to check for any of the solutions, the matrix R is defined or not.

Is there a way to force execution like the following ones?

seq(1/i,i=0..3);

for i from 0 to 3 do
    1/i;
end do;

I think the problem is that, Maple stops immediately when it's undefined. But the following terms could be defined.

Or is there a way to put this check into "solve"? That's

> solve(f(x),x, "some way of making sure that the solution has to be defined for the matrix R" )

Is there a way to achieve this?

Thanks

Hi guys.

I have variables that look like C[i,j,k,l], where i,j,k,l can range between 1..3 for example. Is there a simple way to create a pattern-based assumption for something like C[i,j,k,l] = C[j,i,k,l] for all i,j,k,l ?

Kind regards,

woggy

Maple returns the wrong answer to the following very simple double sum:

T:=(x,m)->sum( sum(binomial(m,k+p),k=1..m-p)*x^p,p=0..m-1);

Maple gives zero for all m>0.  However, the correct answer is a non-zero polynomial in x. 

For example, with Maple:

T(x, 3);
                               0

The correct answer is  x^2+4*x+7.

It would be helpful if others can check this double sum in their version of Maple and report results.

Does Maple have any tool or package that computes the Fourier & Fourier-Bessel series expansions of a given funtion "f(x)" over a specified interval "[a,b]"?

what are the examples of calculating devices and packages?

I set a physical model for my reseach,

equ1 := x^4-5*x

equ2 := 1/x+x^2+3

x is a function of time t and it meets diff(x(t), t) = equi-equ2.

I want to plot the curve of x(t) varying with time

When I use the following command

DEplot({x(0) = 0, diff(x(t), t) = equi-equ2},x(t), t = 0 .. 20)   it shows: Error, (in DEtools/DEplot) called with too few arguments

Who can tell me what is wrong with my calculation? Thanks

Hi I am working on the following problem (See below the line) from the text A Introduction to the Mathematics of Biology (Ch 2 pgs.21-23). I am working in Maple 18 and the code for this problem is from Maple 11 I think. If you look at the last line of code, I ran into a snag. Can anyone help????? We are in a study Gourp and alll are stumped.

_______________________________________________________________________________________________________

Problem # 2 - Find a fit for the cumulative US AIDS data as a polynomial function. Alsoo find an exponential fit for the data. (assume the Data I have inputed is correct if you don't have the book).

AIDS := [97, 206, 406, 700, 1289, 1654, 2576, 3392, 4922, 6343, 8359, 9968, 12990, 14397, 16604, 17124, 19585, 19707, 21392, 20846, 23690, 24610, 26228, 22768];
print(`output redirected...`); # input placeholder
[97, 206, 406, 700, 1289, 1654, 2576, 3392, 4922, 6343, 8359,

9968, 12990, 14397, 16604, 17124, 19585, 19707, 21392, 20846,

23690, 24610, 26228, 22768]

CAC := [seq(sum(AIDS[j]/(1000.0), j = 1 .. i), i = 1 .. 24)];
print(`output redirected...`); # input placeholder
[0.09700000000, 0.3030000000, 0.7090000000, 1.409000000,

2.698000000, 4.352000000, 6.928000000, 10.32000000,

15.24200000, 21.58500000, 29.94400000, 39.91200000,

52.90200000, 67.29900000, 83.90300000, 101.0270000,

120.6120000, 140.3190000, 161.7110000, 182.5570000,

206.2470000, 230.8570000, 257.0850000, 279.8530000]

Time := [seq(1981+(i-1)*(1/2), i = 1 .. 24)];
ln(CAC)=k*lnt+A
LnCAC := map(ln, CAC);
print(`output redirected...`); # input placeholder
[-2.333044300, -1.194022473, -0.3438997525, 0.3428802329,

0.9925107578, 1.470635510, 1.935571171, 2.334083760,

2.724054775, 3.071998629, 3.399328971, 3.686677031,

3.968441145, 4.209145378, 4.429661370, 4.615387808,

4.792578782, 4.943918402, 5.085810791, 5.207062453,

5.329074480, 5.441798471, 5.549406770, 5.634264465]
LnTime := map(ln, [seq((i+1)/(2*(1/10)), i = 1 .. 24)]);
print(`output redirected...`); # input placeholder
[ln(10), ln(15), ln(20), 2 ln(5), ln(30), ln(35), ln(40), ln(45),

ln(50), ln(55), ln(60), ln(65), ln(70), ln(75), ln(80), ln(85),

ln(90), ln(95), 2 ln(10), ln(105), ln(110), ln(115), ln(120),

3 ln(5)]
with(stats);
fit[leastsquare[[x, y], y = k*x+lnA]]([LnTime, LnCAC]);
print(`output redirected...`); # input placeholder
y = 3.293411005 x - 10.12289000
k := op(1, op(1, rhs(%))); LnA := op(2, rhs(`%%`)); A := exp(LnA);
print(`output redirected...`); # input placeholder
Error, invalid input: rhs received exp(LnA), which is not valid for its 1st argument, expr
Error, invalid input: rhs received exp(LnA), which is not valid for its 1st argument, expr
exp(LnA)

Hi,

So, I have the following problem. In my research, I am trying to prove that something is impossible. I do that by encoding it as a set of equations, and asking Maple to solve them, and if it finds that the system has no solutions, I got what I wanted. However, I really wanted to see Maple output this fact explictly, more or less like Mathematica (where I either get an error message, or I get "False" if I use the Reduce command). The reason for this is that I will embed this procedure in a loop that varies some of my starting assumptions, and I want Maple to output that the system has no solutions at each iteration. This is useful since each iteration takes more than half an hour to run, and there are dozens of them, so I want the program to provide me some intermediate feedback, rather than waiting for everything to finish to get an answer.

Is the question clear?

Thanks a lot for the help!

Solve, using 4000 miles for the radius of the earth.                                                                                              

A space shuttle is in circular orbit 150 miles above the surface of the earth. Approximate                                                                             

1. the speed
2. the time required for one revolution.

I keep getting an error about the while. It is calling it an invalid if statement.  It may relate to semi colons, but I cannot find any errors myself. Please help.

InterpolIDW2:=proc(n,rwin,mink,p,R,IR)       

local ix,iy,i,j,i0,i1,j0,j1,k;      

local count,wgt,top,bot,hp;      

print(`n=`,n,`rwin=`,rwin,`mink=`,mink,`p=`,p);      

hp:=p/2;        

for iy from 1 to n do      

for ix from 1 to n do           

if (IR(ix,iy)=0) then                

i0:=max(1,ix-rwin);  i1:=min(n,ix+rwin);                

j0:=max(1,iy-rwin);  j1:=min(n,iy+rwin);                 

k:=add(add(IR(i,j),i=i0..i1),j=j0..j1)                                                                

while (k<mink )do                              

if(i0>1) then                                       

i0:= i0-1;                                       

k:= k+ add(IR(i0,j),j=j0..j1)                              

end if;                              

if(i0<1) then                                        

i0:=i0+1;                                        

k:= k+ add(IR(i0,j),j=j0..j1)                              

end if;                              

if(j0>1) then                                        

  j0:=j0-1;                                       

  k:= k+ add(IR(i,j0),i=i0..i1)                         

  end if;                           

 if(j0<1) then                                

 j0:= j0+1;                                         

k:= k+ add(IR(i,j0),i=i0..i1)                             

end if;                         

end do;                  

top:=0;  bot:=0;                

for j from j0 to j1 do                

for i from i0 to i1 do                      

if (IR(i,j)<>0) then                           

wgt:=1/evalf[14](((ix-i)^(2)+(iy-j)^(2))^(hp));                           

top:=top+wgt*R(i,j);                           

bot:=bot+wgt;                      

end if;                

end do;                

end do;                

R(ix,iy):=top/bot;           

end if;      

end do;      

end do;  

end proc:

I'm trying to solve a series of equations and then graph them. I'm trying to solve for the variables involved:

values := solve({eq1, eq2, eq3, eq4, eq5, eq6, eq7, eq8, eq9,

eq10, eq11, eq12, eq13, eq14, eq15, eq16}, {a, b, b, c, d, e, f,

g, h, i, j, k, l, m, n, o, p});

. . . but it gives me this:

Warning, solving for expressions other than names or functions is not recommended. 

values :=