Hi,

tmp:=seq([phi[5-j]*(1-p[6-j]),phi[5-j]*p[6-j]],j=1..4);
Vector(8,[seq(op(tmp[i]),i=1..4)]);

"tmp" is basically what I want. Some kind of 'paired' terms, indexed in reverse order.

I wonder if there is a better (perhaps more efficient and "direct") way to do it?

Thanks,

casper

Hi all

 
i am solving a system of four coupled differential equations using bvp[midrich]approach..

The code i am using is in the attach files 

command is like

> A1 := dsolve({bc, eq1, eq2, eq3, eq4}, numeric, output = array([seq(0+0.5e-2*i, i = 0 .. 200*bb)]), method = bvp[midrich]);

with(plots)

 I have a system of 4 rational equations involving 6 variables. I want to solve this system for 4 variables. 

Whenver I try to do this I get either an error message about solutions being lost, or maple being unable to allocate enough memory.

The Commands involved are here: 

Hello,

I want to optimal control, i.e. determine command which minimize an onjective function.

So I want know how I can use the Maple function "optimization[NLPSolve]" in a custom block of Maplesim ?

(I do that in Matlab with "fminsearch" and a Matlab function block but I don't manage to do it with maplesim)

Hello,

I use Maple to solve dynamics of multibody systems. In the mechanical system i study, all my kinematic schemes use the parameter gamma and the project i work on, many people use the parameter gamma.

The issue is that i use the parameter gamma(t) depending of the time and if i use this parameter Maple treat this parameter as the Euler constant gamma and consequently make bad simplifications.

Is it possible to add a code lign to desactivate this...

An indefinite integral cancels out one differentiation:

int(diff(w(s),s,s),s)=diff(w(s),s);

But for a definite integral Maple does not recognize that:

int(diff(w(s),s,s),s=0..L) = int(diff(w(s),s,s),s=0..L);

What I would have expected is something like:

int(diff(w(s),s,s),s=0..L) = eval(diff(w(s),s),s=L)-eval(diff(w(s),s),s=0);

I cannot get Maple to recognize this even when assuming real or positive.

What I am really after is this: w(s...

Good morning sir,

I request you to support me How to generate the maple procedure for QR-Method in linear Algebra?

Good morning sir,

I request you to support me How to generate the iterations for Powermethod?

using Maple commands.

Hi,

While using pdsolve for a coupled system of pdes, the maple results doen't match with the desired ones.

Here is the system

restart:

Eq1:=diff(f(eta,tau),eta$3)+(f(eta,tau)+h(eta,tau))*diff(f(eta,tau),eta$2)-diff(f(eta,tau),eta$1)^2

+(1+epsilon*cos(Pi*tau))*theta(eta,tau)=Omega*diff(f(eta,tau),tau);

Eq2:=diff(h(eta,tau),eta$3)+(f(eta,tau)+h(eta,tau))*diff(h(eta,tau),eta$2)-diff(h(eta,tau),eta$1)^2

+c*(1+epsilon*cos(Pi*tau...

Why can't Maple solve this integral equation using intsolve? May be my error entering it?

sys:= diff(y(t),t)=2-t^2/2 - 1/4 * int(y(t1),t1=0..t);
intsolve(sys,y(t));

Error, (in unknown) invalid input: diff received 0, which is not valid for its 2nd argument

I solve it by hand, and wanted to verify the solution. I get y(t)=C* sin(t/2) -4*t and C can be found to be 12
 

Hello,

I just installed maple 15 as well as the latest update. However, when I was trying to plot, it gives me this:

Which had never happened before I reinstall my operation system.

Could anyone help me please? Thanks a lot!

I have a variable with some assumptions: assume(0<=W<=1).

Now, when I make a list and sort it, everything goes fine:

k := [W, (1/3)*W, (1/3)*W, (1/3)*W]:
sort(k, (a,b)->is(a>b))

I get [(1/3)W, (1/3)W, (1/3)W, W] as expected.

Now I try this:

k:= Vector[column]([W, (1/3)*W, (1/3)*W, (1/3)*W]) 
sort(k, (a,b)->is(a>b))

And Maple gets stuck on this computation. I have tried on both Windows and Linux, on several different computers.