Question

Me pueden ayudar con este programa de Python
Estoy haciendo un programa que resuelva sistemas de ecuaciones no lineales con el método Newton-Raphson me marca un error pero no se cual es, si me ayudan se los agradecería.

from numpy import *
from numpy.linalg import inv
x=array([0,0,0])

def f(x):
f=zeros(len(x))
f[1]= 3*x[0]-cos(x[1]*x[2])-5.0
f[2]= 9*x[1]+sqrt((x[0]**2)+sin(x[2])+1.06)+0.9
f[3]= exp(-x[0]*x[1])+20*x[2]+9.0
return (f)
def j(x):
j=zeros((len(x),len(x)))
j[0][0]=3.0
j[0][1]=x[2]*sin(x[1]*x[2])
j[0][2]=x[1]*sin(x[1]*x[2])
j[1][0]=x[0]/sqrt(x[0]**2+sin(x[2])+1.06)
j[1][1]=9.0
j[1][2]=cos(x[2])/(2*sqrt(x[0]**2+sin(x[2])+1.06))
j[2][0]=-x[1]*exp(-x[0]*x[1])
j[2][1]=-x[0]*exp(-x[0]*x[1])
j[2][2]=20.0
return(j)
N=50
i=1
jac=0
dx=0
r=0

while 1 jac=j (x)
inversa=inv(jac)
dx=dot(inversa,f(x))
x=x-dx
r=sqrt(dot(f(x),f(x)))
print(r)
i=i+1

Answer 1

Respuesta:

hola esta hecho en python2? o 3 yo creo que es porque las variables estan declaradas mal

eres cientifico de datos?