MM08 Geometri I Ugeseddel 10 - Institut for Matematik og Datalogi ...
MM08 Geometri I Ugeseddel 10 - Institut for Matematik og Datalogi ...
MM08 Geometri I Ugeseddel 10 - Institut for Matematik og Datalogi ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
3. 7.14.<br />
d) Brug Eulers <strong>for</strong>mel <strong>for</strong> at vise at det langs γ gælder<br />
|τ| = √ −K (Beltrami-Ennepers Sætning).<br />
4. En kurve på en flade siges at være en krumningskurve (“line of curvature“), hvis dens tangent<br />
i enhvert punkt er en hovedretning i dette punkt. Antag at S1 <strong>og</strong> S2 skær hinanden langs en<br />
regulær kurve γ, <strong>og</strong> lad θ(p) være vinklen mellem S1 <strong>og</strong> S2 i punktet p ∈ S1 ∩ S2. Antag at γ er<br />
en krumningskurve på S1. Vis at θ er konstant hvis <strong>og</strong> kun hvis γ <strong>og</strong>så er en krumningskurve<br />
på S2 (Joachimsthals Sætning).<br />
5. Brug maplekoden neden til at vise at alle punkter på fladen z = 3 2 xy − 1 2 x3 er hyperbolske.<br />
Maple<br />
Vi kan bruge Maple til at plotte omdrejningsflader med konstant Gausskrumning, <strong>og</strong>så i tilfælderne<br />
hvor vi ikke kan integrere per hånd. Bemærk brug af Int i stedet <strong>for</strong> int. Vi begynder med<br />
with(plots):setoptions3d(scaling=constrained):<br />
Lad os så prøve med Gausskrumning 0:<br />
sigmaZero:=(u,v) -> [(a*u+b)*cos(v),(a*u+b)*sin(v),sqrt(1-a^2)*u];<br />
a:=0.5; b:=0.5; plot3d(sigmaZero(u,v), u=-1..1, v=0..2*Pi);<br />
positiv Gausskrumning:<br />
sigmaPlus:=(u,v)->[a*cos(u)*cos(v),a*cos(u)*sin(v),<br />
Int(sqrt(1-a^2*sin(s)^2),s=0..u)];<br />
a:=0.5; b:=0.5; plot3d(sigmaPlus(u,v),u=-1..1, v=0..2*Pi);<br />
<strong>og</strong> så negativ Gausskrumning:<br />
sigmaMinus:=(u,v)->[(a*exp(u)+b*exp(-u))*cos(v),<br />
(a*exp(u)+b*exp(-u))*sin(v),Int(sqrt(1-(a*exp(s)-b*exp(-s))^2),s=0..u)];<br />
a:=0.5; b:=0.5; plot3d(sigmaMinus(u,v), u=-1..0, v=0..2*Pi);<br />
Prøv <strong>og</strong>så andre parametre, <strong>og</strong> <strong>for</strong>sikre jer om at fladerne har den specificerede Gausskrumning.<br />
Det er <strong>og</strong>så nemt at skrive rutiner <strong>for</strong> at beregne Gauss- <strong>og</strong> middelkrumning <strong>for</strong> en flade. Vi begynder<br />
med at hente de pakke vi har brug <strong>for</strong>:<br />
with(plots): setoptions3d(scaling=constrained): with(LinearAlgebra):<br />
Derefter definerer vi n<strong>og</strong>le operationer:<br />
DotProd:=proc(u::list,v::list)<br />
convert(zip(‘*‘,u,v),‘+‘);<br />
end proc;<br />
CrossProd:=proc(u::list,v::list)<br />
[u[2]*v[3]-u[3]*v[2],<br />
u[3]*v[1]-u[1]*v[3],<br />
2