Dear Maple users,

My problem is as follows:

I have a factor base [2,3,5,7,11,33,34,35,36,37,38,39,40]

The numbers from 2 till 11 are primes, the rest is not. 

Then I have to factor (H+c1)(H+c2) in numbers of the factor base , where c1 and c2 go from 1 to some pre-defined limit. H=32 in my case.
And then I have to put the powers of the numbers of the factor base in a matrix. For example: (H+1)(H+1)=33² but also (H+1)(H+1)=3²*11².

That will become in matrix form [0 , 2, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0 ] but also (!) [0 , 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0 ].

This is not what I want! I want no double representations....

What I want is that (H+c1)(H+c2) should be represented in primes in the matrix if possible and else just represented as the other numbers.

hope you guys can help me!

I want to print 2+3= in the input and get exactly the same output.

And how can i do it in a program?

Hello people in mapleprimes,

I want to solve the next system of equation for B/A and C/A.

eq1:=A+B=F+G;
eq2:=k*(A-B)=kappa*(F-G);
eq3:=F*exp(I*kappa*a)+G*exp(-I*kappa*a)=C*exp(I*k*a);
eq4:=kappa*F*exp(I*kappa*a)-kappa*G*exp(-I*kappa*a)=k*C*exp(I*k*a);

But, though it is well-known, solve({eq1,eq2,eq3,eq4},{B/A,C/A})
does not work well, as the values I want to solve it for are
expressions: B/A and C/A not variables.

Then, you might thing the next works well.
eq:=subs({B=A/t,C=A/u},{eq1,eq2,eq3,eq4}):
solve(eq,{t,u});

But, this doesn't work well, with the answer was
only the ratio of t and u expressed as the following:

t = t, u = exp(I*k*a)*(exp(-I*kappa*a)*k^2-exp(I*kappa*a)*k^2-exp(-I*kappa*a)*kappa^2+exp(I*kappa*a)*kappa^2)*t/(4*kappa*k*exp(I*kappa*a)*exp(-I*kappa*a))

Isn't there nice way to solve the above system of equation, except that
sol1:=solve({eq3,eq4},{F,G});assign(sol1);
sol2:=solve({eq1,eq2},{A,B});assign(sol2);

Best wishes
taro

test.mw

1step- I want to integrate the (ME*dF) from t=0 to ∞ .

2step- Evaluate Z=<x,y,z> by solving Z=β*MEint.

Download test.mw

After installing the 18.02 update to Maple 18, the inverse Laplace transform no longer works!

I have a long expression with different order derrivatives, that is written in form like that:

-(D[1](f))(x, y)

I'd like to transform it into standard maple form like:

Is there any special procedure to achieve this goal?

Hi experts

In a procedure with no declared parameteres I would like to return (print) the passed arguments (expecting Matrix structures) in a modified form along with the name of the symbol holding the structure passed as arguments.

That is, the procedure just iterates through the _passed arguments in a for-loop to display the name of the the passed argument (the symbol), a colon and then the modified matrix structure.

My problem is that when looping from i to _npassed arguments, refering to _passed[i] gives me the evaluated form.Tat is what I need to modify the structure but not to list the symbol name.

Say I wanted to print a transposed version of my passed matrices. Then I would call

M:=<<1,2>|<3,4>>;

myProc(M);

And the result I want would be

M: 

But I don't know if _passed holds the symbol names or just the evaluated versions of the passed arguments??

And theoretically the passed argument (assuming a matrix) could be the matrix structure put directly in the procedure call, in which case there is no symbol to refer to.

I hope you get my question and can help me out.

Thanks

Simon

Hello

Please help solve this system

 restart; 
B:=1: 
q:=1*10^3: 
l:=1: 
n:=4.7: 
M_F:=z->2*q*l*(z-l)-q*z^2/2: 
M_1:=z->piecewise((z<l), 2*q*l*(z-l)+l-z-q*z^2/2, (z>l), 2*q*l*(z-l)+l-q*z^2/2-l): 
M_2:=z->2*q*l*(z-l)+2*l-z-q*z^2/2: 
one_int:=z->int(B*(M_F(z)+X_1*M_1(z)+X_2*M_2(z))^n*M_1(z),z=0..2*l); 
two_int:=z->int(B*(M_F(z)+X_1*M_1(z)+X_2*M_2(z))^n*M_2(z),z=0..2*l);
value(one_int(z)); 
value(two_int(z)); 
eqs1:={value(one_int(z))=0,value(two_int(z))=0}; 
fsolve(eqs1);

We received an interesting and timely submission to the Maple Application Center this morning that I think people might be interested in.  It's called:

The Comet 67P/Churyumov-Gerasimenko, Rosetta & Philae, by Dr. Ahmed Baroudy. From the abstract:

Our plan is rather a modest one since all we want is to get , by calculations, specific data concerning the comet and its lander.
We shall take a simplified model and consider the comet as a perfect solid sphere to which we can apply Newton's laws.

We want to find:

I- the acceleration on the comet surface ,
II- its radius,
III- its density,
IV- the velocity of Philae just after the 1st bounce off the comet (it has bounced twice),
V- the time for Philae to reach altitude of 1000 m above the comet .

We shall compare our findings with the already known data to see how close our simplified mathematical model findings are to the duck-shaped comet already known results.
It turned out that our calculations for a sphere shaped comet are very close to the already known data.

Click on the link above if you want to take a look.

eithne

In this document, I think that is the relevent section. But I couldnt understand it.

I want something like this

>plot(sin(x),color=yellow);

>plot(sin(x),color=organe);

>plot(sin(x),color=red);

>plot(sin(x),color=green);

>plot(sin(x),color=ABC); # error

>plot(sin(x),color=CBA);# error

optional, but works fine within a given sets of choices.

How do I write a procedure like that?

Thanks!

st6.mw,st7.mwI want to obtain the analytical solution for this PDE by maple       . But I can not do it. Please help me

((d^2)T/dx^2)+ ((d^2)T/dy^2)+ ((d^2)T/dz^2)=-A*Q(x,y,z)

Where   0 <x<a ,  0 <y<b  ,  0 <z<l

With the boundary conditions:

(dT(0,y,z)/dx)=n-T(0,y,z)    (dT(a,y,z)/dx)=n-T(a,y,z)   (dT(x,0,z)/dy)=n-T(x,0,z)

(dT(x,b,z)/dy)=n-T(x,b,z)     (dT(x,y,0)/dz)=n-T(x,y,0)    (dT(x,y,l)/dz)=n-T(x,y,l)

where n is constant and A is set of parameters.

we have number from 1 to 2014. If we underline the number that divide by 3 and 4, how many pattern of two underlines that next to each other as example?(1,2,3,4,5,6,7,8,9,10,...,2014).

Hello all,

I have the following equation:

N*exp(-(1/2)*eta*epsilon*(N*alpha*epsilon*w+2*N*w*C[max]-alpha*epsilon*z-2*Q1*alpha)/(w*N))*S1*upsilon*w-N*S1*upsilon*w+K1^2*alpha*eta*z*epsilon+K1*alpha*eta*z*epsilon*S1 = 0

in which I need to find solution for epsilon (analytical solution) when epsilon>0.  

Thanks,

Dmitry

I have a nonlinear system with 4 equations and 4 unknowns. I am using fsolve. I know that there are multiple solutions for each variable but am only getting one. I need the others. what do I do??

This is my code:

R__1 := Matrix([[1, 0] , [0, 1] ]);

R__2 := Matrix([[1/2, sqrt(3)/2] , [-sqrt(3)/2, 1/2] ]);

R__3 := Matrix([[-1/2, sqrt(3)/2] , [-sqrt(3)/2, -1/2] ]);

R__4 := Matrix([[-1, 0] , [0, -1] ]);

R__5 := Matrix([[-1/2, -sqrt(3)/2] , [sqrt(3)/2, -1/2] ]);

d__1 := Vector( [ 0, 5.4] );

d__2 := Vector( [ 6.4, 4.539] );

d__3 := Vector( [ 11, 4.078] );

d__4 := Vector( [ 15.5, 2.079] );

d__5 := Vector( [ 19, 1.039] );

a := Vector( [ a__x, a__y] );

A__1:=R__1.a+d__1;

A__2:=R__2.a+d__2;

A__3:=R__3.a+d__3;

A__4:=R__4.a+d__4;

A__5:=R__5.a+d__5;

OO:=Vector([O__x,O__y]);

DA1:=A__2.A__2-A__1.A__1-2*(A__2-A__1).OO;

DA2:=A__3.A__3-A__1.A__1-2*(A__3-A__1).OO;

DA3:=A__4.A__4-A__1.A__1-2*(A__4-A__1).OO;

DA4:=A__5.A__5-A__1.A__1-2*(A__5-A__1).OO;

fsolve({DA1,DA2,DA3,DA4},{a__x,a__y,O__x,O__y});

Thanks for any tips you may be able to offer

How do I plot the following equations in Maple ?:

I already tried this: 

According to the given solution the plot should look like that: 

I think I have to tell maple that the function is defined from R^2 -> R, but I don't know how to do this. 

Thanks in advance for your kind help.