2017-04-07 1 views
-2

Donc je travaille sur une calculatrice intégrale en ce moment, et cela fonctionne assez bien, bien que dès que définissant y comme x ** 2 + x * 2 + 2 ça cesse de fonctionner. Qu'est-ce qui incite mon programme à fonctionner avec y = x ** 2 + 2 mais pas avec y = x ** 2 + x * 2 + 2?Calculatrice intégrale ne fonctionne pas lorsque donné formule donnée plus de 1 x dans la formule

import math as math 

x0 = 0                      
x1 = 0                      

def func(y, x):                    
    return eval(y) 

def func2(z, x):                    
    return eval(z) 

def func3(c, d):   
    a = 0.0                     
    for i in range(1, n+1):                 
     x0 = a + (i-1) * dx                  
     Ai = dx * (c + d)/ 2.               
     a = a + Ai                  
    return a 

y = str(raw_input("Function 1: "))              
z = str(raw_input("Function 2: "))           

a = float(input("Left boundary: "))              
b = float(input("Right boundary: "))              
dx = float(input("Trapezoid width: "))              

n = int((b - a)/dx)                

Area2 = func3(func(y, x0), func(y, x1))                    
Area3 = func3(func2(z, x0), func2(z, x1)) 

Area4 = Area2 - Area3 
if Area4 < 0: 
    Area4 = Area3 - Area2 

print "Area = ", Area4                 
+0

Ce n'est pas clair quant à ce dont vous avez besoin, s'il vous plaît réviser. –

+0

Dans quel sens n'est-il pas clair? Je demande à quelqu'un de comprendre pourquoi il me donnera le résultat correct pour y = x \ ** 2 + 2 z = x \ ** 2 mais pas pour y = x \ ** 2 + 2x + 2 z = x \ ** 2 – Meh

+0

@Meh Avez-vous essayé de déboguer votre code? – nemo

Répondre

0
import math as math 

x0 = 0                      
x1 = 0                 

def f(function, x): 
    function = eval(function) 
    return function 

def func3(b, a, func):   
    area = 0.0 
    dx = (b - a)/n             
    for i in range(1, n+1):                 
     x0 = a + (i-1) * dx 
     x1 = a + i*dx 
     Ai = dx * (f(func, x0) + f(func, x1))/ 2.               
     area = area + Ai                  
    return area 

y = str(input("Function 1: "))              
z = str(input("Function 2: "))           

a = float(input("Left boundary: "))              
b = float(input("Right boundary: "))              
n = int(input("Trapezoid width: "))              


Area2 = func3(b, a, y)                    
Area3 = func3(b, a, z) 

Area4 = Area2 - Area3 
if Area4 < 0: 
    Area4 = Area3 - Area2 

print "Area = ", Area4 

Est-ce plus ou moins ce que vous vouliez faire?