Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Lab 3 PDF

Download as pdf or txt
Download as pdf or txt
You are on page 1of 3

In[74]:= Clear[F, y, y0, y1];

borra

In[75]:= y0 = 0; x0 = 0;
F[x_, y_, z_] := 2 z - Cos[x + z] + y Cos[x] (*z=y'*);
coseno coseno
FindRoot[F[x0, y0, y1] ⩵ 0, {y1, 0.5}]
encuentra raíz

Out[77]= {y1 → 0.450184}

In[78]:= F[x0, y0, 0.4501836112948736`]


Out[78]= 1.11022 × 10-16


In[79]:= y0 = 0.450184
Out[79]= 0.450184

In[80]:= Clear[F, y, x, z]; F[x_, y_, z_] := 2 z - Cos[x + z] + y Cos[x] (*z=y'*);


borra coseno coseno
Solve[D[F[x, y[x], y '[x]], x] ⩵ 0, y ''[x]]
resue⋯ deriva
′ - Cos[x] y′ [x] 
Out[81]= y′′ [x] → - Sin[x + y [x]] + Sin[x] y[x]
2 + Sin[x + y′ [x]]

- Sin[x + y′ [x]] + Sin[x] y[x] - Cos[x] y′ [x]


In[82]:= y′′ [x] = sujeta a y[0] = 0, y '[0] = 0.450184;
2 + Sin[x + y′ [x]]

In[82]:= Clear[Y, y, x];


borra
T = 30;
- Sin[x + y′ [x]] + Sin[x] y[x] - Cos[x] y′ [x]
Y[X_] := NDSolvey′′ [x] == &&
resolvedor diferencial numérico 2 + Sin[x + y′ [x]]

y[0] ⩵ 0 && y '[0] ⩵ 0.4501836112948736`, y[x], {x, 0, T}[[1, 1, 2]] //. x → X;


Plot[Evaluate[Y[x]], {x, 0, T}, PlotRange → All, AxesOrigin → {0, 0}]
repr⋯ evalúa rango de rep⋯ todo origen de ejes

0.2

5 10 15 20 25 30

-0.2
Out[83]=

-0.4

-0.6

-0.8

In[84]:= Clear[Y, F, T];


borra
2 lab3.nb

In[85]:= T = 30;

In[86]:= Y[X_];

- Sin[x + y′ [x]] + Sin[x] y[x] - Cos[x] y′ [x]


In[87]:= NDSolvey′′ [x] == && y[0] ⩵ 0 &&
resolvedor diferencial numérico 2 + Sin[x + y′ [x]]

y '[0] ⩵ 0.4501836112948736`, y[x], {x, 0, T}[[1, 1, 2]] //. x → X;


F[x_, y_, z_] := Cos[x] y + 2 z - Cos[x + z];
coseno coseno
lab3.nb 3

In[89]:= R[X_] := F[X, Y[X], Y '[X]];


Table[{X, R[X]}, {X, 0, 10, 0.25}]
tabla

Out[90]= {{0., - Cos[0. + Y′ [0.]] + 1. Y[0.] + 2 Y′ [0.]},


{0.25, - Cos[0.25 + Y′ [0.25]] + 0.968912 Y[0.25] + 2 Y′ [0.25]},
{0.5, - Cos[0.5 + Y′ [0.5]] + 0.877583 Y[0.5] + 2 Y′ [0.5]},
{0.75, - Cos[0.75 + Y′ [0.75]] + 0.731689 Y[0.75] + 2 Y′ [0.75]},
{1., - Cos[1. + Y′ [1.]] + 0.540302 Y[1.] + 2 Y′ [1.]},
{1.25, - Cos[1.25 + Y′ [1.25]] + 0.315322 Y[1.25] + 2 Y′ [1.25]},
{1.5, - Cos[1.5 + Y′ [1.5]] + 0.0707372 Y[1.5] + 2 Y′ [1.5]},
{1.75, - Cos[1.75 + Y′ [1.75]] - 0.178246 Y[1.75] + 2 Y′ [1.75]},
{2., - Cos[2. + Y′ [2.]] - 0.416147 Y[2.] + 2 Y′ [2.]},
{2.25, - Cos[2.25 + Y′ [2.25]] - 0.628174 Y[2.25] + 2 Y′ [2.25]},
{2.5, - Cos[2.5 + Y′ [2.5]] - 0.801144 Y[2.5] + 2 Y′ [2.5]},
{2.75, - Cos[2.75 + Y′ [2.75]] - 0.924302 Y[2.75] + 2 Y′ [2.75]},
{3., - Cos[3. + Y′ [3.]] - 0.989992 Y[3.] + 2 Y′ [3.]},
{3.25, - Cos[3.25 + Y′ [3.25]] - 0.99413 Y[3.25] + 2 Y′ [3.25]},
{3.5, - Cos[3.5 + Y′ [3.5]] - 0.936457 Y[3.5] + 2 Y′ [3.5]},
{3.75, - Cos[3.75 + Y′ [3.75]] - 0.820559 Y[3.75] + 2 Y′ [3.75]},
{4., - Cos[4. + Y′ [4.]] - 0.653644 Y[4.] + 2 Y′ [4.]},
{4.25, - Cos[4.25 + Y′ [4.25]] - 0.446087 Y[4.25] + 2 Y′ [4.25]},
{4.5, - Cos[4.5 + Y′ [4.5]] - 0.210796 Y[4.5] + 2 Y′ [4.5]},
{4.75, - Cos[4.75 + Y′ [4.75]] + 0.0376022 Y[4.75] + 2 Y′ [4.75]},
{5., - Cos[5. + Y′ [5.]] + 0.283662 Y[5.] + 2 Y′ [5.]},
{5.25, - Cos[5.25 + Y′ [5.25]] + 0.512085 Y[5.25] + 2 Y′ [5.25]},
{5.5, - Cos[5.5 + Y′ [5.5]] + 0.70867 Y[5.5] + 2 Y′ [5.5]},
{5.75, - Cos[5.75 + Y′ [5.75]] + 0.861192 Y[5.75] + 2 Y′ [5.75]},
{6., - Cos[6. + Y′ [6.]] + 0.96017 Y[6.] + 2 Y′ [6.]},
{6.25, - Cos[6.25 + Y′ [6.25]] + 0.999449 Y[6.25] + 2 Y′ [6.25]},
{6.5, - Cos[6.5 + Y′ [6.5]] + 0.976588 Y[6.5] + 2 Y′ [6.5]},
{6.75, - Cos[6.75 + Y′ [6.75]] + 0.893006 Y[6.75] + 2 Y′ [6.75]},
{7., - Cos[7. + Y′ [7.]] + 0.753902 Y[7.] + 2 Y′ [7.]},
{7.25, - Cos[7.25 + Y′ [7.25]] + 0.567924 Y[7.25] + 2 Y′ [7.25]},
{7.5, - Cos[7.5 + Y′ [7.5]] + 0.346635 Y[7.5] + 2 Y′ [7.5]},
{7.75, - Cos[7.75 + Y′ [7.75]] + 0.103794 Y[7.75] + 2 Y′ [7.75]},
{8., - Cos[8. + Y′ [8.]] - 0.1455 Y[8.] + 2 Y′ [8.]},
{8.25, - Cos[8.25 + Y′ [8.25]] - 0.385748 Y[8.25] + 2 Y′ [8.25]},
{8.5, - Cos[8.5 + Y′ [8.5]] - 0.602012 Y[8.5] + 2 Y′ [8.5]},
{8.75, - Cos[8.75 + Y′ [8.75]] - 0.780846 Y[8.75] + 2 Y′ [8.75]},
{9., - Cos[9. + Y′ [9.]] - 0.91113 Y[9.] + 2 Y′ [9.]},
{9.25, - Cos[9.25 + Y′ [9.25]] - 0.984765 Y[9.25] + 2 Y′ [9.25]},
{9.5, - Cos[9.5 + Y′ [9.5]] - 0.997172 Y[9.5] + 2 Y′ [9.5]},
{9.75, - Cos[9.75 + Y′ [9.75]] - 0.94758 Y[9.75] + 2 Y′ [9.75]},
{10., - Cos[10. + Y′ [10.]] - 0.839072 Y[10.] + 2 Y′ [10.]}}

You might also like