Hi all,

I have three points in 3d space say A1=[a11, a12, a13]; A2=[a21, a22, a23] and A3=[a31, a32, a33]. I want to fill the triangle formed by these points. How can I do that?

Thanks is advance.

Hi all,

I have a 3D graph and when I right click on the surface, under style, I can see the option of contour. When I click on that I can see some contour associated with the surface. I want maple to show the values for each contour when I observe it from the top . I want to copy and paste some surfaces in a same figure and campare the contours.

Thank you

rather than last time saved file result?

When I open maple script file, it display the result which file created.

however, not the latest result 

it need to move cursor to the end of script and press enter again.

I wish to delete the rows that have imaginery components from my results Matrix. Have tried many variants of for loops etc.

Download matrix_imaginery_elements.mw

I have the following command.

with(StringTools);
message := `Kajian ini mempunyai tiga objektif pertama seperti yang ditunjukkan dalam bahagian 1.11. Objektif tersebut harus`;

m := convert(message, bytes);

block := map(convert, m, binary);
block := map2(nprintf, "%08d", block);
block := map(proc (t) options operator, arrow; [seq(parse(convert(t, string)[i]), i = 1 .. length(convert(t, string)))] end proc, block);

block := [[0, 1, 0, 0, 1, 0, 1, 1], [0, 1, 1, 0, 0, 0, 0, 1], [0, 1, 1, 0, 1, 0, 1, 0], [0, 1, 1, 0, 1, 0, 0, 1], ........]

with(Bits);
for i to l do
for j from 3 to 7 do
block[i][j] := 1-block[i][j];  //used to flip bit in between 3rd to 7th bit in a block
end do;
c_block[i] := block[i];
end do;
c_block1 := [seq(c_block[i], i = 1 .. l)];

Error, assigning to a long list, please use Arrays

May i know how to solve this problem? I need to change some bit in a list but receive error when there is more than 100 elements in a list. Thank you.

I have a Document that I have been putting together.  When I insert a Subsection into a section below the title I get a 1D-math input command symbol.  Is there a way to prevent this from happening?  I have not had this problem before. 

It appears there is something confused in the startup of the document because when I start a new document and this doesn't happen.

It's almost like when I insert the subsection it converts that part to a worksheet?????

i'm using maple in a research but i want to add a recursive function h_m(t) in 2 case : if m is integer positive and not, 
la formule est donnée comme suit : 
  if (mod(m,1) = 0  and m>0) then  h:=proc(m,t)  local  t ;  h[0,t]:=t ;   for  i from -4 to  m  by  2 do  h [m,t]:= h[0, t]-(GAMMA(i/(2)))/(2*GAMMA((i+1)/(2)))*cos(Pi*t)*sin(Pi*t)  od:  fi:  end;
and i wanna to know how to programmate a Gaus Hypegeometric function. Thank You

Dear all,

I am trying to solve the following system of equations by using dsolve, but I get the error:  error, (in RootOf) expression independent of, _Z, could you please help me to solve it. Thank you.

restart;
Digits := 20;
with(plots);
Nr := .1; Nb := .3; Nt := .1; Rb := 0; Lb := 1; Le := 10; Pe := 1; ss := .2; aa := .1; bb := .2; cc := .3; nn := 1.5;
Eq1 := nn.(diff(f(eta), eta))^(nn-1).(diff(f(eta), `$`(eta, 2)))-(nn+1)/(2.*nn+1).eta.(diff(theta(eta), eta)-Nr.(diff(h(eta), eta))-Rb.(diff(g(eta), eta))) = 0;
Eq2 := diff(theta(eta), `$`(eta, 2))+nn/(2.*nn+1).f(eta).(diff(theta(eta), eta))+Nb.(diff(theta(eta), eta)).(diff(h(eta), eta))+Nt.((diff(theta(eta), eta))^2) = 0;
Eq3 := diff(h(eta), `$`(eta, 2))+nn/(2.*nn+1).Le.f(eta).(diff(h(eta), eta))+Nt/Nb.(diff(theta(eta), `$`(eta, 2))) = 0;
Eq4 := diff(g(eta), `$`(eta, 2))+nn/(2.*nn+1).Lb.f(eta).(diff(g(eta), eta))-Pe.((diff(g(eta), eta)).(diff(h(eta), eta))+(diff(h(eta), `$`(eta, 2))).g(eta)) = 0;
etainf := 10;
bcs := f(0) = ss/Le.(D(h))(0), theta(0) = lambda+aa.(D(theta))(0), h(0) = lambda+bb.(D(h))(0), g(0) = lambda+cc.(D(g))(0), (D(f))(etainf) = 0, theta(etainf) = 0, h(etainf) = 0, g(etainf) = 0;
dsys := {Eq1, Eq2, Eq3, Eq4, bcs};
dsol := dsolve(dsys, numeric, continuation = lambda, output = procedurelist);
Error, (in RootOf) expression independent of, _Z

hi

how i can solve nonlinear differential equations with shooting method in maple?ω in equation is unknown...

thanks

eq.mw

Download eq.mw

Hi

We know determinant of a square matrix A[ij] (i,j ∈ {1,2,3}) is equal to the following expression

det(A) = 1/6 * e[ijk] * e[pqr] * A[ip] * A[jq] * A[kr] 

in which e[ijk] is a third order Tensor (Permutation notation or Levi-Civita symbol) and has a simple form as follows:

e[mnr] = 1/2 * (m-n) * (n-r) * (r-m).

The (i,j) minor of A, denoted Mij, is the determinant of the (n − 1)×(n − 1) matrix that results from deleting row i and column j of A. The cofactor matrix of A is the n×n matrix C whose (i, j) entry is the (i, j) cofactor of A,

C[ij]= -1 ^i+j * M[ij]

A^-1=C^T/det(A)

The general form of Levi-Civita symbol is as bellow:

I want to write a program for finding inverse of (NxN) matrix:

N=2 →

restart;
N := 2:
with(LinearAlgebra):
f := (1/2)*(sum(sum(sum(sum((m-n)*(p-q)*A[m, p]*A[n, q], q = 1 .. 2), p = 1 .. 2), n = 1 .. 2), m = 1 .. 2)):
A := Matrix(N, N, proc (i, j) options operator, arrow; evalf((37*i^2+j^3)/(2*i+4*j)) end proc):
f/Determinant(A);

N=3 →

restart;
N := 3:
with(LinearAlgebra):
f := (1/24)*(sum(sum(sum(sum(sum(sum((m-n)*(n-r)*(r-m)*(p-q)*(q-z)*(z-p)*A[m, p]*A[n, q]*A[r, z], m = 1 .. N), n = 1 .. N), r = 1 .. N), p = 1 .. N), q = 1 .. N), z = 1 .. N)):
A := Matrix(N, N, proc (i, j) options operator, arrow; 10*i^2/(20*i+j) end proc):
f/Determinant(A);

The results of above programs are equal to 1 and validation of method is observed.

If we can write the general form of determinant then we can find the inverse of any square non-singular matrices.

Now I try to write the mentioned program.

restart;
with(linalg):
N := 7:
Digits := 40:
e := product(product(signum(a[j]-a[i]), j = i+1 .. N), i = 1 .. N):
ML := product(A[a[k], b[k]], k = 1 .. N):
s[0] := e*subs(`$`(a[q] = b[q], q = 1 .. N), e)*ML:
for i to N do
s[i] := sum(sum(s[i-1], a[i] = 1 .. N), b[i] = 1 .. N)
end do:
A := Matrix(N, N, proc (i, j) options operator, arrow; evalf((3*i+j)/(i+2*j)) end proc): # arbitrary matrix
CN:=simplify(s[N]/det(A));

Therefore det(A)= CN^-1 * e[a1,a2,..an] * e [b1,b2,.., bn] * A[a1,b1] * A[a2,b2] * ... * A[an,bn].

The correction coefficient is CN(for N)/CN(for N-1) = N!/(N-1)! =N.

restart:
with(linalg): N := 4: Digits := 20:
e := product(product(signum(a[j]-a[i]), j = i+1 .. N), i = 1 .. N):
ML := product(A[a[k], b[k]], k = 1 .. N):
s[0] := e*subs(`$`(a[q] = b[q], q = 1 .. N), e)*ML:
for r to N do s[r] := sum(sum(s[r-1], a[r] = 1 .. N), b[r] = 1 .. N) end do:
A := Matrix(N, N, proc (i, j) options operator, arrow; evalf((3*i+2*j)/(i+2*j)) end proc):
DET:=S[N]:
for x to N do for y to N do
e := product(product(signum(a[j]-a[i]), j = i+1 .. N-1), i = 1 .. N-1):
ML := product(AA[a[k], b[k]], k = 1 .. N-1):
S[0, x, y] := e*subs(`$`(a[q] = b[q], q = 1 .. N-1), e)*ML:
for r to N-1 do S[r, x, y] := sum(sum(S[r-1, x, y], a[r] = 1 .. N-1), b[r] = 1 .. N-1) end do:
f[y, x] := (-1)^(x+y)*subs(seq(seq(AA[t, u] = delrows(delcols(A, y .. y), x .. x)[t, u], t = 1 .. N-1), u = 1 .. N-1), S[N-1, x, y])
end do: end do:
Matrix(N, N, f)/(DET)*(24/6);
A^(-1);

CN for N=4 and N=3 is 24 and 6 respectively i.e. CN(4)/CN(3)=24/6.

When I use bellow procedure the error "(in S) bad index into Matrix" is occurred.

Please help me to write this algorithm by using procedure.

Thank you 

restart; with(linalg): Digits := 40: n := 7:
S := proc (N) local e, ML, s, i:
e := product(product(signum(a[j]-a[i]), j = i+1 .. N), i = 1 .. N):
ML := product(A[a[k], b[k]], k = 1 .. N):
s[0] := e*subs(`$`(a[q] = b[q], q = 1 .. N), e)*ML:
for i to N do s[i] := sum(sum(s[i-1], a[i] = 1 .. N), b[i] = 1 .. N) end do
end proc:
A := Matrix(n, n, proc (i, j) options operator, arrow; evalf((3*i+j)/(i+2*j)) end proc): # arbitrary matrix
CN := simplify(S(n)/det(A))

with(Statistics):
X := RandomVariable(Normal(0, 1))

DensityPlot(X,filled=true)

I don't know why the plot doesn't produce a shaded plot.

Hi all,

I seem to be quite stuck on figuring out how to leave certain letters (e.g. planck's constant h) inside the equation without having to assign it as some particular number. 

What I am trying to do is find the value of a when the following equation is at a minimum:

E = (a*(h^2)/2m) + 0.3989422804/sqrt(a)

Here h and m are what I want to set as constants without actually setting them to h := 1 because I want a in terms of h and m. I have already found the derivative dE/da:

((h^2)/2m) - 0.1994711402/a^(3/2)

But I cannot use fsolve to find the value of a at the minimum because it keeps saying that h and m are variables and unsolved for.

Any help would be greatly appreciated.

How to solve following recurrence equation:

a(0)=2;

a(n+1)=a(n)+a(n)^2

I tried,but it doesn't work.

How to find the sequence a_{n ?}

Maple_worsheet.mw

Mariusz Iwaniuk

Below is a custom distribution created based on a function that takes a parameter.

It is possible to create the custom distribution e.g. as D1 and then use it afterwards to find e.g. Mean, but it is not possible to call Mean directly with the creation of the distribution in the call.

Why is that ?

Hi!!

I am trying to plot the above curve:

restart; with(plots)

> f0 := proc (t) options operator, arrow; t, (-1)*3.9*t*(t-1) end proc; 

> IFS := proc (i, x, y) if i = 1 then return (1/2)*y, (1/2)*x end if; if i = 2 then return (1/2)*x, (1/2)*y+1/2 end if; if i = 3 then return (1/2)*x+1/2, (1/2)*y+1/2 end if; if i = 4 then return -(1/2)*y+1, -(1/2)*x+1/2 end if end proc; 

> g := proc (t) local j; for j to 4 do if evalf((1/4)*j-1/4) <= evalf(t) and evalf(t) <= evalf((1/4)*j) then return IFS(j, f0(4*t-j+1)); break end if end do end proc; 

Thus, the instruction  parametricplot(['g'(t),t=0..1]) return the message  

Error, (in plot) incorrect first argument [g(t), t = 0 .. 1]

Some idea or hit to plot this?

Thank you for your time