program is：
for t from 10 by 2 to 100 do h := solve({x*y*z = 6*t^3, x-y-z = 0, x+y+z = 6*t}, {x, y, z}); A[(t-8)*(1/2)...

Hi all,

I'm currently running a parameter space search via Maple scripting. The pseudo-code goes something like this:

fopen text file to write data to;

initialize parameters (variable) x, y, and z;

for x number of iterations

   A:-Simulate(model.msim, with [parameters]);

   Update parameters with new values;

   Write to text file;

end do;

fclose(text file)

I use Maple 15 and I want to take the real and the imaginary part of a simple expression e.g. 5*Dirac(x)+3*I, having already assumed x a real value. However, Maple seems to have problem with the Dirac function. The output looks like this:

> assume(x, 'real');
> Re(5*Dirac(x)+3*I);
                         5 Re(Dirac(x))
> Im(5*Dirac(x)+3*I);
                       3 + 5 Im(Dirac(x))

Any thoughts on how to overcome this are welcomed.

Hi I'm new to maple and have to solve this:

Find the transferfunction

I beleve the answer should be:
H(s) = Y(s)/G(s)  =  1/(s^2+(1/4)s+1)   (correct me if I'm wrong)
but i have to show how to do it in Maple in a scool project

please help me have tried for several hours now, but cant get it right

I want to solve equations,program in maple as follows:

solve({n = b+2*c+q+2*r, (a+b+c)/(p+q+r) = 9, n*b-3.32*a*(a+p+2*n)*10^(-7)*t^(5/2)*exp(-18.27*10^4/t) = 0, n*c-3.32*a*(a+p+2*n)*10^(-7)*t^(5/2)*exp(-32.035*10^4/t) = 0, n*q-3.32*p*(a+p+2*n)*10^(-7)*t^(5/2)*exp(-9.16*10^4/t) = 0, n*r-3.32*p*(a+p+2*n)*10^(-7)*t^(5/2)*exp(-18.77*10^4/t) = 0, a+b+c+p+q+r+n-.73*10^28/t = 0}, {a, b, c, n, p, q, r})

I can't get the results but  "Warning, solutions may have been lost"

Can anyone help me write a procedure for Maple 15 that does the following:

INPUT: number of unknowns and equations n; augmented matrix A=[aij], where 1<=i<=n and 1<=j<=n+1.
OUTPUT: solution x[1], x[2], ... , x[n] or message that the linear system has no unique solution.

Step 1: For i=1,...,n-1 do steps 2-4.  # Elimination process.

Step 2: Let p be the smallest integer with i<=p<=n and A[p,i]≠0.
          ...

I have this

plotS := plots[:-implicitplot3d](S0(2, 1, x0, sigmaH) = S0(c, delta, xT, sigmaL), c = 0 .. 4, delta = 1 .. 3, xT = .8 .. 1)

to draw the graph below. What I would really like is to revert the direction of the xT axis. Any way to do this??

Thank you for your help

I want to write a code to generate Legendre Polynomials using the Gram-Schmidt Process. I have no experience with programing. Here is how the process should work:

Given "(f) = {f[0]=1, f[1]=x , f[2]=x^(2) , f[3]=x^(3) , ... , f[n]=x^(n)} "
       
Find the Legendre Polynomials: "{P[0] , P[1] , P[2] , P[3] , ... , P[n]}"

First,  P[0]=f[0]=1

then, the rest are defined recursively by

P[1]=f[1]-(()/())*P[0] = x

I have a problem with animations using a density plot as background.

The density plot covers the animation.

Here is an example:

plt:=densityplot(-(x^2+y^2),x=-1..1,y=-1..1);
animate(plot,[[sin(r)*cos(t),sin(r)*sin(t),t=0..2*Pi],thickness=10],r=0..Pi,background=plt);

Can this be solved? 

Best regards

  anthei 

I am unable to find where the problem is within a plot. Here is the error message: Error, (in plot) incorrect first argument [(k-1)*r*(cos(t)+cos(k-1)*t/(k-1)), (k-1)*r*(sin*t-sin(k-1)*t/(k-1)), t = 1 .. 1, r = 1 .. 1, k = 3 .. 3]

Is there a symbol or something that i am supposed to be looking for in the error message so i know where the error is occuring?

Firstly applogies for if this question has been asked before.

I am trying to plot a simple epidemic model for where I can change the rates.

My equations are:

de1:= diff(X(t),t)=a*X(t).Y(t);

de2:=diff(Y(t),t)=a*X(t).Y(t) - b*Y(t);

where X is susceptibles and Y is infected.

I would like to be able to produce a graph with the susceptible and infected plotted. Where i can change the initial conditions of X(0) and Y(0) using sliders...

I have two inequalities in the complex variable z, called "ineq" and "ineq2".  When I issue the following plotting command, I can show the region of the complex plane that satisfies both inequalities:  It appears as a black region:

plots[implicitplot](eval(ineq = ineq2, z = x+I*y), x = -8 .. 8, y = -8 .. 8, gridrefine = 3, filled = true)

Now, I would like to generate a similar overlap plot, but for three inequalities,...

Hi,

I just noticed something rather peculiar (well, to me, at least). When I compute the limit

limit((1-1/x)^(-x),x=-infinity)

as is, I get the expected result of exp(1).

However, if I load the RealDomain package prior to computing the limit,

with(RealDomain):

limit((1-1/x)^(-x),x=-infinity)

I get that the limit is undefined. 

Any ideas as to why that is the case? Thanks!

As an intermediate result in a longer calculation I got the following equation:

Isolating p, Maple delivers

Calculated manually,...

Dear all,

Maple gives me the complete solution to the attached PDE. This solution is expressed in terms of KummerU and KummerM. What I want is to obtain the solution for real numbers only, i.e. I want to get rid of the possible imaginary numbers in the solution. This implies that the separation constant, _c1, must be less or equal to 2r. 

How can I do this?

Furthermore, the four arbitrary constants will need to be determined through boundary conditions....