vv

12453 Reputation

19 Badges

10 years, 0 days

MaplePrimes Activity


These are replies submitted by vv

Please try to formulate mathematically the problem (forget about Maple for the moment).
For example:

Find the C^1 functions f(t,x,p) defined on R^3 such that

   f(t,x(t),x'(t)) = x'(t) * (∂f/∂p)(t,x(t),x'(t))

for any C^1 function x(t) defined on R.

 

For this problem the solution would be

f(t,x,p) = p * C(t,x),  where C is an arbitrary C^1 function on R^2.

Now, try to formulate in this manner your problem.

@asa12 

You must define/explain it. Use quantifiers.

@mikemeson 

In the above example the Array was obtained of course by a simple procedure, but you may define it as complicated as you want and then use surfdata. I don't understand your objection.

@tomleslie 

I posted once a modification like yours including the modified proc() and the post was deleded because of "Copyright" problems.
I did not understand why it was considered so, but this was the fact. A patch like mine (even almost  identical to yours) seems to be acceptable.

@miguelbravo 

Rouben's method is fine for parallel projections of plot curves and plot3d surfaces. It will not work for your dodecahedron.
Using images here it is:

@miguelbravo 

Have you read the last line of my answer? Probably not.

@Adam Ledger 

33..127  are the ASCII codes of "regular" characters.

Note that "!" (ASCII 33) also corresponds to an operator (factorial) but it is posfix and has a special treatment

( op(0, x!)  returns factorial).

@Seb1123 

The ode is solved symbolically; using the initial conditions y(0)=0, D(y)(0)=1;
the solution (Y above) depends on two constants _C3, _C4.

Computing Y' (Y1 above) it is easy to see (even without Maple) that the limit at infinity
exists only when _C3 = _C4 = 0.

I would rather suggest something like

deg:=Pi/180:    # constant
45*deg;
          Pi/4
                          
sin(30*deg);
      1/2
                              
arccos(1/2) / deg;
           60

After all the radian is a natural unit while the degree is an arbitrary (but convenient) one, used for historical reasons.

Note that if such option will be implemented, most of the existent programs will fail for a user who sets it to "degree".

So, you think that I had the impression that the very short SQ procedure is a rigorous proof of the celebrated Toeplitz' conjecture :-)

@Markiyan Hirnyk 

@one man 

The worksheet works in Maple 2016/2017.
The attached version should work in Maple 17 too, but I cannot test it.
If you have problems please attach the result.

Tangent_plane_VV1.mw

@optoabhi 

There are a lot of discussions about this on the forum (just search them).

- The Maple engine is text based, so anyway the 2D math is internally converted to 1D.
- The 2D expressions may contain hidden fields, difficult or impossible to detect visually;
   if something goes wrong it will be difficult to debug.
- Most users which write consistent code prefer the Worksheet mode, 1D for input and 2D for output.
  The 2D input is reserved for presentations.
- A beginner may think that the 2D input mode is easier but when starting a more complex codding he/she will find that the old 1D is more reliable.
  Probably for an occasional user the 2D option is OK.
   

@nm 

You can choose symbolic values for x and y, and after that give them special values in order to simplify the result (optional).

Facxy:=proc(F::algebraic, X::name=anything, Y::name=anything)
  F=1/eval(F,[X,Y]).eval(F,Y).eval(F,X), is(F*eval(F,[X,Y])=eval(F,X)*eval(F,Y))
end:

(I have used Rouben's remark to simplify the procedure.)

 Facxy(cos(x + y + 1) + sin(x - 1)*sin(y + 2),   x=a, y=b);

Facxy(1/(x-1)*y,x=a,y=b);

 

 

I'd recommend to try a different programming and presentation style.
For example, for the curve only with animation:

 

NPar :=proc(F::{list,set}, X::{list,set}(`=`), L, N)
# Natural Parametrization; L=length, N=number of points generated
local n:=nops(X), x:=lhs~(X),  d, i,j, J,t,r,s;
if nops(F) <> n-1 then error "The number of functions must be nops(X)-1" fi;
J:=Matrix(n,n-1, (i,j) -> diff(F[j],x[i]));
d:=seq(LinearAlgebra:-Determinant(J[[1..j-1,j+1..n]])*(-1)^j,j=1..n);
d:=subs(x=~x(t), [d]);  
r:=sqrt(add(d^~2));
s:=dsolve({seq( (diff(x[i](t),t) = d[i]/r), i=1..n), X[](0)}
          ,numeric, output=Array([seq(i*L/N,i=0..N)]) );  
s[2,1][..,2..n+1];
end:

###############################

L1 := 1.2:
f1 := (x1-.5)^4+(x2-1)^4+x3^4+2*x1*x2*x3-L1^4:
f2 := (x1-sin((x1^2+x2^2+x3^2)^.5)^2)^2+(x2-sin(x1)^2)^2+(x3-sin(x1)^2)^2-1:

X0:=fsolve({f1,f2,x1-x3});

{x1 = -.4294971683, x2 = -0.7179716957e-1, x3 = -.4294971683}

(1)

K:=NPar([f1,f2], X0, 9.6, 100):

with(plots):

BG:=implicitplot3d(f2, x1 = -1.5 .. 2.5, x2 = -1 .. 3, x3 = -1 .. 2,
           color = green, transparency = .5, numpoints = 3000, style = surface):

animate(pointplot3d, [ 'K'[1..floor(a)], color=red, style=line, thickness=3 ], a=1..101, background=BG, frames=50);

 

 

 

@asa12 

Sorry, I cannot understand the problem. What eigenvector?   sys2 and sys3 are complex-valued so Minimize is out of the question.

First 101 102 103 104 105 106 107 Last Page 103 of 166