Ayuda
Publicado por Alberto (23 intervenciones) el 16/11/2018 00:48:34
Hola necesito ayuda.
Adentro del archivo rar están los códigos y el proyecto.
Mi problema es con el ejercicio uno, tengo mi función escrita, creo que esta bien, y el código sisecnl para resolverlo
Sisecnl, también usa otros dos códigos que están dentro de la carpeta.
El problema lo resuelvo con newton con diferencia finita y usando la opción de EGPP (eliminación gaussiana con pivoteo parcial )
tomando como x0=linspace(1,1,20)
Sisecnl lo ejetuto de la siguiente manera:
en la ventana de comando pongo [x,itnewton]=sisecnl('sumidero','nameder',x0) y el resudado que obtengo no es correcto. El verdadero resultado está en la imagen, este resultado es para una tolerancia de 10^-2 y converge en 5 interacciones.
Mi problema es que no puedo encontrar el error.
Si alguien es tan amable de mirarlo, se lo agradecería mucho.
- codigos.rar(71,8 KB)
Valora esta pregunta
0