MaplePrimes Questions

Does Maple have something similar to c-means clusteirng in Matlab?

http://www.mathworks.com/help/fuzzy/fcm.html

 How would I go about doing something like this in Maple?

The docs for the package geometry say:

 

The command with(geometry,distance) allows the use of the abbreviated form of this command.

I am new to Maple. What is the abbreviated form of the distance command?

aa3 := Normalize(<<aa2[1,1],aa2[2,1],aa2[3,1]>,<aa2[1,2],aa2[2,2],aa2[3,2]>,<aa2[1,3],aa2[2,3],aa2[3,3]>>, Euclidean);

the result is one column

how to normalize to a matrix form?

I have a problem with finding a interval where parametr "a" is changing.
My funcion:

                 exp(x2(a-x)) for xcR.
I did this, and I don't know what's next.

g:=x->exp(x2(a-x));
dg:=diff(g(x),x);
solve(dg=0,x);
                      0,2/3a
g(0);
g(2/3a);

Hi, all. I need to built a user interface for my program. First, I create two label and two text area for both the value of x and y. After that I assign a button to random generate the value for x and y within the range with x <y. I write the code below into click action of the button.

use DocumentTools in

Do(%x=RandomTools[Generate](integer(range=1..7);

Do(%y=RandomTools[Generate](integer(range=1..7);

while %x>= %y do

Do(%x=RandomTools[Generate](integer(range=1..7);

Do(%y=RandomTools[Generate](integer(range=1..7);

end do;

end use;

This code can generate two random number and insert inside the boxes created. When the button is clicked, two random number generated but the error dialog box shown on the screen with message : Error in component button with caption "generate" cannot determine if this expression is true or false: %y=%x. What I want is only the true value of x<y will be inserted to the text boxes. 

Any suggestion? Any help will be highly appreciated. Thank you.

 

For example,

I have original bit :=[10000000]. after encryption by using technique cryptography, I get thecipher text:=[11110011]. Then I need to measure the total flipped bit in the cipher text with the original bit. Have any method can calculate the total flipped bit? Can someone help me? Thank you.

For example, I havd matrix m:=[0 1 0 1 0 1 0], I want to convert it back to m:=[0,1,0,1,0,1,0]. Can someone help to solve? Thank you.

Hi all, 

I have two inter-related expressions of the form,

(i) T=f(w)

(ii) E=g(w)

where the limits of w in both expressions are same. Now I need to plot a graph T Vs E

Any suggestions???

 

``

restart

``

``

(1)

``

interface(displayprecision = 4)

-1

(2)

(3)

NULL

``

Vert_Coef := Matrix(10, 5, {(1, 1) = 0.332e-1, (1, 2) = 0.187e-1, (1, 3) = 0.104e-1, (1, 4) = 0.52e-2, (1, 5) = 0.21e-2, (2, 1) = 0.509e-1, (2, 2) = 0.301e-1, (2, 3) = 0.166e-1, (2, 4) = 0.83e-2, (2, 5) = 0.31e-2, (3, 1) = 0.571e-1, (3, 2) = 0.364e-1, (3, 3) = 0.218e-1, (3, 4) = 0.104e-1, (3, 5) = 0.42e-2, (4, 1) = 0.582e-1, (4, 2) = 0.405e-1, (4, 3) = 0.260e-1, (4, 4) = 0.125e-1, (4, 5) = 0.42e-2, (5, 1) = 0.582e-1, (5, 2) = 0.416e-1, (5, 3) = 0.270e-1, (5, 4) = 0.135e-1, (5, 5) = 0.52e-2, (6, 1) = 0.582e-1, (6, 2) = 0.416e-1, (6, 3) = 0.281e-1, (6, 4) = 0.145e-1, (6, 5) = 0.52e-2, (7, 1) = 0.582e-1, (7, 2) = 0.416e-1, (7, 3) = 0.291e-1, (7, 4) = 0.166e-1, (7, 5) = 0.62e-2, (8, 1) = 0.582e-1, (8, 2) = 0.416e-1, (8, 3) = 0.301e-1, (8, 4) = 0.177e-1, (8, 5) = 0.62e-2, (9, 1) = 0.582e-1, (9, 2) = 0.416e-1, (9, 3) = 0.301e-1, (9, 4) = 0.187e-1, (9, 5) = 0.73e-2, (10, 1) = 0.582e-1, (10, 2) = 0.416e-1, (10, 3) = 0.301e-1, (10, 4) = 0.197e-1, (10, 5) = 0.83e-2})

Matrix(10, 5, {(1, 1) = 0.332e-1, (1, 2) = 0.187e-1, (1, 3) = 0.104e-1, (1, 4) = 0.52e-2, (1, 5) = 0.21e-2, (2, 1) = 0.509e-1, (2, 2) = 0.301e-1, (2, 3) = 0.166e-1, (2, 4) = 0.83e-2, (2, 5) = 0.31e-2, (3, 1) = 0.571e-1, (3, 2) = 0.364e-1, (3, 3) = 0.218e-1, (3, 4) = 0.104e-1, (3, 5) = 0.42e-2, (4, 1) = 0.582e-1, (4, 2) = 0.405e-1, (4, 3) = 0.260e-1, (4, 4) = 0.125e-1, (4, 5) = 0.42e-2, (5, 1) = 0.582e-1, (5, 2) = 0.416e-1, (5, 3) = 0.270e-1, (5, 4) = 0.135e-1, (5, 5) = 0.52e-2, (6, 1) = 0.582e-1, (6, 2) = 0.416e-1, (6, 3) = 0.281e-1, (6, 4) = 0.145e-1, (6, 5) = 0.52e-2, (7, 1) = 0.582e-1, (7, 2) = 0.416e-1, (7, 3) = 0.291e-1, (7, 4) = 0.166e-1, (7, 5) = 0.62e-2, (8, 1) = 0.582e-1, (8, 2) = 0.416e-1, (8, 3) = 0.301e-1, (8, 4) = 0.177e-1, (8, 5) = 0.62e-2, (9, 1) = 0.582e-1, (9, 2) = 0.416e-1, (9, 3) = 0.301e-1, (9, 4) = 0.187e-1, (9, 5) = 0.73e-2, (10, 1) = 0.582e-1, (10, 2) = 0.416e-1, (10, 3) = 0.301e-1, (10, 4) = 0.197e-1, (10, 5) = 0.83e-2})

(4)

``

``

whattype(Vert_Coef)

Matrix

(5)

``

NULL

``

``

``

``

SR := Vector([10.8696, 21.7391, 32.6087, 43.4783, 46.7391, 54.3478, 65.2174, 76.0870, 86.9565, 100.0000], datatype = float)

Vector(10, {(1) = 10.8696, (2) = 21.7391, (3) = 32.6087, (4) = 43.4783, (5) = 46.7391, (6) = 54.3478, (7) = 65.2174, (8) = 76.0870, (9) = 86.9565, (10) = 100.0000})

(6)

E_G := Vector([250, 500, 1000, 2500, 10000], datatype = float)

Vector(5, {(1) = 250.0000, (2) = 500.0000, (3) = 1000.0000, (4) = 2500.0000, (5) = 10000.0000})

(7)

``

``NULL

NULL

``

NULL

``

B_interp := proc (SR, E_G) options operator, arrow; CurveFitting:-ArrayInterpolation([SR, E_G], Vert_Coef, Array(1 .. 1, 1 .. 1, 1 .. 2, [[[SR, E_G]]]), method = linear) end proc

proc (SR, E_G) options operator, arrow; CurveFitting:-ArrayInterpolation([SR, E_G], Vert_Coef, Array(1 .. 1, 1 .. 1, 1 .. 2, [[[SR, E_G]]]), method = linear) end proc

(8)

``

``

f[18*_1] := B_interp(25.0, 1500.0)

Error, (in CurveFitting:-ArrayInterpolation) invalid input: dimensions of independent data do not match the dimensions of dependent data

 

 

Download Matrix_Interpolation.mw

Hello Mapleprime Community,

Would you be kind to help me with the attached Maple worksheet. My goal is to do a 2D interpolation given the value of SR (slenderness ratio) and E_G (Young's modulus / Shear modulus ratio) in the Matrix table Vert_Coeff with header as vector E_G and the other as vector SR. I’m confused as to how to do the interpolation. Please help.

 

Regards,

 

Moses

So, I have the following inputted into Maple 2015, but the plots are not displaying.  I've double checked the syntax, but I don't see anything wrong.  Any thoughts?  Thank you!

 

plot3d(cos(theta)^2, phi = 0 .. 2*Pi, theta = 0 .. Pi, axes = BOX, style = line, grid = [50, 50], view = [-1 .. 1, -1 .. 1, -1 .. 1], coords = spherical);

 

Trying to plot that here actually displays the plot, but it isn't showing up in my Maple worksheet.

hi

how pdsolve 2D couple  non linear differential equations which attached below?

thanks

2D.mw

PD1 := 12.6000000000000*(diff(U(x, theta), x, x, x, x))-7500*(diff(U(x, theta), x, x))+10.2112755389544*(diff(U(x, theta), x, x, theta, theta))+7.22165554279476*(diff(V(x, theta), x, x, x, theta))-3730.19397871630*(diff(V(x, theta), x, theta))+4.57316679658261*(diff(V(x, theta), x, theta, theta, theta))+36.0000000000000*(diff(W(x, theta), x, x, x))-9375.00*(diff(W(x, theta), x))-15.6731205946742*(diff(W(x, theta), x, theta, theta))-1947.26649812618*(diff(U(x, theta), theta, theta))+1.41357763467822*(diff(U(x, theta), theta, theta, theta, theta)):

PD2 := 3.52500000000000*(diff(V(x, theta), x, x, x, x))-3150.00*(diff(V(x, theta), x, x))+5.05278814097746*(diff(V(x, theta), theta, theta, theta, theta))-4972.65371594662*(diff(V(x, theta), theta, theta))+10.2112755389544*(diff(V(x, theta), x, x, theta, theta))+7.22165554279476*(diff(U(x, theta), x, x, x, theta))+4.57316679658261*(diff(U(x, theta), x, theta, theta, theta))-2536.53190552708*(diff(U(x, theta), x, theta))+97.2834589649212*(diff(W(x, theta), x, x, theta))+30.9917088694028*(diff(W(x, theta), theta, theta, theta))-31273.9463175575*(diff(W(x, theta), theta)):

PD3 := -0.168000000000000e-2*(diff(W(x, theta), x, x, x, x, x, x))+16.9000000000000*(diff(W(x, theta), x, x, x, x))-0.319161728473364e-2*(diff(W(x, theta), x, x, x, x, theta, theta))-37.5000000000000*(diff(W(x, theta), x, x))+796.276013656496*(diff(W(x, theta), x, x, theta, theta))-0.202111525639098e-2*(diff(W(x, theta), x, x, theta, theta, theta, theta))-36.0000000000000*(diff(U(x, theta), x, x, x))+9375.00*(diff(U(x, theta), x))+15.6731205946742*(diff(U(x, theta), x, theta, theta))-97.2834589649212*(diff(V(x, theta), x, x, theta))+31273.9463175575*(diff(V(x, theta), theta))-30.9917088694028*(diff(V(x, theta), theta, theta, theta))+2.06625000000000*10^5*W(x, theta)+313.462411893484*(diff(W(x, theta), theta, theta))+6.77715234781900*(diff(W(x, theta), theta, theta, theta, theta))-0.426628729281504e-3*(diff(W(x, theta), theta, theta, theta, theta, theta, theta))-400/((1-h3*f3(x))*ln(10-10*h3*f3(x))^2)-500/(1-h3*f3(x))^4:

BC := {D[1](U)*(0, theta) = 0, D[1](U)*(1, theta) = 0, D[1](U)*(x, 0) = 0, D[1](U)*(x, 1) = 0, D[1](V)*(0, theta) = 0, D[1](V)*(1, theta) = 0, D[1](V)*(x, 0) = 0, D[1](V)*(x, 1) = 0, D[1](W)*(0, theta) = 0, D[1](W)*(1, theta) = 0, D[1](W)*(x, 0) = 0, D[1](W)*(x, 1) = 0, D[2](W)*(0, theta) = 0, D[2](W)*(1, theta) = 0, D[2](W)*(x, 0) = 0, D[2](W)*(x, 1) = 0, U(0, theta) = 0, U(1, theta) = 0, U(x, 0) = 0, U(x, 1) = 0, V(0, theta) = 0, V(1, theta) = 0, V(x, 0) = 0, V(x, 1) = 0, W(0, theta) = 0, W(1, theta) = 0, W(x, 0) = 0, W(x, 1) = 0}

``

 

Download 2D.mw

I need to receive a plot of the next equation:

E:=0.001;

eq := E*(diff(y(x), `$`(x, 2)))+(x^2+1)*(diff(y(x), x))-x*(x-1/2)^2 = -x^2+2.7^x ;

/*(x=0..1, y(0)=-1, y(1)=0)*/

I try to use DEplot in such way DEplot(eq, y(x), x = 0 .. 1, y = -1 .. 0); but i always get a bug like 

Error, (in DEtools/DEplot) cannot produce plot, non-autonomous DE(s) require initial conditions.

 

Is it possible to do the ab-initio method for differentiation in maple?

In working on an answer to a recent question on MaplePrimes:
http://mapleprimes.com/questions/210948-Can-We-Trust-Maple

I noticed that in solving a "simple" system of two equations in two unknowns, potentially undesirable expansion occurs.
The two equations were only simple in the sense that the solving was very trivial, since the two unknowns occurred (almost) already isolated on the left:
eq1:=diff(x(t),t)=rhs1;
eq2:=1.2345*diff(v(t),t)=rhs2; #rhs2 is very complicated in the link given.
So in isolating I tried:
solve({eq1,eq2},{diff(x(t),t),diff(v(t),t)});
The resulting ode system took considerably longer to solve numerically than an unexpanded version.
It is trivial (of course) to do the solving without solve in this case.

Here is an extremely simple example, where the unknowns are already isolated, so that the equations themselves are actually the solution.

restart;
eq1:=a=b*(c+d);
solve(eq1,a); #No expansion
eq2:=e=b*(c+d)+f*(8+k);
solve({eq1,eq2},{a,e}); #Expanded:
                   {a = b*c+b*d, e = b*c+b*d+f*k+8*f}

Can expansion be avoided?

So, I am working in a partial differeital equations course and we are doing image processing. I am taking an image and pixel by pixel applying the heat equation to smudge out noise. I created a for loop for the equation and it works, I need to itterate it over and over again. Here is my code:

for t from 1 to 20 do    

for j from 2 to 149 do

for k from 2 to 149 do  

`tile2__j,k`:=dellT/(h^(2))*(`tile__j+1,k`-4*`tile__j,k`+`tile__j-1,k`+`tile__j,k-1`)+`tile__j,k`  

 end do

end do  

for m from 2 to 149 do

for n from 2 to 149 do  

`tile__m,n`:= `tile2__m,n`  

end do

end do  

t := t+1  

end do

 

Sorry, I could not get the Maple Math editor to work so I just copied my code. Apologize. I keep getting the error "Error, Invalid loop statement termination." Help?

 

First 1143 1144 1145 1146 1147 1148 1149 Last Page 1145 of 2434