Adrià Vilanova Martínez | ce8597a | 2021-09-21 18:26:59 +0200 | [diff] [blame] | 1 | % Problema 1.21.b) |
| 2 | % Geometria Diferencial, FME, UPC |
| 3 | % 2021/03/05 |
| 4 | |
| 5 | % > 1.21.b) Trobeu numèricament amb Octave, i representeu per s \in [-1, 1], |
| 6 | % > una corba beta(s) tal que la seva curvatura sigui 0.5, la seva torsió sigui |
| 7 | % > tau(s) = 1 + s, beta(0) = O i el seu triedre de Frenet per s = 0 sigui la |
| 8 | % > base canònica. |
| 9 | |
| 10 | s0 = -1; |
| 11 | s1 = 1; |
| 12 | P = [0; 0; 0]; |
| 13 | Iden = eye(3); |
| 14 | Y0 = [P; Iden(:)]; |
| 15 | |
| 16 | % Integració del sistema dinàmic |
| 17 | % @TODO: Afegir definicions de funcio_corba i les funcions dependents. |
| 18 | [s, Y] = ode45(@funcio_corba, [s0, s1], Y0); |