I have a problem with solving these two equations.

When I put this into maple -> solve([0.4=a*cos(b*8)+3 and 0.2=-a*sin(b*8)*b],[a,b])

It say: "Warning, solutions may have been lost". I don't know what to do.

Dear guys,

recently, I met a problem. I need to compute numerically the Mathieu function of bessel I type and bessel K type(Ie,Io,Ke,Ko or denoted by Ce(-q),Se(q),Fek,Gek etc.). That is the mathieu function could be expanded in the series of the modified bessel function I or K. However, I find the commands in maple is very limited, and a question is coming: how to compute these functions with maple?

Thanks for your attention.

yxxu

why do i get this 

"Error, (in solve) cannot solve expressions with diff(diff(y(x), x), x) for x"

when i input these: 

> restart; with(DEtools); with(plots); with(linalg);
> ode2 := x^2*(diff(y(x), `$`(x, 2)))+x*(diff(y(x), x))+4*y(x) = -2*x+7; bc2 := y(1) = 7, y(4) = -1;

a := 1; b := 4

ic2 := bc2[1], (D(y))(a) = alpha

constraint := lhs(bc2[2])-rhs(bc2[2]); constraint = 0

Hello,

I have a big problem With maple, it's a bit complicated, but I will try to detail. it's been months since I still looking for the solution :

I have a profit function expressed as follows:

Hi,

I am trying to implement some numerical algorithm for an optimal control problem (still).

It involves double cycle. In the inner one two DE systems are solved sequentially and the result is used to update values of some objective function.

This part uses remeber option thanks to pagan  and is ok.

However, in the outer cycle I use the values of this objective function to obtain solution for yet another DE and update ICs for the inner cycle.

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