MaplePrimes Questions

Hello everyone !

I have a problem when I want to calculate the following multiple integration numerically:

>evalf(Int(exp(sum(x[i],i=1..6)^2),[seq(x[i]=-1..1,i=1..6)]));
  value(%);

It doesn't work. But when I replace sum(x[i],i=1..6)^2 with sum(x[i],i=1..6), it works. Is there any feasible solution to my problem ?

Thank you for reading !

 

Dear Friends

In differential expressions(See Maple file) how to find coefficiets of dependent variable "u(x,t)" and "v(x,t)" and of their differentials ? There is command "dcoeffs(function)", but that work for single dependent variable only but in our case there are two dependent variables in consideration. There are other options like "indets", "specindex" but those do not work.

 


with(PDEtools):

DepVars; -1; [u(x, t), v(x, t), r[1](t), r[2](t), s[1](t), s[2](t), p[1](t), p[2](t), alpha[1](x, t), beta[1](x, t), beta[2](x, t), delta[1](x, t), delta[2](x, t)]

[u(x, t), v(x, t), r[1](t), r[2](t), s[1](t), s[2](t), p[1](t), p[2](t), alpha[1](x, t), beta[1](x, t), beta[2](x, t), delta[1](x, t), delta[2](x, t)]

(1)

alias(u = u(x, t), v = v(x, t), r[1] = r[1](t), r[2] = r[2](t), s[1] = s[1](t), s[2] = s[2](t), p[1] = p[1](t), p[2] = p[2](t), alpha[1] = alpha[1](x, t), beta[1] = beta[1](x, t), beta[2] = beta[2](x, t), delta[1] = delta[1](x, t), delta[2] = delta[2](x, t))

u, v, r[1], r[2], s[1], s[2], p[1], p[2], alpha[1], beta[1], beta[2], delta[1], delta[2]

(2)

(diff(r[1], t))*(-s[1]*u*(diff(u, x))-p[1]*((diff(u, x))*v+u*(diff(v, x)))-alpha[1]*(diff(u, x))-beta[1]*u-delta[1])/r[1]+r[1]*(diff(alpha[1]*(diff(u, x))+beta[1]*u+delta[1], x, x))+(diff(s[1], t))*u*(diff(u, x))+s[1]*(alpha[1]*(diff(u, x))+beta[1]*u+delta[1])*(diff(u, x))+s[1]*u*(diff(alpha[1]*(diff(u, x))+beta[1]*u+delta[1], x))+(diff(p[1], t))*((diff(u, x))*v+u*(diff(v, x)))+p[1]*((diff(alpha[1]*(diff(u, x))+beta[1]*u+delta[1], x))*v+(diff(u, x))*(alpha[1]*(diff(v, x))+beta[2]*v+delta[2])+(alpha[1]*(diff(u, x))+beta[1]*u+delta[1])*(diff(v, x))+u*(diff(alpha[1]*(diff(v, x))+beta[2]*v+delta[2], x)))+(diff(alpha[1], t))*(diff(u, x))+alpha[1]*(diff(alpha[1]*(diff(u, x))+beta[1]*u+delta[1], x))+(diff(beta[1], t))*u+beta[1]*(alpha[1]*(diff(u, x))+beta[1]*u+delta[1])+diff(delta[1], t)

(diff(r[1], t))*(-s[1]*u*(diff(u, x))-p[1]*((diff(u, x))*v+u*(diff(v, x)))-alpha[1]*(diff(u, x))-beta[1]*u-delta[1])/r[1]+r[1]*((diff(diff(alpha[1], x), x))*(diff(u, x))+2*(diff(alpha[1], x))*(diff(diff(u, x), x))+alpha[1]*(diff(diff(diff(u, x), x), x))+(diff(diff(beta[1], x), x))*u+2*(diff(beta[1], x))*(diff(u, x))+beta[1]*(diff(diff(u, x), x))+diff(diff(delta[1], x), x))+(diff(s[1], t))*u*(diff(u, x))+s[1]*(alpha[1]*(diff(u, x))+beta[1]*u+delta[1])*(diff(u, x))+s[1]*u*((diff(alpha[1], x))*(diff(u, x))+alpha[1]*(diff(diff(u, x), x))+(diff(beta[1], x))*u+beta[1]*(diff(u, x))+diff(delta[1], x))+(diff(p[1], t))*((diff(u, x))*v+u*(diff(v, x)))+p[1]*(((diff(alpha[1], x))*(diff(u, x))+alpha[1]*(diff(diff(u, x), x))+(diff(beta[1], x))*u+beta[1]*(diff(u, x))+diff(delta[1], x))*v+(diff(u, x))*(alpha[1]*(diff(v, x))+beta[2]*v+delta[2])+(alpha[1]*(diff(u, x))+beta[1]*u+delta[1])*(diff(v, x))+u*((diff(alpha[1], x))*(diff(v, x))+alpha[1]*(diff(diff(v, x), x))+(diff(beta[2], x))*v+beta[2]*(diff(v, x))+diff(delta[2], x)))+(diff(alpha[1], t))*(diff(u, x))+alpha[1]*((diff(alpha[1], x))*(diff(u, x))+alpha[1]*(diff(diff(u, x), x))+(diff(beta[1], x))*u+beta[1]*(diff(u, x))+diff(delta[1], x))+(diff(beta[1], t))*u+beta[1]*(alpha[1]*(diff(u, x))+beta[1]*u+delta[1])+diff(delta[1], t)

(3)

In above differential expressions how to find coefficiets of dependent variable "u(x,t)" and "v(x,t)" and of their differentials ? There is command "dcoeffs(expr,u(x,t))", but that work for single dependent variable only but in our case there are two dependent variables in consideration. There are other options like "indets", "specindex" but those do not work.

``


Download Coefficients_in_differential_expression.mw

Hello

I have an SEIR model.

Equation 5 is for disease death but I would like to plot the cumulative numbers of disease death which will be the integral of Equation 5. I added the integral inside odeplot but it is not working. Any idea  about  how to compute the integral ?

Maple code is attached

Thank you

code.mw

Found a strange behaviour in Mapke 2015 of the sqrt-function after loading the GRTensor package:

the square-root of a non-square integer, e.g. sqrt(5), does not terminate. 5^(1/2) instead works fine.

Can be reproduced with Maple 18, but not with Maple 11.

I consider this a serious bug, as it makes any expressions containing such roots useless.

As it worked with Maple 11 I am inclined to see it as your fault.

 

Let's say I've got somehow a result of the product:

exp(k*t)*cos(t*omega)

as series expression and I want to factor out the series of

exp(k*t)

in order to get

cos(t*omega)

How can I do that with Maple?

In PDEtools, suppose to I wish assign zero value to certain first order partial derivative such that higher order derivatives automatically vanish in subsequent excutions, how I can do that?

 

with(PDEtools):

alias(u = u(x, y, t))

u

(1)

Suppose we wish following derivative equal to zero,

diff(u, x) = 0

diff(u, x) = 0

(2)

If we use ":=" for value assignment we will get error. Under above assuption how can we make following derivatives zero?

diff(u, x, x); 1; diff(u, x, y); 1; diff(u, x, x, x, y)

diff(diff(u, x), x)

 

diff(diff(u, x), y)

 

diff(diff(diff(diff(u, x), x), x), y)

(3)

``

 

Download Assiging_Derivative.mw

I tried to load my document containing some notes, but then I got the message "There were problems during the loading process, Your worksheet may become incomplete", and as the message said my worksheet were incomplete. Is there a way to restore the document? I have tried following this and added the line it suggested:

http://www.maplesoft.com/support/help/Maple/view.aspx?path=worksheetmaybeincomplete

But it didn't work.

I have attached the file.

Noter.mw

I never expected that the reflected light direction of sun from moon in the sky would be so dificult to imagine ...

at the following article mentioned :

we derive an equation for the magnitude of the moon tilt illusion that can be applied to all con gurations of sun and moon in the sky.

THE MOON TILT ILLUSION

 

since the calculations contains many steps and high level mathematical formula , there is no way rather to recourse to maple (powerful math assistant )

I hope there was adaptations between a lots of functions and predefined schema of maple and this problem so that the calculations and visualization facilitated several times ?

your effort will be a graet present for all the people of the world that look to the moon crescent everytime !

 

hi.i trust that attached equation has more answer but fsolve only gain some of them!!! how i can gain another that i know value of them?

another root  that i known, are : 0.165237712988657e-1    and     .103583272213766    and    .290071279318035

thanks 

root.mw

Hi!

I am comptuing the eigenvalues and the characteristic polynomial of a 8 by 8 symmetric matrix, say M. Thus, we define the matrix M, and compute its charast. plynm. by

 

 

and its eigenvalues with the command

 

 

Well, Maple returns the charast. polynm. an dthe eigenvalues. But, if we compute p(E[k]), for k=1,...,8, thats is, the values of the polynomial p(x) in the eingenvalues, Maple not turns cero!!! I'm really confused ... anyone know what could be happening?

 

Maple attached file with this example. Thank very much for your help!!

 

Download exam_eigenvalues.mw

hi .how i can calculate eigenvector associated with the eigenvalue of the matrix.for example according attached file below

what are  eigenvector associated with the eigenvalue of matrix q which  determined as (2646.408147, 3142.030259, 6621.757707) respectively??

thanks...

eign.mw 

I would like to learn to use Maple to develop applications of math, physics, astronomy ecc.. . The problem is that I do not know where to start. Could you help me ? Thank you.

Blustar

I have some preferences for viewing. Specially, I prefer atomic variables to be showed with a different color. So I check mark the related option for this in the view menu. But every time I close and reopen the Maple, it just restore my setting for atomic variables. Suprisingly, this does not happen for other viewing options in the view menu! What should I do? Why is it like this?

plots[implicitplot3d](max(-x+y+z, x-y+z, x+y-z) = 1.0, x = 0 .. 1, y = 0 .. 1, z = 0 .. 1);

The help page for max does not explain or show an example of max(sequence of expressions)= a constant. 

hi.after calculate Determinant of matrix  and gain value omega'' ω'' by fsolve rule ,when substuting result (ω) in matrix (q) and calculate Determinant again, this value is not zero!!!! may i use LUDecomposition?determinan.mw

First 1174 1175 1176 1177 1178 1179 1180 Last Page 1176 of 2429