[ Foro de Python ]

ProgramacionSimbolica Y Analitica

08-Jun-2018 18:41
Invitado (Rodrigo)
1 Respuestas

Tengo Este codigo y necesito este while " El que esta comentado al final" que  tengo que convertir de lenguaje analitico a lenguaje simbolico en python

Podrian darme  consejos ?

 
# -*- coding: utf-8 -*-
 
from mpmath import *
 
def Bisect(Func, a, b):
	eps = 1e-10
	itmax = 100
 
	x = a; fa = Func(x)
	if (fabs(fa) == 0e0): return (x,0)
	x = b; fb = Func(x)
	if (fabs(fb) == 0e0): return (x,0)
	if (fa*fb > 0): return (x,1)
 
	for it in range (1,itmax+1):
		x = 0.5e0 * (a + b)
		fx = Func(x)
		if (fa*fx > 0): a = x
		else: b = x
		if (((b-a) <=eps*fabs(x)) or (fabs(fx) <= eps)): return (x,0)
	print("Bisect: max. no. of iterations exceeded !"); return (x,2)
 
def main():
	f=lambda x: (x**4-5*x**2+4)*sin(2*x)
	y=findroot(f, [-3.5,3.5], solver="bisect", maxsteps=100)
	print(y)
	print(Bisect(f,1,2))
	xmin = -3.5
	xmax = 3.5
	h= 0.1
 
	#a = xmin
	#while (a < xmax):
		#b = a+h
		#(x,ierr) = Bisect(func,a,b)
		#if((ierr == 0) and (x !=b)):
			#print (x = {0:8.5f} in ({1:6.2f},{2:6.2f}) f(x) = {3:7.0}.
		#format (x,a,b,func(x)))
			#a=b
plot([f], [-3.5,3.5])
 
 
if __name__== "__main__":
	main()
 


 


17-Jun-2018 00:32
Nacho Cabanes (+36)

No entiendo. Todo suena "bastante a Python". ¿Has probado simplemente a quitar los comentarios y ver si alguna línea concreta no se interpretar correctamente?






Si ya eres usuario del sistema, puedes contestar desde tu cuenta y así ganar prestigio.

Si sólo eres un visitante, puedes optar por...