I have following lists

a,a,a,b,c
b,b,b,c,e
a,c,e,e,e
a,a,a,b,c
b,b,b,c,e
a,c,e,e,e
a,a,a,b,c
b,b,b,c,e
a,c,e,e,e
a,a,a,b,c

I want to write the bare rewrite rules for these so that if any list is given as input these rules produce the next one. It is obvious every third is same as first.

I want to get rid of the unnecessary output lines after the evaluation of the following ocde

restart; with(linalg);
a := 2; b := .29; d := 1.85;
for h from .5 by .1 to 1 do
eq1 := x*(-b*x^2-x+1);
eq2 := y*((a*x*x)/(b*y^2)-d-h*y);
S := solve({eq1, eq2}, {x, y});
SS := solve(subs(S[3], {omega^4+(h*y+x)*omega^2+h^2*x-y}), {omega});
tau := simplify(subs(S[3], subs(SS[3], (b^2*h*y+a*x)/omega)))
end do

Hi all,

Say i have two Vectors:

V1:=Vector(2,[a,b]);

V2:=Vector(4,[d,b,e,a]);

and how do I get

V3:=Vector([d,e]);

?? I currectly use a crude way using 'set' operations:

V3:=convert(V2,set) minus convert(V1,set);

V3:=convert(V3,Vector); # Does not work straightford, so

V3:=convert(V3,list):
V3:=convert(V3,Vector); # Which works,

Also, is there a strightforward way to which elements are the same and their indices?

Let us consider the group of elements {(a,b,c,d): a::integer,b::integer,c::integer,d::integer,a+b+c+d=0, a=c mod 12} with the component-wise addition. How to represent it in Maple? Is it possible? What I know is

restart;

with(GroupTheory):

CyclicGroup(infinity);     

GroupTheory:-CyclicGroup(infinity, form = "fpgroup")

G := DirectProduct(CyclicGroup(infinity), CyclicGroup(infinity), CyclicGroup(infinity), CyclicGroup(infinity));

when solving using output=listprocedure i face this message :

Sys:= {seq}(eq[k], k= 1..18):

fsolve(Sys, {seq}(a[k], k= 1..18),output=listprocedure):

Error, (in fsolve) output = listprocedure is an invalid option

what should i do? how can i list my output to use in latter commands? for example i want to use the output of the problem to plot an expression that the outputs are the coefficients of that expression?

so i got this
> low := 170; high := 195; j := 0; byy := 5; plott := vector((high-low)/byy+1);
with(plots);
for i from low by byy to high do
proc (X) options operator, arrow; X(t) end proc;
proc (Y) options operator, arrow; Y(t) end proc;
ode := diff(X(t), t) = evalf(subs(x = X(t), y = Y(t), subs(vvx = VXX[i-1], vvx))), diff(Y(t), t) = evalf(subs(x = X(t), y = Y(t), subs(vvy = VYY[i-1], vvy)));
ics := X(0) = Cordx[i-1], Y(0) = Cordy[i-1];
ds := dsolve([ode, ics...

I need three orthogonal unit vectors, k, i, and j. I give the first two vectors, k=[sqrt(2)-1,-1,1] and i=[1,sqrt(2)-1,0], then use crossproduct to produce the third vector j=crossprod(k,i) . But it fails to produce the correct vector since dotprod(j,k) is not 0. 

I have a two-variable matrixfunction A(x,y) with entries A_i,j(x,y).
I want to use mtaylor for each entries and the result to store in a matrix B.
Is it possible to do using the map command?
If not, what methods do exist?
(Instead of typing mtaylor for each entires.)

Sandor

restart:

dsys := {diff(x(t),t)=y(t),diff(y(t),t)=-x(t),x(0)=1,y(0)=0}:

dsn := dsolve(dsys,numeric):

dsn(1);          [t = 1., x(t) = HFloat(0.5403023317785669),  y(t) = HFloat(-0.8414711011553073)]

rhs(dsn(1)[2]);                   HFloat(0.5403023317785669)

After solving  this system I wish to obtain a numeric value to the integral:

int(x(t...

what should i do to write matirx of this form in maple  for nth order

[1        4        1                                    ]

[B_-1,0   B_0,0   B_1,0                                   ]

DMS_13.mw

this code executes fine in  M13. but in M15 theres a problem. doesn't like the local.

i wish maple had a compatability/update feature like mathematica...

Hi,

I recently have a problem when I try to use dsolve/numeric to solve an IVP. Suppose we have diff(y,x) = w(y) and y(0) =  y0. We can solve this IVP numerically by dsolve/numeric. In Matlab, we can obtain y and  the derivatives yp using the command [y,yp] = deval(sol,t) at any time t within the range we pre-defined. Basically, matlab ode solvers return a solution  skeleton first, then evaluate y and yp of the other points by a piecewise interpolation which is called...

with insequence=true command, first, the first option is showed, then second, then third and so on.. but i want to show all options from in grpah, then animate in sequence.

hi  projet_TW.mw
 procedure_final_(1.mw

could someone help me.thanks

Hi,

     I'm reading in a file as

eigvect:=ImportMatrix("C:\\Users\\Sorted_Eigenvectors_1.txt", source=csv):

but I want to generalize this to 

eigvect:=ImportMatrix("C:\\Users\\Sorted_Eigenvectors_i.txt", source=csv):

where i is a positive integer, I tried

eigvect:=ImportMatrix("C:\\Users\\Sorted_Eigenvectors_%d.txt",i, source=csv):

but this is incorrect. What is the correct syntax?

Thanks!