vv

13400 Reputation

20 Badges

9 years, 143 days

MaplePrimes Activity


These are answers submitted by vv

Use any of the following:

eval(x,ans[2]);
                  
rhs(ans[2][]);
                  
op([2,1,2],ans);

 

With the default method and precision, Maple approximates wrongly some of the integrals (or does not approximate at all).
To get the correct results, change for example  A[a]:=...

to

A[a] := int(f(x)*g(x,a),x=-infinity..a-1,numeric,digits=20,epsilon=10^(-10));

 

I think that there are good alternatives.

- For collected terms in an expression see ?LargeExpressions

- The results in your worksheet are easier to read if we use:

alias(alpha=RootOf(_Z^2+5));
alias(beta=RootOf(9*_Z^6-48*_Z^5+30*_Z^3+25));

[this way we have full control on abbreviations]

y(z) = F(z) + int(sqrt(z - sigma)*y(sigma), sigma = 0..z);
intsolve(%, y(z), method = Laplace);

ff:=proc(a)
local aa,sa:=convert(a,set),T,u;
aa:=zip(`[]`,a,[$1..nops(a)]);
for u in sa do T[u]:=NULL od;
for u in aa do T[u[1]]:=T[u[1]],u[2] od;
[seq([T[u]],u=sa)];
end:

ff([[x+y+z],[2*x+y+z],[x-y+z],[2*x+y+z]]);
           
            [[3], [1], [2, 4]]

Seems to be a little faster and easier to read.

It's difficult to understand what you are trying to do.

It seems that you want to use the implicit equation of the spiral r = theta  (given in polar coordinates) and to obtain the points of intersection with some line.

1. You forgot  the third equation:  r = sqrt(x^2 + y^2)

2. eq1 is wrong; LHS should be r  (not r^2) . But RHS is wrong too (see below).

The correct implicut equation of the spiral is

(F(x,y) =)     frem(sqrt(x^2 + y^2) - arctan(y,x), 2.*Pi) = 0

But you should not use it; with polar coordinates is much better. Note that e.g. implicitplot  is not able to draw correctly  F(x,y)=0.

If you want to intersect the spiral with a line (say a*x+b*y+c=0), just solve the system

{r = theta,  a*r*cos(theta) + b*r*sin(theta) + c = 0}

It will have an infinity of solutions; you can use solve followed by evalf(allvalues(%))  or RootFinding.

 

 

 

 

 

 

You must define the points first.

with(geometry):
point(F1,0,1), point(F2,4,1);
ellipse(e1,['foci'=[F1, F2], 'MajorAxis' = 8],[x,y]);

draw(e1);


If you can still open other files (.mw ?), it seems that some of your files are corrupted.

Upload one of them here. Or at least inspect its size, complete path, open it with Notepad (if the OS is Windows) and copy here the first few lines.

Logic:-Satisfy is new in Maple 2016. The variable p was forgotten from the local sequence.

As a workaround you can start with
local p
in your worksheet.

restart;
with(Logic):
local p;
p:=5;
Satisfy(x or y); p;
                     {x = true, y = false}
                               5


Edit. Now p is local to your program and the variable altered by Satisfy is :-p

 

 

To inspect the coordinates use plottools[getdata]

A:=plots:-implicitplot(x^2+y^2=1, x=-1..1, y=-1..1):
plottools:-getdata(A)[3];  # double click the summarized Matrix

convert(%,list,nested);    # if you really want a list

 

simplify(eval(x-y+4*z, [x = (u-v)^2, y =  u^2-3*v^2, z = (1/2)*v*(u-2*v)]));
      0

So, the surface is contained in the plane x-y+4*z=0.

(Of course x-y+4*z can be obtained via Groebner.)

r:=[(u-v)^2, u^2-3*v^2, (1/2)*v*(u-2*v)]:
with(LinearAlgebra):
E:=Vector(diff(r,u))^+ . Vector(diff(r,u)):
F:=Vector(diff(r,u))^+ . Vector(diff(r,v)):
G:=Vector(diff(r,v))^+ . Vector(diff(r,v)):

A := Matrix([[-diff(E,v,v)/2+diff(F,u,v)-diff(G,u,u)/2, diff(E,u)/2, diff(F,u)-diff(E,v)/2],
                     [diff(F,v)-diff(G,u)/2, E, F],
                     [diff(G,v)/2, F, G] ]):
B := Matrix([[0, diff(E,v)/2, diff(G,u)/2], [diff(E,v)/2, E, F], [diff(G,u)/2, F, G]]):
K := simplify((Determinant(A)-Determinant(B))/(E*G-F^2)^2);

                               0

You can use maxima. It is a free, complete and nice CAS which was ported to Android. It has many "common points" with Maple.
Note that on a phone or even on a tablet it is not comfortable to type longer code.

Mathematically this question is related to https://en.wikipedia.org/wiki/Kolmogorov%E2%80%93Arnold_representation_theorem

In your particular case, F(x,y) can be expressed as f(x/y) iff F is homogeneous of degree 0, which is not the case for F(x,y) = x+y-1.
Generally, F must have special properties in order to be expressed via functions of a single variable; e.g. if F has radial symmetry (i.e. F(u)=F(v) for ||u||=||v||) then F(u) can be expressed as f(r), where r = ||u||.

Your function F(n,N) = n+N-1 can be trivially expressed as z-1  for z=n+N, but obviously this does not help.

    floor((n+1)/2)

First 99 100 101 102 103 104 105 Last Page 101 of 117