// // Funcion que del problema del perro siguiendo al dueno // // * Parametros del problema son: // * Velocidad del perro (Vp) = 10ms // * El dueño se mueve siguiendo una circunferencia de radio 100 metros y centro en (0,0) // * Posicion inicial del dueño (XD, YD) = (0,50) // * Posicion del dueño en el instante t es (XD, YD) = (100*sin(omega*t),100*cos(omega*t) ) // * La velocidad del dueño VD = 5ms // function V = derivada_perro(t,Z) Vp = 5; VD = 2; radio = 15; omega = VD/radio; // Posicion del perro xp = Z(1); yp = Z(2); // Posicion del dueno XD = radio*sin(omega*t); YD = radio*cos(omega*t); Vx = Vp*(XD - xp)/sqrt((XD - xp)^2 + (YD-yp)^2); Vy = Vp*(YD - yp)/sqrt((XD - xp)^2 + (YD-yp)^2); V = [Vx Vy]; endfunction