restart;with(geom3d):

point(M,1,2,1):

plane(P,x+y-3=0,[x,y,z]):

plane(Q,2*x-2*y+z-1=0,[x,y,z]):

line(d,[1+2*t,1-t,t],t):

n:=NormalVector(P):

line(Delta,[M,n],t):

Eq:=Equation(Delta):

intersection(N,d,P):

coordinates(N):

decom:=proc(n::posint,m::posint)
> local a,c,d::posint,e,f,i,j,k,l,s,t,u;s:=[];t:={};u:={};
> if m=n then return [[1$m]];fi;
>  for i from 1 to nops(decom(n-1,m)) do
>   a[i]:=op(i,decom(n-1,m));# a is a list
>   c[i]:=[op({op(a[i])})];
>   #d[i,j]:=coun(a[i],op(j,c[i]));#c:unrepeated numbers in a. d:number of repeatance
>    for k from 1 to nops(c[i]) do
>     e[i,k]:=aadd1(c[i],k);

decom:=proc(n::posint,m::posint)
> local a,c,d::posint,e,f,i,j,k,l,s,t,u;s:=[];t:={};u:={};
> if m=n then return [[1$m]];fi;
>  for i from 1 to nops(decom(n-1,m)) do
>   a[i]:=op(i,decom(n-1,m));# a is a list
>   c[i]:=[op({op(a[i])})];
>   #d[i,j]:=coun(a[i],op(j,c[i]));#c:unrepeated numbers in a. d:number of repeatance
>    for k from 1 to nops(c[i]) do
>     e[i,k]:=aadd1(c[i],k);

For example, in the list [1,1,1,1,2,4,6,1,9], how many 1s are there?

Is there a procedure in Maple? Or how can I manage to do it?

See the attached worksheet: myDet.mw

Why does the determinant of A  contain  x, x^3, x^5, x^7 and x^9 ?

Maple15/Linux/Fedora16_x86_64

i am solving for a system of df eq's, but the eq's have these little tricky integrals inside them.....so it either retards operation or no solutoin can be found...They look like...

B_12:=(10*a_2*a_3)*(u/( ((10^2+u)^(3/2)) * ((a_2^2+u)^(3/2)) * ((a_3^2+u)) )):

B_122:=(10*a_2*a_3)*(u/( ((10^2+u)^(3/2)) * ((a_2^2+u)^(5/2)) * ((a_3^2+u)) )):

B_123:=(10*a_2*a_3)*(u/( ((10^2+u)^(3/2)) * ((a_2^2+u)^(3/2)) * ((a_3^2+u)^(3/2)) )):

B_23:=(10*a_2*a_3...

My problem looks as follows................pts1FIXED comes from dsolved solution......it plots a empty figure...and i checked that l(the data pairs looks fine)..

--------------------------------------

Energy:=array(1..450);

Time:=array(1..450);

for i from 1 to 450 do 

	

	Energy[i]:=-1*((32*G*Inertia^2*e^2)/(5*c^2))*pts1FIXED(i);

end do:

l := [[t,Energy[t]] $t=1..450]:

plot(l,style=point,symbol=circle,color=blue);

I need some help...i need a solver using the Runge-Kutta method adaptive size control, to solve a sytem of 5 diff eq's, all of them are interlinked...

And 

Say that dsolve numeric produces a results, called Solution, Now ex  solution(0)=1, and in general solution(t) can be evaluated, now, i want to use this, in E(t)=2*solution(t), in a plot, thus...

plot(E(t),t=0..100,....), but returns the error....

Ok, I am not a very experienced user, but I am mystified by this.  The essence of the problem is that when entering a range of values in an integration, I get a different answer if the range is given as 0..1 to if it is given as 0.0..1 or 0..1.0. Here is the program

### dgig
dgig := proc(psi, chi, lambda, t)
    (psi/chi)^(lambda/2)*t^(lambda - 1)*exp(-(psi*t+chi/t)/2)/\
    (2*BesselK(lambda,sqrt(chi*psi))):
end: ## End of dgig

Suppose I have a file of text that I would like to read into a Maple (16) worksheet as a string. (I have in mind encrypting the string using the RSA encryption scheme.)

If I paste the text directly into the worksheet and surround it by ",  I have no problem, but if I try to read, for example, a file some_text.txt using a command like

     read("C:\\Users\\MyName\\Desktop\\some_text.txt");

I've created a matrix in which each row has two values. I'm using Add(i,i=[matrixName][row element number])/2 to find the mean of each row, which I'd like to transfer to a corresponding element number in an array or a vector, but so far I can only get it to work on single rows or (if I use only Add(i,i=[matrixName]/2), understandably, the mean of the whole matrix. Is there a way to do it wholesale?

Hi! I need help with my project I dont know where to even start!

Question:

Consider a simple polymer , which consist of long chains of simple molecules attached to each other, e.g:

....-CH₂-CH₂-CH₂-...

If the interaction energy between consecutive sections ca be ignored or is low enough, the chain can be bend around at varios places to form varios configurations. One can start by  fixing the first section...

Hi everyone 

Is there a method in maple to plot a rectangular matrix.

Amr

See my code:

----------------

restart;

with(plots); with(PolynomialTools); with(LinearAlgebra);

f := (x, y) -> ((2*x^2+y^2*(x^2*(-2+sqrt(1/x^4+4/(x^2*y^2)+4/y^4+4))-1))/(4*x^2-2*y^2+4))^(1/2) ;

g := diff(f(x, y), x);

plot([g(.1, y), g(1, y), g(3, y)], t = 0 .. 10);

-----------------

Here is error message:

Warning, unable to evaluate the functions to numeric values in the region; see the plotting...