Tag: NUMERICAL-METHODS
Estoy tratando de escribir un código de MATLAB para el método de Euler hacia adelante, pero no creo que la salida sea completamente correcta. Este es mi código: function [t ,u] = Euler(f,tspan ,u0,n) T = [tspan(1) :n: tspan(2)]; u = zeros(1,n); u(n) = u0; h = (tspan(2) – tspan(1))/n; . . . Read more
Estoy tratando de implementar el método de la secante usando Octave, y en particular hice una función llamada “secant(“función”,”x0″,”x1″,”tolerancia”)”. Usé el bucle while para calcular la raíz de una función y creo que el problema está en el bucle. Pero el resultado que obtengo no es el que espero. La . . . Read more
Tengo dos puntos en 3D, por ejemplo a = (100, 100, 10) y b = (0, 100, 60), y me gustaría ajustar una línea a través de esos puntos. Sé que la ecuación de la línea en 3D puede tener diferentes formas: Forma vectorial: (x,y,z)=(x0,y0,z0)+t(a,b,c) Forma paramétrica: x=x0+ta y=y0+tb z=z0+tc . . . Read more
Estoy intentando encontrar ex sin usar math.h. Mi código da respuestas incorrectas cuando x es mayor o menor que ~±20. Intenté cambiar todos los tipos double por tipos long double, pero dio basura en la entrada. Mi código es: #include <stdio.h> double fabs1(double x) { if(x >= 0){ return x; . . . Read more
Estoy intentando utilizar algunas de las funciones en MATLAB para resolver numéricamente un par de EDOs de segundo orden acopladas del siguiente formulario: \ddot{x} = f(x,y,\dot{x},\dot{y}) \ddot{y} = f(x,y,\dot{x},\dot{y}). Puedo lograr que funcione con solo una EDO de segundo orden, pero el código que estoy intentando no funciona para un . . . Read more