Hi,

   The error showed up whenever I tried to do the Maximum Likelihood Estimation with Global Optimization Toolbox.My objective (likelihood ) function f is defined as follows(A1 to A5 are intermediate variables; parameters  are x[1] to x[12]）:

   f := 0;

   for j from 1 by 1 to 21 do

       A1 := exp(x[1]+x[2]*c[j]);

       A2 := exp(x[3]+x[4]*c[j]);

How do I input the integral?

int((x-m^2)/(x+1),x=0..m^2).

Please help me. Thank you.

I'm using Maple to simulate particle collision and interaction (absorption and scattering) and keep tracking its energy until each particle either get absorbed or leaked (Energy less than predefined minimum energy). I also use Maple statistical pacakge. I also need to call MATLAB link to generate histogram in a logarithmically spaced bin, a feature that I could not find anywhere in Maple and I think it is one of the reasons for calculation performance issue.

The problem...

I am new to Maple so apppologies if this is a very basic question. I have a dual sum loop. Sigma i=1to 6(sigma j=1to 6)A_iB_j.

That works fine. I need to include the condition i+j<=8. So when when i=3 j can only count to 5, i=4 j count to 4, i=5 j count to 3. etc.

At present I am just using the work sheet. Haven't got into writing code yet.

Would appreciate any assistance,

I am using a link with Matlab and I want to call upon matlab to do determinants of matrices for me in a procedure.

For example say i wanted to get the determinant of

a:=([1,z],[1,z^2]), when I put this into maple with Matlab[det](a) I get error in Matlab:setvar unable to store 'Re(z)' when data type=float[8].

Any idea how to fix this error?

Any help would be much appreciated, thanks!

Hi All,

I am using maple 15 for a mechanical modeling, which has some very complex trigonomitry manipulations. The

script runs fine in maple 9.5.1 under window XP. But when I try to load the script in maple 15 under window 7,

it keeps loading forever. I noticed the memory usage in window manager to approach 1GB. (I have 4 GB in my PC).

The 9.5.1 uses much small portion of the memory. Does anyone see this problem before?

Hi,

I am new to Maple and asked a question a few days a go and got a great answer (thanks for that!). I am getting more familiar with the software.

However, I am struggling to plot a 'stability diagram' with the Mathieu solutions. 

In an attempt to solve the equations from (link):

http://www...

I have some data that requires a logarithmic y-axis for its histogram, but I want the y-axis to start at y=10^0. In my histogram, the y-axis starts near 10^-2, seemingly because I have 0-points in my data. How can I restrict the y-axis to start at 10^0? 

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,

I have been trying, without success, to reproduce the plot consisting of the family of curves for the equation v(t), as can be seen on this web page : http://tutorial.math.lamar.edu/Classes/DE/Linear.aspx.

Any help appreciated.

thanks,
Andrew

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

Hi,

Can someone please tell me why I am getting this error (in bold):

> ODE := diff(y(x), x)+3*y(x) = x^2
                                     
> ics := y(0) = 1

> dsolve({ode,ics})

Error, (in dsolve) not an ODE system with respect to the unknowns {y(0)}

restart; with(plots); Eq1 := diff(f(eta), `$`(eta, 3))-(diff(f(eta), eta))^2+m*f(eta)*(diff(f(eta), `$`(eta, 2)))-(M^2+lambda*P)*(diff(f(eta), eta)) = 0; Eq2 := diff(theta(eta), `$`(eta, 2))+m*pr*f(eta)*(diff(theta(eta), eta))-lambda1*pr*theta(eta)*(diff(f(eta), eta)) = 0; bcs1 := f(0) = s, (D(f))(0) = -1, (D(f))(10) = 1; pr := .72; M := 1; m := 2; lambda := 1; lambda1 := 1; bcs2 := theta(10) = 0, theta(0) = 1; L := [0.5e-1, .1, .2, .4, .6, .8, 1, 5, 10, 20]; for k to 10 do R := dsolve...