Maple 18 Questions and Posts

These are Posts and Questions associated with the product, Maple 18

How to resolve of this ODE System Error, the proble is "Error, (in dsolve/numeric/process_input) input system must be an ODE system, got independent variables {t, [S, E, In, R]}

Download TKO.mw

FourTrans.mwI am wondering if I can calculate this Fourier Transfrom analytically using Maple:

Hi! I need to simplify a polynomial over an arbitrary field F on a radical variable ν, such that for some power j, ν is in F.

the polynomial takes the form:

k0+a1ν1/j+...+anνn/j

but, since νj is in F, this divides every element νi/j into a cotient group of order j, so this can be rewritten as a radical extension of F, in the form: 

k1+(a1+...)ν1/j+...+(ak+...)ν(j-1)/j

where ki, ai are in F. I feel like this is a very straightforward technique for handwritten algebra, but i can't see a command for this on Maple 18. Nor the simplify(..., radicals) Or the combine(...,radicals) seems to help here. Maybe there's an special command for this? Must i do it myself? Please, any help is aprecciated.

Hi maple users 

I am working with the PDE solver.

i am receiving a following error

"unable to compute solution for t>HFloat(0.0):" 

Kindly do the needful how to rectify this error.

dumm.mw

How can we determine the highest (or peak) value for every hump and at what value of (x, t) in the plot?

peak.mw

How to find the dispersion relation (3) in Maple? I try it by equating the coefficients of cos(kx + wt) and its derivative to zero, but couldn't find it.

restart

with(LinearAlgebra):

alias(E = E(x, t), P = P(x, t), N = N(x, t), Q = Q(x, t));

E, P, N, Q

(1)

d := Matrix([[c*(diff(E, t))+diff(E, x), alpha*P, 0, 0], [N[0]*E, diff(P, t), -sigma[2]*Q[0]*N/N[0], sigma[2]*omega[0]*Q], [0, sigma[1]*sigma[2]*Q[0]*P/N[0], diff(N, t), 0], [0, -omega[0]*P, 0, diff(Q, t)]]) = 0

d := (Matrix(4, 4, {(1, 1) = c*(diff(E(x, t), t))+diff(E(x, t), x), (1, 2) = alpha*P(x, t), (1, 3) = 0, (1, 4) = 0, (2, 1) = N(x, t)[0]*E(x, t), (2, 2) = diff(P(x, t), t), (2, 3) = -sigma[2]*Q(x, t)[0]*N(x, t)/N(x, t)[0], (2, 4) = sigma[2]*omega[0]*Q(x, t), (3, 1) = 0, (3, 2) = sigma[1]*sigma[2]*Q(x, t)[0]*P(x, t)/N(x, t)[0], (3, 3) = diff(N(x, t), t), (3, 4) = 0, (4, 1) = 0, (4, 2) = -omega[0]*P(x, t), (4, 3) = 0, (4, 4) = diff(Q(x, t), t)})) = 0

(2)

re := -c*omega^3-kappa*omega^2+(c*omega[0]*(N[0]^2*sigma[2]+Q[0]^2*sigma[1])/N[0]^2-alpha*N[0])*omega+kappa*omega[0]*(N[0]^2*sigma[2]+Q[0]^2*sigma[1])/N[0]^2 = 0;

-c*omega^3-kappa*omega^2+(c*omega[0]*(Q[0]^2*sigma[1]+N[0]^2*sigma[2])/N[0]^2-alpha*N[0])*omega+kappa*omega[0]*(Q[0]^2*sigma[1]+N[0]^2*sigma[2])/N[0]^2 = 0

(3)

``

``

Download d_relation.mw

restart;
alias(u = u(x, y, z, t), f = f(x, y, z, t));
                              u, f
u := (diff(f, y, x, x, x, x))*f^4-(diff(f, x, x, x, x))*(diff(f, y))*f^3-4*(diff(f, y, x, x, x))*(diff(f, x))*f^3+8*(diff(f, x, x, x))*(diff(f, x))*(diff(f, y))*f^2-4*(diff(f, x, x, x))*(diff(f, y, x))*f^3+12*(diff(f, y, x, x))*(diff(f, x))^2*f^2-36*(diff(f, x, x))*(diff(f, x))^2*(diff(f, y))*f+24*(diff(f, x, x))*(diff(f, x))*(diff(f, y, x))*f^2-6*(diff(f, x, x))*(diff(f, y, x, x))*f^3+6*(diff(f, x, x))^2*(diff(f, y))*f^2-24*(diff(f, x))^3*(diff(f, y, x))*f+24*(diff(f, x))^4*(diff(f, y))+2*(diff(f, y, x, t))*f^4-2*(diff(f, y, x))*(diff(f, t))*f^3-2*(diff(f, x, t))*(diff(f, y))*f^3+4*(diff(f, x))*(diff(f, y))*(diff(f, t))*f^2-2*(diff(f, x))*(diff(f, y, t))*f^3-3*(diff(f, z, x, x))*f^4+3*(diff(f, x, x))*(diff(f, z))*f^3+6*(diff(f, x))*(diff(f, z, x))*f^3-6*(diff(f, x))^2*(diff(f, z))*f^2;

Please how  can I use getdata to extract only 200x2 of 208x2 matrix from maple to excel.

For instance I have:

Q:= ( seq( seq( plottools:-getdata(ans1[s1,3])[j,3],j=1..3), s1=1));

It returns  208x2,  200x2 200x2 matrices and I can't extract it using 

`<|>`( seq( seq( plottools:-getdata(ans1[s1,3])[j,3],j=1..3), s1=1));

because the matrices have different dimension. Please, how can make 208x2 matrix to 200x2?

Hi, please how do I delete only the 0*I of the complex expression: 0.*I - 3?

I tried this:  remove(`=`, [0.*I-3, 5, 7, -2], 0.*I);

But couldn't work.

How to linearize eq (6) by neglecting all higher order terms i.e., epsilon[1]^2, epsilon[2]^2, epsilon[1]*epsilon[2]... etc? How to do it in maple?

Calculation.mw

How to find values of w for which determinant of A is zero?

determinant.mw

I want to change the colors of a contour such that the region above zero is represented by a different color instead of yellow color,  and that the colors for conts <0 are one set of colors (maybe yellow to red), and the area for which conts >0 is one color which is very different from the others (so maybe white).

Case4Contour.mw

Hi 

I wanna export data from maple to excel and I tried this:

ExcelTools:-Export(op([1, 1], plots:-display( convert( ans1[12..14,3], list))));

I was expecting 4 columns including z vector and 3 columns from ans1[12..14,3] but got only 2 columns.Case1_Revised_080922.mw

Please any help?

Please, how do I use surface or contour to plot discontinous function?

z = f(x,y) := x-.8193318913*sin(x)*cos(x)/(cos(x)^2+sinh(y)^2)-(0.2931044702e-2*(0.7500000000e-3*x^2-0.7500000000e-3*y^2+0.1622336517e-1))*sin(x)*cos(x)/(((0.7500000000e-3*x^2-0.7500000000e-3*y^2+0.1622336517e-1)^2+0.2250000000e-5*x^2*y^2)*(cos(x)^2+sinh(y)^2))-0.4396567053e-5*x*y*sinh(y)*cosh(y)/(((0.7500000000e-3*x^2-0.7500000000e-3*y^2+0.1622336517e-1)^2+0.2250000000e-5*x^2*y^2)*(cos(x)^2+sinh(y)^2)):

The matrix data[gridji,4] in my code has float undefined. Please, how do I handle it?

1 2 3 4 5 6 7 Last Page 1 of 80