Thomas Dean

312 Reputation

10 Badges

20 years, 18 days

MaplePrimes Activity


These are questions asked by Thomas Dean

I can not find a cumulative sum procedure in maple2016.

I want to compute the cumulative sum of a list or vector.

A:=[1, 2, 3, 4];

cumsum(A); ## or whatever the maple function is... 

[1,3,6,10]

Tom Dean

I have a long list of two element lists, for example,

A := [ [2,3], [4,5] ,[6,7]];

I want apply '/' to the elements of each sublist.

The result will be

R := [2/3, 4/5, 6/7];

I do not need the list, I could use vectors or matrices.

Is it possible to do this other than by iteration?

Tom Dean

When I was editing the head of the question (? instead of .), its body disappeared. Please, insert it again.

Regard,

Markiyan Hirnyk

Maple 15.

I have a set of equations I can solve manually, but, solve fails.

restart;
eq1 := tgtX[1] = 0;
eq2 := tgtY[2] - y[2]    = m[2]*(tgtX[2]-x[2]);
eq3 := tgtY[3] - y[3]    = m[3]*(tgtX[3]-x[3]);
eq4 := tgtY[4] - y[4]    = m[4]*(tgtX[4]-x[4]);
eq5 := tgtY[2] - tgtY[1] = vy*t[2];
eq6 := tgtY[3] - tgtY[1] = vy*t[3];
eq7 := tgtY[4] - tgtY[1] = vy*t[4];
eq8 := tgtX[2]           = vx*t[2];
eq9 := tgtX[3]           = vx*t[3];
eq10:= tgtX[4]           = vx*t[4];
#
# solve the equations
eqs  := {eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9,eq10};

solvx := solve(eq10,vx);
solvy := solve(eq7,vy);
sol1  := subs(vx=solvx,{eq8,eq9});
sol2  := subs(vy=solvy,{eq5,eq6});

soln  := subs(sol1,{eq2,eq3});

soln  := subs(tgtY[2]=solve(sol2[1],tgtY[2]),soln);
soln  := subs(tgtY[3]=solve(sol2[2],tgtY[3]),soln);
soln  := subs(tgtY[1]=solve(soln[1],tgtY[1]),soln[2]);
soln  := solve({eq4,soln},{tgtX[4],tgtY[4]});

# this returns empty solution
solve(eqs,{tgtX[4],tgtY[4]});

Any ideas?

Tom Dean

I have been here before...  My head is full of cotton, as usual.

 

f := proc (t) 2*t^3+9*t^2-60*t+1 end proc;
deq:=diff(f(t),t);
df := t->deq;  ## this is most likely one of my problems.
plot(df(t),t=-10..10);
isneg := x -> if is(df(x) < 0) then df(x) else 0 end if;  ## and, another??
plot(isneg(t),t=-10..10);
Plt1:=plot(isneg(t),t=-10..10);
plots:-display(Plt1);
isneg(0.12017234);

The if statement is not fully evauated.

I am missing something.  What?

Tom Dean

First 12 13 14 15 Page 14 of 15