acer

32490 Reputation

29 Badges

20 years, 8 days
Ontario, Canada

Social Networks and Content at Maplesoft.com

MaplePrimes Activity


These are replies submitted by acer

@Ramakrishnan Sorry, but I don't understand what you're trying to say.

@Ramakrishnan When I execute your worksheet in my Maple 2018 it runs fine, and t1 gets assigned a numeric value in the first line, and it all seems in order. Please try again. Maybe you're doing something weird.

@Ramakrishnan As I mentioned, the nature of the animation wasn't you main difficulty or issue. You wanted to know how to make it play, programmatically. So I added the SetProperty call and added it to your code (and for which you seem to have thanked Tom).

@wlferguson19 I don't understand what you mean. The pair [x,y] = [sin(3t),cos(3t)] parametrizes a circle, which doesn't intersect itself (leaving alone the originating point, or points following overlap, from discussion).

@wlferguson19 I don't understand what you mean. The period of the path is Pi. It repeats every Pi unit of the parameter. If you make the parameter go from 0 to 2*Pi then the x- and y-coordinates of the intersection points will just be repeated (as will the curve's path).

restart;

eq:=[sin(6*t),cos(2*t)]:

All:={eliminate(Equate(eq,eval(eq,t=s)),[t,s])}:
Temp:=map[2](op,1,All):
Temp:=remove(u->nops(eval({s,t},u))=1, Temp):

LK:=[seq(eval(eq,uu)=uu, uu=Temp)]:
T:=table(LK):
Ans:=[seq(T[uu],uu={seq(lhs(u),u=LK)})];

[{s = Pi, t = 0}, {s = (2/3)*Pi, t = (1/3)*Pi}, {s = (5/6)*Pi, t = (1/6)*Pi}]

BG:=plots:-display(
  plot(map(a->eval(eq,a),Ans),style=point,color=blue,
       symbolsize=25,symbol=solidcircle),
  plot([op(eq),t=0..Pi],color=red),
  gridlines=false
):

plots[animate]( plot, [[[sin(6*A),cos(2*A)]],
                       style=point,symbolsize=25,
                       symbol=solidcircle,color="Green",
                       gridlines=false],
                A=0..2*Pi, frames=200, background=BG );

 

Download intersect_2Pi2.mw

@nepomukk The sin command uses the result of sign on the argument, but it doesn't offer anything but the default lexicographic ordering of variables to determine that sign.

It is possible to force the issue, with a forced sort and a custom call to sign.

restart;

Chg:=module()
  export F;
  local ModuleApply;
  ModuleApply:=proc(ee,fcn::identical(sin,tan,csc,cot),
                    nm::name)
                 subsindets(ee,'specfunc'(fcn),F,nm);
  end proc;
  F:=proc(u::specfunc(anything,{sin,tan,csc,cot}),nm::name)
       local new,o,r,s,t,dum;
       if op(u)::polynom(anything,nm) then
         try
           s:=sign(op(u),[nm,op(indets([op(u)],
                                       And(name,Non(constant)))
                                minus {nm})]);
           if op(u)::`*` then
             (t,o):=selectremove(type,[op(op(u))],
                                 And(polynom(':-anything',nm),
                                     satisfies(p->degree(p,nm)>0)));
             # Alternatively, sort every element of t.
             new:=`*`(sort(map(`*`,op(1,t)/s),'order'=':-plex'(x)),
                      op(2..,t),op(o));
           else
             new:=sort(op(u)/s,'order'=':-plex'(x));
           end if;
           assign(subs(dum=op(0,u),'dum'(new)),
                  subs(dum=op(0,u),'dum'(new)));
           assign(subs(dum=op(0,u),'dum'(op(u))),
                  s*subs(dum=op(0,u),'dum'(new)));
           return s*op(0,u)(new);
         catch:
           return u;
         end try;
       else return u;
       end if;
  end proc:
end module:

 

# Or some larger expression.
#
foo := sin(k1*(c*t-x));

sin(k1*(c*t-x))

Chg(foo, sin, x);

-sin((x-c*t)*k1)

sin(k1*(c*t-x));

-sin((x-c*t)*k1)

-sin((x-c*t)*k1);

-sin((x-c*t)*k1)

# Check it can be called twice, without error.
#
Chg(foo, sin, x);

-sin((x-c*t)*k1)

 

Download signsin.mw

@minhthien2016 

ee := int(tan(x)^9, x);

(1/8)*tan(x)^8-(1/6)*tan(x)^6+(1/4)*tan(x)^4-(1/2)*tan(x)^2+(1/2)*ln(1+tan(x)^2)

convert(simplify(combine(subsindets(ee,
                                    specfunc(ln),
                                    simplify),
                         ln),power),abs)
assuming cos(x)::real;

(1/8)*tan(x)^8-(1/6)*tan(x)^6+(1/4)*tan(x)^4-(1/2)*tan(x)^2-ln(abs(cos(x)))

 

Download intsimp.mw

@awass This might be a clue.

When editing normally within an Execution Group, pressing the <Enter> key executes the statements and pressing <Shift><Enter> inserts a new line and moves the cursor down.

Your problematic examples are in Execution Groups that appear in the .mw file as two adjacent separate XML Equation substructures, which I don't recall seeing before to represent a single statement.

Normally, when one pastes in code containing newlines (LF, \n, U+000A) then all is ok. For example, in the last attempt in your attached worksheet there is such a newline after the Q:=proc(x) and before the sin(x)

But I suspect that something else is getting pasted in, after the sin(x), that may be causing interpretation by the GUI as two separate statements. (I don't mean a hidden semi-colon, as that seems to come as a consequence of further conversion to 1D, not the original cause.) And this seems to cause the GUI to send the first and second portion to the kernel separately.

The first portion being like,

Q:=proc(x)
  sin(x)

and the second portion being,

end proc;

There was a post a month (or several) ago, in which a Mac OSX user was having trouble due to an unexpected character in a text file. I am wondering whether your issue might be similar. For example, something to do with a difference between a carriage return and a newline.

I have never seen anyone report this problem before. On the one hand, there is some evidence that were some changes to the 2D parser in Maple 2018, but on the other hand this situation doesn't seem a general problem. As a guess, I suspect that it might be a Mac OSX specific issue, with something (CR, LF, etc) being misinterpreted and causing an unintended statement split.

If my guess is right, then diagnosing the problem might require examining the original plaintext (if you copied from that, with the mouse), or examining the file, or knowing which key strokes were using in writing it.

@Joe Riel Oh, right, sure, I would believe that. I just meant that the whole topic relates to kernel-based parameter processing in general, (including, but not only) what happens in the debugger.

@José Goulart I figured that you were using an older version, and that for some reason you wanted to keep it a secret.

I've now changed the header on the Question, to denote the Produce as Maple 18. (You could do that too, in future.)

@Joe Riel I changed it to an Answer, as I too was preparing to mention _npassed .

@José Goulart Your last shows 2D Input (which just happens to be in red), in an Execution Group. That's not 1D Input. (You would also have to use F5 or the toggle to switch to 1D input mode, if it's not your default.)

The following is 2D Input. Here I'm using the prefix form of the inert `%*` operator, rather than infix form as was previously showed with 1D input.

restart

kernelopts(version)

`Maple 2018.0, X86 64 LINUX, Mar 9 2018, Build ID 1298750`

interface(imaginaryunit = I)

Ke := `%*`(2*I*E/l^3, `<,>`(`<|>`(6, 3*l, -6, 3*l), `<|>`(3*l, 2*l^2, -3*l, l^2), `<|>`(-6, -3*l, 6, -3*l), `<|>`(3*l, l^2, -3*l, 2*l^2)))

`%*`(2*I*E/l^3, Matrix(%id = 18446883719466239510))

value(Ke)

Matrix(%id = 18446883719466230238)

InertForm:-Display(Ke, inert = false)

0, "%1 is not a command in the %2 package", _Hold, Typesetting

 

Download inert_star.mw

@awass I am not sure that Mapleprimes allows a ( .maple file name extension) workbook to be uploaded directly. But .zip files work.

However, I'd be interested in seeing the original .mw , rather than a workbook .maple file which is a container for .mw Worksheets/Documents.

It would also be interesting to see what you get with the Maple 2018.2 update.

It sounds like you are copying plaintext and pasting it in (where the source is a text file?). But, sorry, the exact steps aren't clear to me.

@Kitonum The result from rsolve could be simplified prior to unapplying, and the result from this new operator could also be evala'd (simplified). It seems a bit faster to use add on the binomial expansion, say if m::posint .

If the operator "Tm" is going to be called multiple times with the same arguments then it could also get, say, option remember, system (sure, I realize that's the recursive procedure originally had it for a finesse on that).

It's hard to compete with the performance of the original recursive procedure, though, which is very fast. So this is mostly in fun.

restart;

Tf := unapply('evala'(evala(rsolve({Tm(t,m)=2*t*Tm(t,m-1)-Tm(t,m-2),
                                    Tm(t,0)=1, Tm(t,1)=t}, Tm(t,m)))),
              [t,m]);

proc (t, m) options operator, arrow; evala((1/2)*(t-(t^2-1)^(1/2))^m+(1/2)*(t+(t^2-1)^(1/2))^m) end proc

Tf(3,4);

577

Alt := (t,m::posint) -> add((t^2-1)^k*binomial(m, 2*k)*t^(m-2*k), k = 0 .. m):

Alt(3,4);

577

Tr := unapply('evala'(evala(rsolve({Tm(t,m)=2*t*Tm(t,m-1)-Tm(t,m-2),
                                    Tm(t,0)=1, Tm(t,1)=t}, Tm(t,m)))),
              [t,m], proc_options=[remember, system]):

Tr(3,4);

577

Altr := proc(t,m::posint)
          option remember, system;
          add((t^2-1)^k*binomial(m, 2*k)*t^(m-2*k), k = 0 .. m);
        end proc:

Altr(3,4);

577

 

Download Tm.mw

@awass It is not helpful to not upload a worksheet that exhibits the problem, and to not tell us exactly which Maple version (major and minor number) you are using.

First 233 234 235 236 237 238 239 Last Page 235 of 595