Hello,

I wanted to get the following function integrated with Maple 17:

I used the following command to get a numerical solution for my choise of limits:

Unfortunately, I just get the Integral itself back. Also, taking pi/2 as the upper limit for x does improve anything. I also tried the AllSolutions option but it does not help as well as first trying to integrate only over y and leaving x as it is in order to get at least one integration done.

I hope someone can help. I do not necessarily need to solve it for myself, so if someone gets the result and posts it here, I'm okay with that. Anyway, being able to calculate it myself would be even better.

Hi

I call the MeijerG function in the matlab: MeijerG([[-.3], []], [[.8, 1.3, -.8, -.3, -1.3], []], 1.), and it return an error as follows:

Error, (in evalf/MeijerG) the function is not defined: corresponding GAMMA poles must not coincide.

But when I change the first parameter from -0.3 to -1.3, I can get the result.

So, can you help to explain it? how can I fixed it?

I need to make two graphs.  1st take the equation y= ln3x+3 and graph it and it's derivative.  Then find the equation of the tangent line and the normal line to the given function at x=e   Then I need to create a single graph containing the function, the tangent line, and it's normal line.  Any pointers or help would be greatly appreciated.  I can graph the equation and it's derivative ok but I don't know how to make a graph with multiple stuff on it.  Thanks

Hi there,

I was trying to simulate the behaviour of a one-variable, discrete-time function having three parameters.

The function reads

M[n+1]=(b·theta^(m)·R[n])/(theta^(m)+R[n]^m)

defined for b, theta, m > 0

Say I want to simulate the function for the following values

b := [seq(1 .. 10, 1)];
theta := [seq(1 .. 5, 1)];
m := [seq(1 .. 2, 1)];

for n from 1 to 10.

I guess I need to build a 4-dimension array. But I was not able to find the right way to do this: should I use the Array strucutre? and if yes, how would I do it? As far as I've read, indexing would be an issue: should I create special indexing functions?

If I'm not wrong Maple matrices (Matrix) are just 2D.

Once the values of M computed, I would like to generate the corresponding plots, varying one parameter while the other two are fixed, and drawing the different M's in the same plot.

How can I achieve this?

Furthermore, if I would like to generate all possible combinations, I guess I would need to insert the solution given to the above question within a for loop. Will Maple display all plots or will it overlap/overwrite the preceding plot if used within a loop?

Thanks,

jon

Hi,

I have coded up a vector that is of my interest. The code runs witout any problem and gives me exactly what I want.

newtest.mw

All I want to know, is that if there are more efficient way to do so?

Any tricks, or better use of a particular function that I wasn't aware?

The only tiny bit of unsatisfactory is that, the (1-w) term is at the first term of the addition, is that possible to move it to the last term? Which is more conventional to read.

Thanks,

casper

Hi guys

I have an if statement inside a function. This function then get's called in a sum. However, the result is not correct as the if statement checks wether the argument equals 0, but since the sum passes only the summation index name, it doesn't work as intended.

How can I tell Maple, that it has to evaluate the index before passing it on? MWE appended.

Cheers

NOh

Download MWE.mwMWE.mw

I got a problem in using dsolve.
In my real question, some functions are quite complicated.
so here is a simple example.
F1:=x(t)^2;F2:=piecewise(t>=0,y(t)^3,t>=0.1,exp(y(t)));
eq1:=diff(x(t),t$2)=F1;eq2:=diff(y(t),t$2)=subs(x=y,F1)-F2;
ic1:=x(0)=1.2,D(x)(0)=0;ic2:=y(0)=MM(tf),D(y)(0)=NN(tf);
#tf is the point where x(tf)=30.
dsolve({eq1,eq2,ic1,ic2,x(tf)=30},numeric);

above command can't give an answer.

how to use dsolve solve this problem?

any ideas is appreciated.

Hello Dears

I have this equation

                                          (Napla)^4 * F(x,y) + k^2 *  (Napla)^2 * F(x,y)= 0,      (1)

which may be written as a non-homgeneous Helmholtz equation as

                                          (Napla)^2 * F(x,y) + k^2 * F(x,y)= g(x,y),                (2)

where the function g(x,y) is a harmonic function and (Napla)^2 is the laplace's operator in two dimension.

Can Maple solve equation (1), it will be better. If not may be solve equation (2).

I have an integral that maple can not solve but I can solve it by hand. How can I add this to maple integration database?

f:=int(r^2*BesselJ(0,a*r)*BesselI(1,b*r),r)

Please see the file below.

Integral.mw

Many many thanks! :)

Hi

How can I substitute a function in a matrix? This is what I tried but it's not working at all:

lign1:=x-y+2*z=1:
lign2:=(-2)*x+y+z=0:
lign3:=(-4)*x+y+7*z=2:
lign4:=3*x-2*y+z=1:
T := GenerateMatrix([lign1, lign2, lign3, lign4],[x, y, z], augmented);

Then I'd like to substitute another function with lign1:

lign5:=x+y+z=6:
T2:=subs(lign1=lign5, T);

But T2 = T and that wasn't really the point.

Sorry if the answer is obvious no amount of searching has helped me so far. I hope this question isn't too ridiculous.

Thanks

kappa := Vector(7, [1,w[1]*(1-phi+phi*(1-1/(1+exp(-mu[p]-tau[p3]))))+(1-w[1])*
(1-phi+phi*(1-1/(1+exp(-mu[p]-tau[p3]-eta[p2])))),w[1]*phi/(1+exp(-mu[p]-tau[
p3]))+(1-w[1])*phi/(1+exp(-mu[p]-tau[p3]-eta[p2])),w[1]*(1-phi+phi*(1-1/(1+exp
(-mu[p])))*(1-phi)+phi^2*(1-1/(1+exp(-mu[p])))*(1-1/(1+exp(-mu[p]-tau[p3]))))+
(1-w[1])*(1-phi+phi*(1-1/(1+exp(-mu[p]-eta[p2])))*(1-phi)+phi^2*(1-1/(1+exp(-
mu[p]-eta[p2])))*(1-1/(1+exp(-mu[p]-tau[p3]-eta[p2])))),w[1]*phi^2*(1-1/(1+exp
(-mu[p])))/(1+exp(-mu[p]-tau[p3]))+(1-w[1])*phi^2*(1-1/(1+exp(-mu[p]-eta[p2]))
)/(1+exp(-mu[p]-tau[p3]-eta[p2])),w[1]*(phi/(1+exp(-mu[p]))*(1-phi)+phi^2/(1+
exp(-mu[p]))*(1-1/(1+exp(-mu[p]-tau[p3]))))+(1-w[1])*(phi/(1+exp(-mu[p]-eta[p2
]))*(1-phi)+phi^2/(1+exp(-mu[p]-eta[p2]))*(1-1/(1+exp(-mu[p]-tau[p3]-eta[p2]))
)),w[1]*phi^2/(1+exp(-mu[p]))/(1+exp(-mu[p]-tau[p3]))+(1-w[1])*phi^2/(1+exp(-
mu[p]-eta[p2]))/(1+exp(-mu[p]-tau[p3]-eta[p2]))]);

Download kappa.txt

Here is the expression, I am trying to simplify, given a set of rules. NEW_Cole.mw

I have tried different substitutions, using simplify with side rules, applyrule, eval, subs, algsubs.

But none seem to be working as the way I want them to be.

Is there a better way?

Thanks!

The Embedded Components are containers that currently use industries for modeling complex systems to find viable solutions in real time and thus avoid huge wait times and overload our computer; by this paper should show you how to implement a dynamic worksheet through Embedded Components in Maple; it goes from finding solutions to ordinary differential equations partial; which interact with the researcher using different parameters.
Using graphical programming will find immediate solutions to selected problems in science and engineering criteria of variability and boundary conditions evolving development with buttons on multiple actions.

cimac_2014.pdf

(in spanish)

Solutions_of_Differential_Equations_with_Embedded_Components.mw

Lenin Araujo Castillo

Physics Pure

Computer Science

I am trying the following:

restart; 
with(Physics); 
with(Library); 
Setup(mathematicalnotation = true); 
Setup(op = A);
assume(n::integer);

type(A^n, Physics:-`^`(PhysicsType:-ExtendedQuantumOperator, integer));
type(A^n, PhysicsType:-ExtendedQuantumOperator^integer);

the result is

false
false

However, I expected to get true in both cases. How to construct correct type for A^n?
Thank you.

Is it possible to evaluate a function at multiple points described by an array or something of that sort and have Maple return the evaluations as an array. I need approximations of a function at various values of its argument so it would be nice to do it with a single command.

Thanks

How do I make "simplify" function to recognize complex exponentials. Maple created some exponentiol functions with a square root of negative expression in the denominator of the exponent. It refuses to simplify it. If I manually edit it, by taking negative sign under the root and adding immaginary number, the simplification works. How do I simplify it without manual edit?

The Eq 33 is manually edited line. Simplify function correctly uses cos function.

The Eq 34 is original expression. I tried to tell maple to assume real numbers hoping that it will wotk, but it did not.