Thank you in advance for your time.

I can search for any phrase or symbol in my Maple Worksheet but "sign symbol".

I have a very big result that I am not sure how many minus it has. When I search for - sign, no match found. I can see some minus in the worksheet but Maple can not find them.

Please help me.

I really appreciate.

if there is a funciton like f(i,j)=i+2*j
and I have already print out a 3D picture use the commend
"plot3d(f,i=1..10,j=1..10);"
the 3D picture just like the link : http://imgur.com/eLYhYEs
how can I convert the datas in 3D picture to text file?
the type may be like as follow I think.
-------------------------------------
(i,j)=f(i,j)
(1,1)=3
(1,2)=5
(1,3)=7...
-------------------------------------
I had try
"writedata("C:\\desktop\\a.txt",3Ddatas)"
but it seems can only output one variable parameter in a function, what else can I do?
thanks a lot.

Dear all;

I need you to understand this problem...

when i plot the function, using the graph i see that the function is above the x-axis but when I compute some values of this function I get a negative values....like
evalf(y(99.6));  is a negative value, but in the graph it is possible... I don't  undertand the problem...

restart:
with(plots):

# funciton

y:=x->-4.1123583570*10^281*exp(-(2/3)*x^(3/2))/(x^(1/4)*sqrt(Pi))+1.6554662320*10^(-289)*exp((2/3)*x^(3/2))/(x^(1/4)*sqrt(Pi))+(16/153)*x^(7/6)*sqrt(Pi)*exp((2/3)*x^(3/2))+Pi*((1/2)*exp(-(2/3)*x^(3/2))*(-1+exp((2/3)*x^(2/3)))/(x^(1/4)*Pi)-(16/153)*x^(7/6)*exp((2/3)*x^(3/2))/sqrt(Pi)):

#I plot this function in the interval (a,b)

a:=99;b:=100; # interval (a,b)
 forget(evalf): Digits:=20:
P1:=plottools:-transform((x,y)->[x+a,y])(plot(expand(y(x+a)),x=0..1,color=blue)):
forget(evalf): Digits:=4000:
P2:=plot(ysol, a..b, style=point, adaptive=false, numpoints=25, symbol=solidcircle, symbolsize=20, color=blue):
Digits:=20:
plots:-display(P1,P2);
evalf(y(99.6)); 

 Thank you in advantage for your remarks

how_do_i_get_floating_point_answer.mw

Download how_do_i_get_floating_point_answer.mw

Ramakrishnan V

rukmini_ramki@hotmail.com

Dear all;

Thank you for helping me, to plot two function on the same graph.

restart;
with(plots):

ode := diff(y(x), x, x) = x*y(x)+x;
a:=2; b:=3;  # we work in the interval (a,b)
ics := y(a) = 0, y(b) = 1;
sol:=dsolve({ode,ics}, numeric):  # Solution of my ode

# First figure: solution of my ode in the interval (2,3)

odeplot( sol,[x, y(x)], x=a..b);  # I plot the solution in the interval (a,b)

# Here I define a second function

z:=x->-4.9354831550*exp(-(2/3)*x^(3/2))/(x^(1/4)*sqrt(Pi))+0.752447265e-1*exp((2/3)*x^(3/2))/(x^(1/4)*sqrt(Pi))+(16/153)*x^(7/6)*sqrt(Pi)*exp((2/3)*x^(3/2))+Pi*((1/2)*exp(-(2/3)*x^(3/2))*(-1+exp((2/3)*x^(2/3)))/(x^(1/4)*Pi)-(16/153)*x^(7/6)*exp((2/3)*x^(3/2))/sqrt(Pi));

# Here a code to plot this function z

forget(evalf):
Digits:=20:
P1:=plottools:-transform((x,z)->[x+a,z])(plot(expand(z(x+a)),x=0..1)):
forget(evalf): Digits:=4000:
P2:=plot(z, a..b, style=point,adaptive=false, numpoints=25):
plots:-display(P1,P2);  # Second figure

My question: I want one graph that shows the solution of my ode and the funciton z. ( i.e plot the two figures in only one graph).

Thanks for helping me.

A strange verical line in output of dsolve.Should be a comma in Jacobi SN function and comma's in invert Jacobi SN function.

Possible_Bug.mw

Can anyone help me confirm this?

Hello,

How to change the length underline in latex ?

For example in this environment :

$\underline{\mathcal{G}}$

Thank you,

Gérard.

I tried doing LinearSolve with my own matrices, thought maybe there was a mistake there so I copied the example from Maple's website and it won't execute. Instead of printing an output matrix, it just prints LinearSolve(input matrix 1, input matrix 2). I feel this is a very simple syntax mistake, but I can't figure it out!

Hello,

While doing some tests, we encountered an issue with LinearAlgebra[ConditionNumber], possibly caused by another routine in the LinearAlgebra package (MatrixInverse?). I include an example worksheet bug-maple.mw for easy testing.

To reproduce, you should start by creating this matrix:

restart;
with( LinearAlgebra ):
Digits := 50:
n := 4;
M := [ seq( cos( Pi*j/n ), j=0..n ) ];
V := VandermondeMatrix( M );

Now compute the condition number (infinity-norm) of V:

c1 := ConditionNumber(V);
c1 := simplify(%);
evalf(c1); # this gives exactly 40

We can still compute this exactly and get 40.

This is the same result as with linalg[cond]:

c2 := linalg[cond](V);
c2 := simplify(%);
evalf(c2);

Now increase the value of n to 5 and redo the test.
The exact result should be 120, but LinearAlgebra[ConditionNumber] gives the wrong result!
On the other hand, linalg[cond] gives the correct, exact result (after simplification).

n := 5; # this gives a wrong result :-(

We could use Norm(V) * Norm(Inv(v)) since this should be equal to the condition number.
But then it seems there's something bad inside the MatrixInverse routine.
Again, linalg[inverse] gives the correct, exact result.
So compare the two results:

n := 4;
simplify(MatrixInverse(V));
linalg[inverse](V);

Finally, to make matters even worse, ConditionNumber gives completely wrong results when using a float-matrix, while linalg[cond] does not:

n := 3;
V := VandermondeMatrix( evalf(M) );
c1 := ConditionNumber( V );
c2 := linalg[cond]( V );

For n=3, we get c1=10.6 and c2=16.
For n=4, we get c2=5 and c2=40.
For n=5, we get c1=68.07 and c2=119.999.

Is this what it looks like? With it being a really, really ugly bug?

-- Many thanks for your reply,

Franky

This is my code

am:=proc(dx,Nminus,Nplus)
global k_r,n,tau,a2,m;
k_r:=1/2;
a2:=a/2;
k=0;
u:=[]:
for n from Nminus to Nplus do
k:=k+1:
u(k):=max(exp((0.5*(k_r+1)^2)*0)*(exp(0.5*n*dx*(k_r+1))-exp(0.5*n*dx*(k_r-1))),0);
u:=[op(u),u(k)];
end do;
end proc;

I keep getting the error.

Error, (in am) invalid left hand side in assignment.

I would like to fill the array u with [u(1),u(2),...,u(Nplus-Nminus)]

Please help!

Dear friend,

Recently I noticed, that numerical integration returns different values for the same function.

For example the code

restart;
evalf(int((exp(x)*(4420*cos(4)*sin(4)-544*cos(4)^2+148147*exp(-1)-4225*cos(4)-215203)/(71825*exp(1)-71825*exp(-1))-exp(-x)*(4420*cos(4)*sin(4)-544*cos(4)^2+148147*exp(1)-4225*cos(4)-215203)/(71825*exp(1)-71825*exp(-1))+(32/4225)*cos(4*x)^2+(1/71825)*(4225+(2210*x-6630)*sin(4*x))*cos(4*x)+x^2+8434/4225)^2, x = 0 .. 1));

each time returns values

0.0005951015934
0.0005950850548
0.0005950974588
0.0005950960805
0.0005951297843 etc.

Maybe, evalf uses a stochastic algorithm for integration?

So I am trying to apply the Optimization function to a rather complicated problem I am currently working on, and having some trouble getting maple to cooporate. This is the part of the code that is giving me the error, with error included. All of the variables imputtet when calling OptimizeSpring are constants. If there is any other details I should add, please say so and I will add them promptly :) Thanks alot!

I suspect that it might be related to this suggestion, but I am not sure how to apply it http://www.mapleprimes.com/posts/95113-Optimization-Workaround-For-Operator-Form

Hey guys! Can anyone help me with solving one of these Differential Equations ?
Thanks.

!

[Here's the original Question, restored as best I could.--Carl Love]

[The question was, of course, How to avoid the error?--Carl Love]

Download subgrouplattice_(1).mw

Hi

how to say 

i from 1 to 4 in maple

for example

if I have this eqauation

and I need to sub i in this equation to get 4 equation

Many thank