I had a new homework in maple but there is two questions that i have a hard time figuring them out. If somewone could help me it would be nice.
1. It says that E=L{v1,v2,v3,v4} (all vectors) where v1= (-5,3,2,0,1), v2=(0,-1,1,3,0), v3=(3,1,0,0,2) and v4=(1,4,4,2,0). Discribe E, which is a subset of R5
(and the teacher would like a description like this exemple: E={(x,y) ∈R2 ∣ 2x-3y=2}
2. A matrix: A := Matrix(5, 5, {(1, 1) = 1, (1, 2) = 0, (1, 3) = 0, (1, 4) = 0, (1, 5) = 0, (2, 1) = 0, (2, 2) = 4, (2, 3) = 0, (2, 4) = 0, (2, 5) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = 52, (3, 4) = 0, (3, 5) = 0, (4, 1) = 0, (4, 2) = 0, (4, 3) = 0, (4, 4) = 7, (4, 5) = 0, (5, 1) = 0, (5, 2) = 0, (5, 3) = 0, (5, 4) = 0, (5, 5) = 3})

how do i expand [ln(1+x)]^2??? and find partial derivatives of functions

I'm generating fairly complex shapes and dumping the x,y,z coordinates of each point to a text file. The shapes have a high degree of spherical symmetry, so they usually look like a torus or ellipsoid.
I've been using pointplot3d to plot the shapes using 'style=line'. The result is something that looks close enough to a wireframe plot.
However, for my final thesis I will be submitting, I'd like to show the shapes opaquely (with hidden lines). Basically, the type of plot generated by:
`plot3d(sin(phi)^2, theta=0..2*Pi, phi=0..Pi,coords=spherical,style=wireframeopaque,color=black);`

Hello everybody,
does anyone know how to integrate the output of BsplineCurve in maple 8? I can plot it over a given interval, but if I try to perform the corresponding definite integral, I get the following error message:
Error, (in int) wrong number (or type) of arguments
as opposed to the output of the Spline routine, which can be integrated without problems.
Thank you for your help.

Hi everybody,
is there any Maple package to deal with algebraic curves? The problem I have to solve is the following: my thesis advisor gave me a polynomial equation in two variables, which I know it is a Riemann surface of genus 0. As far as I know, there is a theorem which says that this surface is birationally equivalent to an elliptic curve. There should be also an algorithm which receives the equation as input and outputs the rational parametrization of the Riemann surface. Moreover, this algorithm should be available in a Maple package (my advisor told me this... the problem is that he does not remember the name of the package, nor the function!!!).

I have a small Maple library (created with march and savelib) that I want to distribute to some people.
It seems to work fine, but if I zip it (with the idea of mailing it), the resulting zipped directory under windows seems to throw up security warnings and gets rejected by my outgoing mail client for being "oversize" although it is in fact small.
Are Maple libraries so created just standard files or are they special in some way?
I'd be grateful for any insight.
Toby

Hello !
I need to create a graph of f(a,b,c) = 1+b/a*c,
where "b" goes from 30 to 120,
"c" from 0 and 30,
and "a" would be from "b" to 220. (but if that part is too complicated, "a" going from 120 to 220 would be okay)
The way I visualize this graph is like a typical 3d graph, with no animation, but with a special color mapping.
Let's say my variable "a" would be on the x (horizontal) axis
variable "b" on the y (depth) axis
variable "c" on the z (height) axis
then color of each point would be mapped to the result, f(a,b,c)
How can this be done ?
Thanks you very much !

I need help...
I'm running the code attached below which has procedure sdef(). After obtaining results from it, I do a couple of substitutions Esc:=subs(b=b(v),Esc); Esc:=subs(v=w,Esc); at the very end -- but it turns out that NOT ALL b is substituted for b(v) and not all v is substituted for w in Esc. Been trying to understand what's going on for a day and could not...
The code:
> restart;
> interface(showassumed=0);
>
> sdef:=proc(x)
> global b,v;
> local m, S, S1, S2, EStmp1, EStmp2, betasol, beta;
> m:=1;
>
> S2:=t*(x+v)*beta; assume(beta>=0,beta

Hello,
perhaps my question is really stupid but I cannot find any hints on the help pages (Maple 10):
How can I find the lower and upper bound of a range, i.e, "a" and "b" from "a..b"?
The following does work
lower:=proc(r::range) min($r) end proc ;
but is, of course, very inefficient for large ranges and only works for integer bounds.

Considering updating my old package

COSVAM, which treats complex-octonionic-valued scalars, vectors, and matrices, first of all I want a much better type check of whether an expression is indeed a complex-octonionic scalar.
A complex-octonionic-valued scalar is a linear combination over the complex numbers of the eight units of the octonions 1,e1,...,e7, say. So all the following expressions are valid examples of complex-octonionic-valued scalars

expr1 := 10;
expr2 := e1;
expr3 := e1*I;
expr4 := e1*cos(t)+e2*sin(t);

I have one awful, long way to solve the problem below, but I am sure there is also an elegant solution; unfortunately, it currently escapes me.
Consider a rigid sphere of radisu 1 and center at 0 that rotates about its center. The angular velocity is omega(t) = Vector([ cos(t), sin(t), sqrt(3) ]). Does the path of the point starting at (0,0,1) ever reach the same point at a later time?

A quick question to help with some Classical Mechanics homework.
If the determinant of a matrix is zero does that mean that it cannot be diagonalized ?
Please feel free to elaborate. I truly appreciate everyone's help on this site.
In case your interested the text is "Classical Mechanics" 3rd ed. by Goldstein, Poole, and Safko.
Chapter 4, specifically problem 4.8.
Thanks,
v/r,

is there a way to get maple to perform simple logarithm manipulations i.e.
ln(a*b^c)=ln(a)+c*ln(b) ? I've noodled around with the convert utility with no luck.
I am trying to do a power law regression via the normal equation "showing my work". I have the answer using CurveFitting so I'm just trying to fill in the steps thanks

From an exercise in "Game Physics", we have the following piecewise position vector r. We wish to show that it is well-defined (as well as twice differentiable) at 0. The goal of the exercise is then to show that the Normal vector is not even continuous. So we try:
```
p1,p2 := <t,t^3,0>, <t,0,t^3>;
r := piecewise( t<0, p1, p2);
limit(r,t=0,left);
```

and unfortunately that limit gives
RTABLE(149559852,MATRIX([[t], [t^3], [0]]),Vector[column])
as a result!
The obvious thing to try next is to "push in" the piecewise into the components. But the

Hi,
I have a small problem. I want to export code to Matlab. The code is mostly a long expression with some named constants. The constants are set to some values in Maple and in Maple only the x-variables are therefore unknown in the expression. However, in the exported code, some of the constants are not replaced by their values. Does someone have a clue what this depends on?
Best regards
Johan