0% found this document useful (0 votes)
12 views19 pages

Tarea4 Matlab

This document contains calculations of limits, derivatives and Taylor polynomials using MATLAB. Several functions are analyzed around different points and their Taylor polynomials are graphed alongside the original functions to show the approximation improves as the degree increases.
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
12 views19 pages

Tarea4 Matlab

This document contains calculations of limits, derivatives and Taylor polynomials using MATLAB. Several functions are analyzed around different points and their Taylor polynomials are graphed alongside the original functions to show the approximation improves as the degree increases.
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 19

Tamara Gascón Moreno

Calculo

Ejercicio 1)
Calcula con MatLab los límites siguientes, tomados en la Tarea 3:
1
a) lim ( ) ln⁡(𝑥)
𝑥→∞ 𝑥

>> syms x

>> f1=(tan(1/x)*log(x))

f1 =

tan(1/x)*log(x)

>> limit(f1,x,Inf)

ans =

b) lim 𝑠𝑒𝑛(𝑥)ln⁡(1 − 𝑒 −𝑥 )
𝑥→∞

>> syms x

>> f2=(sin(x)*log(1-exp(- x)))

f2 =

log(1 - exp(-x))*sin(x)

>> limit(f2,x,Inf)

ans =

𝑠ℎ(𝑥 4 )
c) lim 2𝑥 4 +5𝑥7
𝑥→0

>> syms x

>> f3=(sinh(x^4)/(2*x^4+5*x^7))

f3 =

sinh(x^4)/(5*x^7 + 2*x^4)

>> limit(f3,x,0,'right')

ans =

1/2
Tamara Gascón Moreno
Calculo

1
d) lim+ 𝑥 ln(𝑥)+1
𝑥→0

>> syms x

>> f4=x^(1/log(x)+x)

f4 =

x^(x + 1/log(x))

>> f4=log(x^(1/log(x)+1))

f4 =

log(x^(1/log(x) + 1))

>> f4=1/(log(x)+1)*log(x)

f4 =

log(x)/(log(x) + 1)

>> exp(limit(f4,x,0,'right'))

ans =

exp(1)

√𝑥 2 −3𝑥+4 4𝑥
e) lim ( 𝑥
)
𝑥→∞

%ln L = 1 -> eln L = e1 -> L=e

>> syms x

>> f5=((sqrt(x^2-3*x+4)/x)-1)*(4*x)

f5 =

4*x*((x^2 - 3*x + 4)^(1/2)/x - 1)

%Donde L=
Tamara Gascón Moreno
Calculo

>> limit(f5,x,Inf)

ans =

-6

1
𝑒 𝑥 +2
f) lim− 𝑥 2 +1
𝑥→0

>> f6=(exp(1/x)+2)/x^2+1

f6 =

(exp(1/x) + 2)/x^2 + 1

>> limit(f6,x,0,'left')

ans =

Inf

El resultado tiene que ser 2. Lo calculamos por partes:

>> Lim1=limit(exp(1/x)+2,x,0,'left')

Lim1 =

>> Lim2= limit(x^2+1,x,0,'left')

Lim2 =

Lim1/Lim2

ans =

√𝑒 𝑥2 −1·𝑡𝑎𝑛3 (3𝑥)·ln⁡(1+2𝑥)
g) lim+ 𝑥
𝑥→0 𝑎𝑠𝑒𝑛( )·(cos(2𝑥 2 )−1
2
√𝑒 𝑥2 −1·𝑡𝑎𝑛3 (3𝑥)·ln⁡(1+2𝑥)
h) lim− 𝑥
𝑥→0 𝑎𝑠𝑒𝑛( )·(cos(2𝑥 2 )−1
2

Aplicamos infinitésimos equivalentes


Tamara Gascón Moreno
Calculo

2
𝑒 𝑥 ~⁡𝑥 2
𝑡𝑎𝑛3 (3𝑥)~(3𝑥)3
𝑙𝑛(1 + 2𝑥) ~2𝑥
𝑥 𝑥
𝑎𝑠𝑒𝑛 ( ) ~ ( )
2 2
2
2)
2𝑥 2
𝑐𝑜𝑠(2𝑥 − 1~ ( )
2

>> syms x

>> f7=(abs(x)*(3*x)^3*(2*x))/((x/2)*(-((2*x^2)^2)/2))

f7 =

-(54*abs(x))/x

%Calculamos por la derecha

>> limit(f7,x,0,'right')

ans =

-54

%Calculamos por la izquierda

>> limit(f7,x,0,'left')

ans =

54

|𝑥|
Calculamos -54 𝑙𝑖𝑚 ( 𝑥 )
𝑥→0

>> limit(abs(x)/x,x,0)

ans =

NaN

No existe

Ejercicio 2)
Sea la función f(x) que te corresponda según el primer dígito comenzando por la derecha de tu
documento de identidad (DNI, tarjeta de residencia o pasaporte).

Usando MatLab:

a) Calcula la derivada de la primera, segunda y tercera de la función f(x).


Tamara Gascón Moreno
Calculo

(2 + 𝑥)𝑒 𝑥+1

>> d1=simplify(diff((2 + x)*(exp(x + 1))))

d1 =

exp(x + 1)*(x + 3)

%(3+x)ex+1
Tamara Gascón Moreno
Calculo

>> d2=simplify(diff(d1))

d2 =

exp(x + 1)*(x + 4)

%(4+x)ex+1

>> d3=simplify(diff(d2))

d3 =

exp(x + 1)*(x + 5)

%(5+x)ex+1

b) Evalúa dichas derivadas en x=0 y x=a (Si no se pudieran evaluar en x=0, calcula su
límite para x  0. Si tampoco existiera éste , calcula el límite para x  0+

>> x=0

x=

>> subs(d1)

ans =

3*exp(1)

>> subs(d2)

ans =

4*exp(1)

>> subs(d3)

ans =

5*exp(1)

c) Calcula los tres primeros términos no nulos del desarrollo de Maclaurin de f(x). Si no
existiera el desarrollo de Maclaurin, indica por qué.

>> syms x

>> f = (2+x)*exp(x+1);

>> ts = taylor(f, x, 0, 'Order', 3)

ts =
Tamara Gascón Moreno
Calculo

2*exp(1)*x^2 + 3*exp(1)*x + 2*exp(1)

d) Calcula el polinomio de Taylor de grados 1,2 y 3 (n= 2,3,y 4) para x ≈a y represéntalos


junto con f(x). (Adapta los rangos de los ejes para que se vea bien la mejora de la
aproximación al acrecentar el grado).

>> syms x

>> f = (2+x)*exp(x+1);

>> ezplot(f,[-5 5])

>> ezplot(f,[0 4])

>> ts = taylor(f, x, -1, 'Order', 4)

ts =

2*x + (3*(x + 1)^2)/2 + (2*(x + 1)^3)/3 + 3

>> ezplot(ts,[-5 5])


Tamara Gascón Moreno
Calculo

Ejercicio 3)
Calcula con MatLab los polinomios de orden 2 de todas las funciones que se indican,
expresándolos de la forma a2x2+a1x+a0 las siguientes funciones (puedes utilizar el comando
simplify y representa cada polinomio junto con la función correspondiente, de modo que
se visualice en qué tramo se aproxima bien y en qué tramo empieza a distanciarse.

0) √1 + 𝑥 − 1⁡𝑝𝑎𝑟𝑎⁡𝑥 ≃ 0

>> syms x

>> f1=(sqrt(1+x)-1)

f1 =

(x + 1)^(1/2) – 1

>> ezplot(f1,[-2 2 -2 2])

>> ts1 = taylor(f1, x,0, 'Order', 3)

ts1 =

- x^2/8 + x/2

>> ezplot(ts1,[-2 2 -2 2])


Tamara Gascón Moreno
Calculo

𝜋
1) 𝑠𝑒𝑛 (𝑥 + ) 𝑝𝑎𝑟𝑎⁡𝑡𝑜𝑑𝑜⁡𝑥 ≃ 0
4

>> syms x

>> f2=sin(x+(pi/4))

f2 =

sin(pi/4 + x)

>> ezplot(f2,[-2 2 -2 2])

>> ezplot(f2,[-10 10 -10 10])


Tamara Gascón Moreno
Calculo

>> ts2 = taylor(f2, x,0, 'Order', 3)

ts2 =

(2^(1/2)*x)/2 + 2^(1/2)/2 - (2^(1/2)*x^2)/4

>> ezplot(ts2,[-2 2 -2 2])

2) √𝑥⁡𝑝𝑎𝑟𝑎⁡𝑥 ≃ 2

>> f3=sqrt(x)

f3 =

x^(1/2)

>> ezplot(f3,[-0.5 2 -0.5 2])


Tamara Gascón Moreno
Calculo

>> ts3 = taylor(f3,x,2, 'Order', 3)

ts3 =

2^(1/2) + (2^(1/2)*(x - 2))/4 - (2^(1/2)*(x - 2)^2)/32

>> ezplot(ts3,[-5 5 -5 5])

3) e𝑥 ⁡𝑝𝑎𝑟𝑎⁡𝑥 ≃ 1

>> f4=exp(x)

f4 =

exp(x)

>> ezplot(f4,[-5 5 -5 5])

>> ts4 = taylor(f4,x,1, 'Order', 3)

ts4 =

exp(1) + exp(1)*(x - 1) + (exp(1)*(x - 1)^2)/2

>> ezplot(ts4,[-5 5 -5 5])


Tamara Gascón Moreno
Calculo

1
4) 𝑥 2⁡
𝑝𝑎𝑟𝑎⁡𝑥 ≃1

>> f5=(1/x^2)

f5 =

1/x^2

>> ezplot(f5,[-5 5 -5 5])

>> ts5 = taylor(f5,x,1, 'Order', 3)

ts5 =

3*(x - 1)^2 - 2*x + 3

>> ezplot(ts5,[-5 5 -5 5])


Tamara Gascón Moreno
Calculo

3
5) √x⁡𝑝𝑎𝑟𝑎⁡𝑥 ≃ 1

>> f6=(x^(1/3))

f6 =

x^(1/3)

>> ezplot(f6,[-2 5 -2 5])

>> ts6 = taylor(f6,x,1, 'Order', 3)

ts6 =

x/3 - (x - 1)^2/9 + 2/3

>> ezplot(ts6,[-5 5 -5 5])


Tamara Gascón Moreno
Calculo

𝜋
6) cos(𝑥) 𝑝𝑎𝑟𝑎⁡𝑥 ≃ − 6

>> f7=(cos(x))

f7 =

cos(x)

>> ezplot(f7,[-5 5 -5 5])

>> ts7 = taylor(f7,x,(-pi/6), 'Order', 3)

ts7 =

pi/12 + x/2 + 3^(1/2)/2 - (3^(1/2)*(pi/6 + x)^2)/4

>> ezplot(ts7,[-5 5 -5 5])


Tamara Gascón Moreno
Calculo

7) 𝑠ℎ(𝑥)𝑝𝑎𝑟𝑎⁡𝑥 ≃ ln(2)

>> f8=(sinh(x))

f8 =

sinh(x)

>> ezplot(f8,[-5 5 -5 5])

>> ts8 = taylor(f8,x,(log(2)), 'Order', 3)

ts8 =

sinh(6243314768165359/9007199254740992) +
(sinh(6243314768165359/9007199254740992)*(x -
6243314768165359/9007199254740992)^2)/2 +
cosh(6243314768165359/9007199254740992)*(x - 6243314768165359/9007199254740992)

>> ezplot(ts8,[-5 5 -5 5])


Tamara Gascón Moreno
Calculo

8) ln⁡(𝑥)⁡𝑝𝑎𝑟𝑎⁡𝑥 ≃ 𝑒

>> f9=(log(x))

f9 =

log(x)

>> ezplot(f9,[-5 5 -5 5])

>> ts9 = taylor(f9,x,exp(1), 'Order', 3)

ts9 =

(1125899906842624*x)/3060513257434037 +
log(3060513257434037/1125899906842624) - (633825300114114700748351602688*(x -
3060513257434037/1125899906842624)^2)/9366741398929500034245406117369 - 1

>> ezplot(ts9,[-5 5 -5 5])


Tamara Gascón Moreno
Calculo

9) 𝑥⁡𝑒 −𝑥 ⁡𝑝𝑎𝑟𝑎⁡𝑥 ≃ 1

>> f10=(x*exp(-x))

f10 =

x*exp(-x)

>> ezplot(f10,[-5 5 -5 5])

>> ts10 = taylor(f10,x,1, 'Order', 3)

ts10 =

exp(-1) - (exp(-1)*(x - 1)^2)/2

>> ezplot(ts10,[-5 5 -5 5])


Tamara Gascón Moreno
Calculo

Ejercicio 4)
Deduce analíticamente (“a mano”) los dos desarrollos de Taylor del ejercicio anterior que te
correspondan según el segundo y tercer dígito comenzando por la derecha de tu documento
de identidad. En caso de que sean iguales haz el siguiente. Adjunta la resolución manuscrita
escaneada.

Ejercicio 5)
Calcula con MatLab a) la integral indefinida y b) la integral definida que te correspondan según
el 4º dígito comenzando por la derecha de tu documento de identidad.

a) 0∫(1 − 𝑥) ln(1 − 𝑥) 𝑑𝑥

>> syms x % Definimos la variable x como simbólica

>> fi1 = (1 - x)*log(1 - x) % Definimos la función que queremos integrar.

fi1 =

-log(1 - x)*(x - 1)

>> Fi1=int(fi1) % Integramos

Fi1 =

x*(log(1 - x) - 1/2) - x^2*(log(1 - x)/2 - 1/4) - log(x - 1)/2

>> simplify(diff(Fi1)) % Comprobamos resultado

ans =

-log(1 - x)*(x - 1)

1
b) ∫0 (1 − 𝑥) ln(1 − 𝑥) 𝑑𝑥

>> syms x % Definimos la variable x como simbólica

>> a=0; b=1;

>> Fi2=int(fi1,a,b)

Fi2 =

-1/4

Calcula con Matlab las integrales irracionales a) y b) que te correspondan según el 5º dígito
comenzando por la derecha de tu documento de identidad.

𝑥
a) ∫ √1−𝑥 𝑑𝑥
Tamara Gascón Moreno
Calculo

𝑥
b) ∫ √𝑥+1 𝑑𝑥

You might also like