## stability of the scheme...

i want to find the stability of this equation, but there is seem to have some problems..can somebody help me..

y := A*(1/x+x*exp(-2*sqrt(-1)*b))+4*(exp(h)-1)^2*(2*exp(-sqrt(-1)*b)-3*(exp(h)-1)^2*x^(-exp(h)+1)*exp(sqrt(-1)*b*(-exp(h)+1-1))+3*x^(-exp(h)+1)*exp(sqrt(-1)*b*(-exp(h)+1-1)))/(3*(1-r))-exp(-2*sqrt(-1)*b)/x-x;

subs(A = (1+r)/(1-r), %);

subs(r = (1/3)*(exp(h)-1)^2, %);

subs(b = m*Pi*(exp(h)-1), %);

subs(m = 1, %);

subs(h = 0.5e-1, %);

> ans := solve(%, x);
Warning, solutions may have been lost
ans:=

> r1 := ans[1];
Error, invalid subscript selector
> r2 := ans[2];
Error, invalid subscript selector

## plot multiple for or do loop display...

Dear  i want to plot the graphs for different values of R[2,5,6.5]   here is the codes: Either ido loop or array i have to use or other procedure

The graph display should be at R=2, one graph, R=5 another and R=6.5

onemore graph

restart:
with(DEtools):
eta:=1+k*x+epsilon*sin(2*Pi*x):
A1(x):=-exp(-alpha*x)*J^2*R/(2*eta^3+6*xi*J^2*eta^2):

psi0:=A1(x)*y^3:
psi:=delta*psi0:
V:=-diff(psi,x):
delta:=0.1:
epsilon:=0.01:
alpha:=1:
xi:=0.001:
k:=0.1:
As:= [0, 2, 4]:
x:=0.2
plot([seq]([eval(V, J= A), y, y= 0..eval(eta, J= A)], A= As),title=" velocity at R=2.0 ",labels= ["v", "y"],color= [green, red, blue], linestyle= [solid, dash, dot],legend= [seq](J = A, A= As), axes=boxed);

## Error, (in dsolve/numeric/bvp/convertsys) unable t...

Hye, have a good day everyone..

can anyone know why I got this error? :(second_grade.mw

Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system

Here I have attached the file. thank you in advance :)

## How do I check the monotonicity of a function?...

This might be a ridiculous question, but what is the easiest way to check the monotonicity of a function? I didn´t find anything helpful in the Maple Help, so maybe one of you could give me a nudge here :)

## how to use mathematica code convert to maple code...

I find mapleform software   but I dont know how to use , and does anyone know another method to convert.  I know use

> with(MmaTranslator);
> MmaToMaple();

I can  automatically translated .nb files.  but this have

Error, (in readline) file or directory does not exist
Error, (in readline) file or directory does not exist
Error, (in readline) file or directory does not exist
Error, (in readline) file or directory does not exist

another Error, missing operator or `;`

also I try to use

> with(MmaTranslator);
> FromMmaNotebook(Mma_notebook_filename, options);

I still dont know how is works?  can you explains for me and show me some example  , here is my example files.plot.zip

Hello everyone!

I'm trying to understand how to work with PDEs in Maple by solving a basic PDE such as the wave equation, along with conditions at t = 0 in terms of arbitrary functions. I stumbled upon the very same example, included in the Maple Online help at this page under the part "Three Textbook Examples":

I tried to copy and paste the code and execute it, but to no avail: it seems like pdsolve doesn't give me any other result than a blank space, unlike in the example reported in the Online help. Can anyone tell me if (and where) I am doing wrong? Is it possible that my version of Maple (Maple 13) may not be compatible with the procedures required by this example?

A warm thank you to whoever will help me!

## Kernel connection lost...

Hello,

I'm trying to calculate and export a huge "Strategy-Matrix". The calculation works in a smaller version with less possible strategies. For this version I use a specialised simulation-PC with enough memory (that was my problem using my own PC). But now I've got this error massage ("Kernel connection has been lost. You should save this worksheet and restart Maple. Executing commands in Maple requires a connection to the Maple kernel. Firewalls have been known to cause problems with kernel connections in Maple. Please ensure that any firewall software is configured to allow Maple to create connections to the kernel. Consult the FAQ for more information.")

The firewall should not be the problem. Below the code. num_Strat is 63001.

m_Gi:= Matrix(num_Strat, num_Strat):

for t from 1 by 1 to num_Strat
do
Digits:= 5:
v_h2:= evalf((-(24/17)* v_pEi[t]^2)*v_h1 + (v_pEi[t]*(-(128/51)*v_pKi[t]+980/51))*v_h1+ v_pEi[t]*((40/51)*v_pKj+(15/17)*v_pEj)-((110/51)*v_pKi[t]^2)*v_h1 + v_pKi[t]*((55/51)*v_pKj+1220/51*v_h1+(10/17)*v_pEj)):
m_Gi(..,t):= v_h2:
end do:
ExportMatrix(TestMatrixNK, m_Gi):

## Print "Hello" in Green?...

How do we print "Hello" in green in a .mw file in Maple 13?

Thanks!

mapleatha

## A Simple Example of a maple.ini File and Its Locat...

I cannot make Maple 13 read the maple.ini file no matter where it is located.

I am using Windows 7, 64 bit.

I thought that ...maple 13\bin.win would do it.

The file contains a few commands of mine. Three of these commands

"read" three procedure files A.m, B.m C.m. respectively, and the other two

are

with(LinearAlgebra):
with(DEtools):

I would appreciate it.

Thank you!

mapleatha

## There is no Startup Code Menu Item in Maple 13...

There is no menu item called "Startup Code"

in the "Edit" menu item of my Maple 13.

Thanks!

## Contourplot and surface-plot/plot[surfdata] in sam...

I have a surfdata-plot which - by interpolation - goes through different circles with different inclinations with respect to each other. The surface is color-coded with zhue from violet to red. In order to see the valleys behind peaks I would like to have a filled 2D-plot of my zhue-surfdata colors with contours wherever the surface passes z-coordinate integers located somewhere below the surface (e.g. at z=-7). Unfortunately, I haven't found examples of contourplots for surfdata plots so far.

 >
 >
 (1)
 >
 (2)
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 >
 (3)
 >
 >
 >

## Display a contour and 3d plot in same plot...

When trying to plot "Example1" from this page www.rhitt.com/courses/227/su99/Maple/html/surfaces2.html I get an error message that Maple cannot make a plot structure from an object with the transformation proc( (x,y) -> [x,y,-.3] end proc) . How can I plot this example?

## Define a vector for arbitrary values ...

Helo Dears, I need help to generate a vector whose detail given below

`for `K:=1 and m:=3 i need [[[psi[1,0]],[psi[1,1]],[psi[1,2]]]] and for K:=1 and m:=4 i need [[[psi[1,0]],[psi[1,1]],[psi[1,2]],[psi[1,3]]]]

`similarly for `K:=2 and m:=3 i need [[[psi[1,0]],[psi[1,1]],[psi[1,2]],[psi[2,0]],[psi[2,1]],[psi[2,2]]]] and for K:=1 and m:=4 i need [[[psi[1,0]],[psi[1,1]],[psi[1,2]],[psi[1,3]],[psi[2,0]],[psi[2,1]],[psi[2,2]],[psi[2,3]]]]

Simply I need to generate a general vector for arbitrary values of k and M as discus in above.

New_Wavelet_with_Operational_Matrix.mw

## Command for a square matrix...

Dear friends! Hope you will be fine. I want to generate a general form of the matrix (of order 2M+1 by 2M+1) shown below in Maple

[[[0,0,&ctdot;,0,-1/(M Pi),0,&ctdot;,0,1/(M Pi)],[0,0,&ctdot;,0,-1/((M-1) Pi),0,&ctdot;,1/((M-1) Pi),0],[&vellip;,&vellip;,&dtdot;,&vellip;,&vellip;,&vellip;, &utdot;,&vellip;,&vellip;],[0,0,&ctdot;,0,-1/(Pi),1/(Pi),&ctdot;,0,0],[1/(Pi),1/(Pi),&ctdot;,1/(Pi),1,1/(Pi),&ctdot;,1/(Pi),1/(Pi)],[0,0,&ctdot;,1/(Pi),1/(Pi),0,&ctdot;,0,0],[&vellip;,&vellip;,&utdot;,&vellip;,&vellip;,&vellip;,&dtdot;,&vellip;,&vellip;],[0,1/((M-1) Pi),&ctdot;,0,1/((M-1) Pi),0,&ctdot;,0,0],[1/(M Pi),0,&ctdot;,0,1/(M Pi),0,&ctdot;,0,0]]]

for M=4 we get

Matrix(9, 9, {(1, 1) = 0, (1, 2) = 0, (1, 3) = 0, (1, 4) = 0, (1, 5) = -(1/4)/Pi, (1, 6) = 0, (1, 7) = 0, (1, 8) = 0, (1, 9) = (1/4)/Pi, (2, 1) = 0, (2, 2) = 0, (2, 3) = 0, (2, 4) = 0, (2, 5) = -(1/3)/Pi, (2, 6) = 0, (2, 7) = 0, (2, 8) = (1/3)/Pi, (2, 9) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = 0, (3, 4) = 0, (3, 5) = -(1/2)/Pi, (3, 6) = 0, (3, 7) = (1/2)/Pi, (3, 8) = 0, (3, 9) = 0, (4, 1) = 0, (4, 2) = 0, (4, 3) = 0, (4, 4) = 0, (4, 5) = -1/Pi, (4, 6) = 1/Pi, (4, 7) = 0, (4, 8) = 0, (4, 9) = 0, (5, 1) = 1/Pi, (5, 2) = 1/Pi, (5, 3) = 1/Pi, (5, 4) = 1/Pi, (5, 5) = 1, (5, 6) = 1/Pi, (5, 7) = 1/Pi, (5, 8) = 1/Pi, (5, 9) = 1/Pi, (6, 1) = 0, (6, 2) = 0, (6, 3) = 0, (6, 4) = 1/Pi, (6, 5) = 1/Pi, (6, 6) = 0, (6, 7) = 0, (6, 8) = 0, (6, 9) = 0, (7, 1) = 0, (7, 2) = 0, (7, 3) = (1/2)/Pi, (7, 4) = 0, (7, 5) = (1/2)/Pi, (7, 6) = 0, (7, 7) = 0, (7, 8) = 0, (7, 9) = 0, (8, 1) = 0, (8, 2) = (1/3)/Pi, (8, 3) = 0, (8, 4) = 0, (8, 5) = (1/3)/Pi, (8, 6) = 0, (8, 7) = 0, (8, 8) = 0, (8, 9) = 0, (9, 1) = (1/4)/Pi, (9, 2) = 0, (9, 3) = 0, (9, 4) = 0, (9, 5) = (1/4)/Pi, (9, 6) = 0, (9, 7) = 0, (9, 8) = 0, (9, 9) = 0})

I want the general for of this matrix for any value of M

9_by_9.mw