integral en el sage

296 days ago by habacuc@tecpabellon

playerito
EJERCICIO I
para obtener el resultado de cuanto camino para adelante y para atras tenemos que separar la primer integral en dos partes y luego los resultados que te den se restan y buelbes a el resultado inicial
#integral original d=(x^(2)-(2*x)-8) area=integral(d,1,6) print float(area) 
       
-3.33333333333
-3.33333333333
EL VALOR ABSOLUTO ES POR QUE NO PUEDE ABER DISTANCIAS NEGATIVAS
#esta integral te da el resultado de todo lo que se movio el cangrejo........ d=abs((x^(2)-(2*x)-8)) total=integral(d,1,6) print float(total) 
       
32.6666718426
32.6666718426
#integral 1..... d=abs((x^(2)-(2*x)-8)) total=integral(d,1,4) print float(total) 
       
18.0
18.0
#integral 2........ #var ('t') d=abs((x^(2)-(2*x)-8)) total=integral(d,4,6) print float(total) 
       
14.6666666667
14.6666666667
EJERCICIO II
Zf= lambda x:x^2+1 a=0 b=2 n=10 Dx=(b-a)/n rsum=sum([f(a+i*Dx)*Dx for i in range(10)]) print float(rsum) 
       
4.28
4.28
plot(f(x),x,a+0,b+2.5)a) y=x², -1,1 
       
f= lambda x:x^2+1 a=0 b=2 n=10 Dx=(b-a)/n rsum=sum([f(a+i*Dx)*Dx for i in range(10)]) print float(rsum) 
       
2627000.0
2627000.0
f= lambda x:x^2+1 a=0 b=2 n=10 Dx=(b-a)/n rsum=sum([f(a+(i+.05)*Dx)*Dx for i in range(10)]) print float(rsum) 
       
4.3162
4.3162
f= lambda x:x^2+1 a=0 b=2 n=10 Dx=(b-a)/n rsum=sum([f(a+(i+1)*Dx)*Dx for i in range(10)]) print float(rsum) 
       
5.08
5.08
EJERCICIO III
a=-1; b=1 f= lambda x:x^2 Lb=[[b,f(b)],[b,0],[a,0],[a,f(a)]] Lf=[[i/10,f(i/10)]for i in xrange(10*a,10*b)] P=polygon(Lb+Lf, rgbcolor=("purple")) Q=plot(f(x),x,a-0.5,b+0.5) show(P+Q) 
       
a=0; b=2 f= lambda x:x^3-5*x^2+1 Lb=[[b,f(b)],[b,0],[a,0],[a,f(a)]] Lf=[[i/10,f(i/10)]for i in xrange(10*a,10*b)] P=polygon(Lb+Lf, rgbcolor=("green")) Q=plot(f(x),x,a-0.5,b+0.5) show(P+Q) 
       
a=0; b=2 f= lambda x:x^2-1 Lb=[[b,f(b)],[b,0],[a,0],[a,f(a)]] Lf=[[i/10,f(i/10)]for i in xrange(10*a,10*b)] P=polygon(Lb+Lf, rgbcolor=("yellow")) Q=plot(f(x),x,a-0.5,b+0.5) show(P+Q) 
       
a=-1; b=2 f=lambda x:-x^3+3*x-1 Lb=[[b,f(b)],[b,0],[a,0],[a,f(a)]] Lf=[[i/10,f(i/10)]for i in xrange(10*a,10*b)] P=polygon(Lb+Lf, rgbcolor=("red")) Q=plot(f(x),x,a-0.5,b+0.5) show(P+Q) 
       
N(integral((x^(2)),(x,-1,1))) 
       
0.666666666666667
0.666666666666667
N(integral((x^(3)-5*(x)^2+1),(x,-0,2))) 
       
-7.33333333333333
-7.33333333333333
N(integral((x^(2)-1),(x,0,2))) 
       
0.666666666666667
0.666666666666667
N(integral((-x^(3)+3*x-1),(x,-1,2))) 
       
-2.25000000000000
-2.25000000000000