Here are 14 problems to challenge your computer algebra system.
If you can solve any of the following then post your solution with
the following extra information if known.
-
Given that the arc lengths of the 3 principal perimeters of an ellipsoid
are 0.989, 0.853 and 0.712 meters, find the surface area in square meters
of the ellipsoid to 33 significant digits.
Note: A principal perimeter is the curve that is the intersection of the
ellipsoid and a plane that contains 2 of the 3 principal semi-axis.
-
Find a conic that touches the following 5 given conics?
11689*x^2-9600*x*y+2169*y^2-2204198*x+908304*y+103910536=0,
11281*x^2-10800*x*y+4036*y^2-1981082*x+1110832*y+91469653=0,
221*x^2-264*x*y+144*y^2+11368*x-1056*y+269456=0,
17425*x^2-7560*x*y+17056*y^2-1635370*x-554456*y+50763289=0,
58825*x^2-6160*x*y+41764*y^2+4470090*x-1814952*y+98454609=0 .
If you can't find the coefficients exactly then can you find a numerical
approximation to them to 17 significant digits.
Note: By touching , I mean that the pair of conics have a common tangent
line at an intersection point.
-
What is the shortest path along the surface of the ellipsoid with equation
x^2/17^2 + y^2/8^2 + z^2/3^2 = 1 ,
starting at point 1. (x=7, y=5, z>0) and going to point 2. (x=-9, y=-3, z<0).
Also find the arc length of the shortest path. Find the exact answer in terms
of known special functions (if possible) and also find numerical
approximations to 17 significant digits.
-
What is the shortest path along the surface of the torus with equation
(x^2+y^2+21+z^2)^2-100*x^2-100*y^2 = 0 ,
starting at point 1. with coordinates, (x=5,y=0,z=2) and finishing at
point 2. with coordinates, (-4,2,(-41+20*5^(1/2))^(1/2)) .
Also find the arc length of the shortest path? Find the exact answer,
in terms of known special functions (if possible) and also find numerical
approximations to 17 significant digits.
-
Optimize the Campbell, Baker , Hausdorff formula to order 12.
The formula is :
ln(exp(t*A) @ exp(t*B)) = (A+B)*t + 1/2*c(A,B)*t^2
+ (1/12*c(A,c(A,B)) - 1/12*c(B,c(A,B)))*t^3 + O(t^4), (to order 3)
where t is a real parameter and A and B are square matrices, and "c" is
the commutator defined by: c(A,B) = A @ B - B @ A, and @ denotes matrix
multiplication. It is known that all higher order terms
(higher than first order in "t") can be expressed in
terms of linear combinations of nested commutators. By optimize I mean
minimize the number of commutators. For example , for the above given
CBH series to third order can be optimized as follows:
s[1] = A+B;
s[2] = 1/2*c(A,B);
s[3] = 1/6*c(A-B,s[2]);
CBH = (((s[3]*t+s[2])*t)+s[1])*t;
Cost is only 2 commutators plus a few other operations.
-
Does this equation x^5 - 5*x^3 + 5*x -5 = 0, have a solution in terms
of radicals. If it does then find the radical expression?
-
Construct the regular polygon with 2^32-1 sides.
If this is too difficult then try the regular 17'gon and then the
regular 257'gon and then the regular 2^16+1 gon.
By construct in the weak sense I mean express the coordinates of the the
polygon in terms of square roots and the other 4 basic arithmetic operations
( add, subtract, multiply and divide).
In the strong sense , I mean write out a list of ruler and compass moves
that any good "student" could follow. You are given 2 points on a sheet of
paper and that is your unit distance and you must construct the n'gon
in a circle of unit radius, where one of your given points is the center of
the unit circle and the other one must be a vertex of the regular n'gon.
A possible example of a ruler and compass move:
step number 1897. Draw a circle with center the point given in step
number 643 and radius given in step number 783 . This intersects the
line given in step number 351 in 2 points. The intersection point
closest to the point from step number 498 , will be called the point
determined by this step number 1897.
Try to minimize the number of steps.
-
Is x rational, where x is the real root of F(-.2,.2;.5;x)=5/6 ,
and 0<=x<=1 ?
Note: F denotes the Gauss hypergeometric function and F(-.2,.2;.5;x)
means that the numerator parameters are -.2 and .2 , the denominator
parameter is .5 and the argument is x .
-
Plateau problem. Find the surface with minimal surface area and
boundary given by the following space curve.
x(t) = 23*cos(t)+4*sin(2*t)+cos(4*t),
y(t) = 11*sin(t)-3*cos(2*t)-sin(3*t),
z(t) = 5*sin(t) -2*sin(3*t)+cos(5*t)
for 0<=t<=2*Pi
Find the surface and its surface area to 17 significant digits.
-
Find the maximum of J(10000,x) for 0<=x<=infinity and where it
occurs , both to 17 significant digits.
Note: J denotes the Bessel function of the first kind , and J(10000,x)
means that the order is 10000 and the argument is x .
-
Let u(x,y) be the solution to the following second order linear
partial differential equation, (Poisson's equation):
diff(u(x,y),x,x) + diff(u(x,y),y,y) = -2 ,
on the interier of the L-shaped region with the following boundary
conditions, it is zero everywhere on the following L-shaped boundary:
join (0,0) to (3,0) to (3,1) to (1,1) to (1,2) to (0,2) to (0,0) in all
cases with straight line segments.
What is the maximum value of u(x,y) and where does it occur?
Find numerical approximations to 33 significant digits.
Note: diff(u(x,y),x,x) denotes the second partial derivative of
u(x,y) with respect to x twice.
-
What is the hyper-surface area of the following 4-dimensional
hyper-ellipsoid:
w^2/3^2 + x^2/7^2 + y^2/11^2 + z^2/19^2 = 1 ,
to 17 significant digits?
Can the triple integral be reduced to a double integral?
a single integral? a known special function?
-
What is an efficient method to evaluate a high degree Bernstein
polynomial at a given point, given the Bernstein polynomial in the
form of a array of function values.
Example: Given the degree n=20000,
and an array "a" indexed from 0 to n , such that
a[i] = f(i/n) for i from 0 to n, to a precision of 16 significant digits.
We are also given the function f(x) = cos(x) .
The point "v" is .4567890123456789 .
Evaluate the Bernstein polynomial for f(x)=cos(x) of degree n=20000
at the point "v". We would like only 10 significant digits for the
evaluation. Can this be done without forming the Bernstein polynomial?
Can your computer algebra system write a program to do this that
could be executed in some other language like C or fortran, that
could run without causing exponent overflow or exponent underflow error
messages, and that only needs 16 digits for all intermediate operations?
Note: the Bernstein polynomial of degree n for function f(x) is defined by:
b(x) = sum( f(i/n)*binomial(n,i)*x^i*(1-x)^(n-i), for i from 0 to n).
-
Find the solution of the following second order linear differential
equation.
(9*x+1)^5*((x-1)^2+1/37)^4*diff(y(x),x$2) + (3*x+1)*diff(y(x),x)
+ (x^5-3*x^2+8)*y(x) = 0,
with the following initial conditions:
y(0)=1, D(y)(0)=0 .
Find y(10) and D(y)(10) to 33 significant digits.
Note: diff(y(x),x$2) denotes the second derivative of y(x)
with respect to x, and D(y)(0) denotes the derivative of y(x)
with respect to x, evaluated at x=0 .