blob: 1ecd3c1976cd59fe1bc23a436db7e2a3c823eadc [file] [log] [blame]
% Problema 1.21.b)
% Geometria Diferencial, FME, UPC
% 2021/03/05
% > 1.21.b) Trobeu numèricament amb Octave, i representeu per s \in [-1, 1],
% > una corba beta(s) tal que la seva curvatura sigui 0.5, la seva torsió sigui
% > tau(s) = 1 + s, beta(0) = O i el seu triedre de Frenet per s = 0 sigui la
% > base canònica.
s0 = -1;
s1 = 1;
P = [0; 0; 0];
Iden = eye(3);
Y0 = [P; Iden(:)];
% Integració del sistema dinàmic
% @TODO: Afegir definicions de funcio_corba i les funcions dependents.
[s, Y] = ode45(@funcio_corba, [s0, s1], Y0);