1
« en: Jueves 16 de Abril de 2009, 08:52 »
Hola que tal, gracias por responder.
Mira el problema que presenté aquí ya lo pude solucionar afortunadamente, pero me surgió otro al tratar de realizar un programa que efectuara el método de Newton-Raphson para sistemas de ecuaciones no linealas, el codigo es:
function [Raiz, Iteraciones] = Newton_Rapson_Sistemas(Gx, Gx1, x, y)
It = 1;
It_limit = 30;
Tol = 0.001;
X = sym('x');
Y = sym('y');
Uxy = feval(Gx, X, Y);
Vxy = feval(Gx1, X, Y);
dUdx = diff(Uxy, X);
dUdy = diff(Uxy, Y);
dVdx = diff(Vxy, X);
dVdy = diff(Vxy, Y);
for q = It : It_limit
dU_dx = dUdx(x,y)
end (Este "end" es improvisado para ir probando el programa cada vez que agrego un comando, el programa no está terminado)
Al correrlo me da el siguiente error:
??? Subscript indices must either be real positive integers or logicals.
Error in ==> Newton_Rapson_Sistemas at 14
dUdx(x,y)
El error al parecer es porque los valores que yo ingreso son x=1.5 y=3.5, y me pide que sean enteros, pero lo que yo deseo que suceda es que a la variable dUdx se le sustituyan los valores de "x" e "y", la variable dUdx contiene lo siguiente:
dUdx = 2*x+y
Por lo tanto lo que yo espero es que dUdx = 6.5.
Espero haberme explicado bien al respecto, si necesitan que aclare algo con gusto lo haré.
Saludos y muchas gracias!
Mostrar Mensajes
