@Rouben Rostamian  

these are randomly pick from combinations, no meaning. I am just writing scripts to select equations from all combinations. Thanks. Any common initial and boundary conditions. Suggest?

i want build my own version Calabria yau equation and apply hodge deformation

@nm 

follow example n := 7;

sol is empty. how to use this function?

@Carl Love 

thank you , i know how to do now.

@Markiyan Hirnyk 

in page 137 of an introduction to groebner bases

how to elimininate the redundant solution (y^2-x*z, 0, -x^2+y*w)

from 3 of them?

T := lexdeg([cc1,cc2],[x,y,z]);
intermsof1 := y;
intermsof2 := -z;
GB := Basis([cc1-intermsof1, cc2-intermsof2],T);
result := NormalForm(y^2-x*z, GB, T);

T := lexdeg([cc1,cc2],[x,y,z]);
intermsof1 := -x;
intermsof2 := y;
GB := Basis([cc1-intermsof1, cc2-intermsof2],T);
result := NormalForm(0, GB, T);

T := lexdeg([cc1,cc2],[x,y,z]);
intermsof1 := w;
intermsof2 := -x;
GB := Basis([cc1-intermsof1, cc2-intermsof2],T);
result := NormalForm(-x^2+y*w, GB, T);

@Markiyan Hirnyk

when i try thrid time again,

e1 is correct, it is y,

but e2 is wrong, it should be x in e2 + x, but it return e2 + z

it is wrong in maple 18.

with(Groebner); T := lexdeg([x, y, z], [e1, e2]); intermsof1 := y; intermsof2 := -z; GB := Basis([e1-intermsof1, e2-intermsof2], T); result := NormalForm(-x*z+y^2, GB, T)

[e2 + z, y - e1]
   2
e1 + e2 x

@Markiyan Hirnyk 

after discover my code , reuse e1 e2 e3 again, i rewrite into another variables

and search my previous post.

i try again, no expected coefficient return

T := lexdeg([cc1,cc2],[x,y,z]);
intermsof1 := y;
intermsof2 := -z;
GB := Basis([cc1-intermsof1, cc2-intermsof2],T);
result := NormalForm(y^2-x*z, GB, T);

[cc2 + z, cc1 - y]
f

with(Groebner):
T := lexdeg([cc1],[x,y,z]);
intermsof1 := y;
intermsof2 := -z;
GB := Basis({expand(cc1*intermsof1), expand((1-cc1)*intermsof2)}, T);

[cc2 + z, cc1 - y]
  2
y - x z

@Markiyan Hirnyk 

https://onedrive.live.com/redir?resid=E0ED7271C68BE47C!1016&authkey=!AHon_gXBNyNJsn8&v=3&ithint=photo%2cpng

expect [y,x]

y^2-x*z = y*(y) - z*(x)

but it can not return correct result, where is wrong?

is it a bug?

after change, still got error

> with(Groebner); T := lexdeg([x, y, z], [e1, e2]); intermsof1 := y; intermsof2 := -z; GB := Basis([e1-intermsof1, e2-intermsof2], 'tord', T); result := NormalForm(y^2-x*z, GB, 'tord', T); result := NormalForm(y^2-x*z, GB, T);
lexdeg([x, y, z],

[Vector[column](%id = 87065248), Vector[column](%id = 84467752)])
y
-z
[y + Vector[column](%id = 71871288), Vector[column](%id = 72754708)

+ Vector[column](%id = 85921608) + Vector[column](%id = 86610144) + z]
Error, invalid input: Groebner:-NormalForm expects its 3rd argument, tord, to be of type Or(ShortMonomialOrder, MonomialOrder), but received tord
Error, invalid input: Groebner:-NormalForm expects its 3rd argument, tord, to be of type Or(ShortMonomialOrder, MonomialOrder), but received lexdeg([x, y, z], [Vector(3, {(1) = 1, (2) = 0, (3) = 0}), Vector(3, {(1) = 0, (2) = 1, (3) = 0})])
> result := NormalForm(y^2-x*z, GB, tord, T);
Error, (in Groebner:-NormalForm) invalid argument for sign, lcoeff or tcoeff

@Markiyan Hirnyk 

What surprise about asking questions? please go back to the question and my question do not have behavior problem.

i just had paranoia for correct way.

@Markiyan Hirnyk 

all my gmail accounts are banned to ask questions in stackoverflow and math stackexchange.

Please understand my situation below

Sorry, we are no longer accepting questions from this account. See the Help Center to learn more.

@tomleslie 

 i succeed 

to prevent others who pick matured fruits and edit my code and then 

say that i can not do privately,

i show my correct code.

@tomleslie 

if input matrix is Matrix([[1,3,3],[3,3,3],[3,3,3]]);

output will be

m := Matrix([[a,a*b*c,a*b*c],[a*b*c,a*b*c,a*b*c],[a*b*c,a*b*c,a*b*c]]);
m := Matrix([[b,a*b*c,a*b*c],[a*b*c,a*b*c,a*b*c],[a*b*c,a*b*c,a*b*c]]);
m := Matrix([[c,a*b*c,a*b*c],[a*b*c,a*b*c,a*b*c],[a*b*c,a*b*c,a*b*c]]);

if input matrix is Matrix([[1,2,3],[3,3,3],[3,3,3]]);

output will be

m := Matrix([[a,a*b,a*b*c],[a*b*c,a*b*c,a*b*c],[a*b*c,a*b*c,a*b*c]]);
m := Matrix([[a,a*c,a*b*c],[a*b*c,a*b*c,a*b*c],[a*b*c,a*b*c,a*b*c]]);
m := Matrix([[a,b*c,a*b*c],[a*b*c,a*b*c,a*b*c],[a*b*c,a*b*c,a*b*c]]);

m := Matrix([[b,a*b,a*b*c],[a*b*c,a*b*c,a*b*c],[a*b*c,a*b*c,a*b*c]]);
m := Matrix([[b,a*c,a*b*c],[a*b*c,a*b*c,a*b*c],[a*b*c,a*b*c,a*b*c]]);
m := Matrix([[b,b*c,a*b*c],[a*b*c,a*b*c,a*b*c],[a*b*c,a*b*c,a*b*c]]);

m := Matrix([[c,a*b,a*b*c],[a*b*c,a*b*c,a*b*c],[a*b*c,a*b*c,a*b*c]]);
m := Matrix([[c,a*c,a*b*c],[a*b*c,a*b*c,a*b*c],[a*b*c,a*b*c,a*b*c]]);
m := Matrix([[c,b*c,a*b*c],[a*b*c,a*b*c,a*b*c],[a*b*c,a*b*c,a*b*c]]);

restart;
G := (L::list(name),n::posint)->op((`*`@op)~(combinat[choose](L,n))):
G( [q,r,s], 3 );
with(LinearAlgebra):
aa := Matrix([[1,2,3],[3,3,3],[3,3,3]]);
num1 := 1;
for i from 1 to 3 do
for j from 1 to 3 do
num1 := num1 * nops([G( [a,b,c], aa[i,j] )]);
od:
od:
numcol := nops(aa);
m := Matrix([[0,0,0],[0,0,0],[0,0,0]]);
delta := proc(m,bb,ii2,jj2,poin)
finish := 0;
for ii from 1 to numcol do
for jj from 1 to numcol do
num1 := num1 - nops(alist);
if poin <= 3 and finish = 0 and ii = ii2 and jj = jj2 then
m[ii,jj]:=bb[poin];
finish := 1;
end if:
od:
od:
return m;
end proc:
initit := proc(m,aa)
for ii from 1 to numcol do
for jj from 1 to numcol do
m[ii,jj] := [G([a,b,c], aa[ii,jj] )][1];
od:
od:
return m;
end proc:

Gen := proc(m,aa,kk2,ii2,jj2)
for ii from 1 to numcol do
for jj from 1 to numcol do
if ii > ii2 and jj > jj2 then
alist := [G([a,b,c], aa[ii,jj] )];
for kk from 1 to numcol do
if kk > kk2 then
m := delta(m,alist,ii,jj,kk);
print(m);
Gen(m,aa,kk+1,ii,jj);
end if:
od:
end if:
od:
od:
end proc:
inited := initit(m,aa);
kk := 1;
Gen(inited,aa,kk,0,0);

Error, (in Gen) invalid left hand side in assignment

@acer 

i feel very expert maple code though i do not understand

@Markiyan Hirnyk 

Hi , Markiyan Hirnyk, i have a long time no question to ask here,

what detail do you want to know to discuss with you?

can you send private email with maple prime to me to discuss, i have a dream to do with maple

@mattcanderson1 

with(StringTools);
with(FileTools);
with(Matlab);
openlink();
evalM(sprintf("a := [1, 2, 3, 4, 5];b := [0, 3, 4, 2, 1]; plot(a, b);"));
closelink();

in fact i mean draw matlab's graph in maple

@tomleslie 

i want to read data in maple then pass to matlab to calculate k means fom maple and then get the result from matlab and then plot in maple to use circle to classify data in maple, add one more field in data to indicate which group it belong to