Book a Demo!
CoCalc Logo Icon
StoreFeaturesDocsShareSupportNewsAboutPoliciesSign UpSign In
Download

Jupyter notebook Untitled1.ipynb

298 views
Tweet
Share
Share
Kernel: Python 3 (Ubuntu Linux)
%matplotlib inline
from pandas import DataFrame
from IPython.display import display,Math,clear_output
from numpy import array,linspace,zeros
from sympy import latex
import matplotlib.pyplot as plt
from IPython.display import HTML

HTML('''<script>
code_show=true; 
function code_toggle() {
 if (code_show){
     $('div.input').hide();
 } else {
     $('div.input').show();
 }
 code_show = !code_show
} 
$( document ).ready(code_toggle);
</script>
<form action="javascript:code_toggle()"><input type="submit" value="Click here to toggle on/off the raw code."></form>''')

задача 1(118 а)

xy dx−(1+x2)dy=0xy dx=(1+x2)dyx dx1+x2=dyy∫x dx1+x2=∫dyyпроизводя интегрирование  ∫x dx1+x2=12∫dx21+x2=12ln(1+x2)  ∫dyy=ln(y)+Cполучимln(y)+C=12ln(1+x2)или собирая логарифмы (где C=−ln(C1)):ln(y)=ln(C1(1+x2))y=C1(1+x2)xy\ d x- (1+x^2) d y=0\\ x y\ d x = (1+x^2) d y \\ \frac{x\ dx}{1+x^2}=\frac{dy}{y}\\ \int \frac{x\ dx}{1+x^2}=\int \frac{dy}{y}\\ производя\ интегрирование \\ \ \ \int \frac{x\ dx}{1+x^2}=\frac{1}{2} \int \frac{dx^2}{1+x^2}=\frac{1}{2} ln(1+x^2)\\ \ \ \int \frac{dy}{y} =ln(y)+C \\ получим \\ ln(y)+C=\frac{1}{2} ln(1+x^2) \\ или\ собирая\ логарифмы\ (где\ C=-ln(C_{1})): \\ ln(y)=ln(C_{1} \sqrt(1+x^2))\\ y=C_{1} \sqrt(1+x^2)

Задача 2 (118 б)

y′=2xy+y2, y(2)=1, h=0.2по формуле:yi=yi−1+(xi−xi−1) f(xi−1,yi−1) начиная с x0=2,y0=1 получимy'= \frac{2 x}{y}+\frac{y}{2},\ y(2)=1,\ h=0.2\\ по\ формуле:\\ y_{i}=y_{i-1}+(x_{i}-x_{i-1})\ f(x_{i-1},y_{i-1})\ начиная\ с\ x_{0}=2, y_{0}=1\ получим \\ 

x=linspace(2,3,6)
y=zeros((6))
f=lambda x,y: (2*x/y+y/2.)
index=["i=0"]
columns=[]
y[0]=1
for i in range(len(y))[1:]:    
    y[i]=y[i-1]+(x[i]-x[i-1])*(f(x[i-1],y[i-1]))
    index.append("i=%s"%(i))
#print(y)
d=DataFrame(array([x,y]),columns=index,index=[r'$'+latex("x_{i}")+'$',r'$'+latex("y_{i}")+'$'])#
display(d)
plt.plot(x,y,'o')
[<matplotlib.lines.Line2D at 0x7f4b69079eb8>]
Image in a Jupyter notebook

Задача 3( 118 в)

(x2+y−yex) dx+(x+2y−ex) dy=0 (1)т.к. (x2+y−yex)y=1−ex=(x+2y−ex)x,то выражение (1) являетсяполным дифференциалом некоторой функции F(x,y) т.е:dF(x,y)=Fx(x,y)dx+Fy(x,y)dyгде  Fx(x,y)=(x2+y−yex)  Fy(x,y)=(x+2y−ex)(т.е. выполняется Fxy(x,y)=Fyx(x,y) необходимое усолвие дляполного дифференциала)а т.к. dF(x,y)=0 тоF(x,y)=const (2)−неявное решение уравнения (1)для нахождения F(x,y) нужно интегрировать одну из частных производныхпо соответсвующему аргуметну. Например по x получим:F(x,y)=∫Fx(x,y)dx=∫(x2+y−yex)dx=x33+xy−yex+C(y)   (3)здесь С(у)−постоянная интегрирования (не зависит от x но зависит от y).Для нахождения C(y) надо продифференцировать (3) по y и прировнять кFy(x,y)=(x+2y−ex)получим дифф уравнение:C′(y)+x−ex=(x+2y−ex)C′(y)=2ydC(y)dy=2ydC(y)=2ydy∫dC(y)=∫2ydyC(y)=y2+CПодставляя C(y) в (3) найдем F(x,y) с точностью до константы:F(x,y)=x33+xy−yex+y2+Cа затем, воспользовавшись (2) получим решение в неявной форме, с точностью до константы const:x33+xy−yex+y2=const(x^2+y-y e^x)\ d x +(x+2 y-e^x)\ d y=0\ (1)\\ т.к.\ (x^2+y-y e^x)_{y}=1-e^x=(x+2 y-e^x)_{x},то\ выражение\ (1)\ является \\ полным\ дифференциалом\ некоторой\ функции\ F(x,y)\ т.е:\\ d F(x,y)=F_{x}(x,y) d x + F_{y}(x,y) d y \\ где\\ \ \ F_{x}(x,y)=(x^2+y-y e^x)\\ \ \ F_{y}(x,y)=(x+2 y-e^x)\\ (т.е.\ выполняется\ F_{x y}(x,y)=F_{y x}(x,y)\ необходимое\ усолвие\ для\\ полного\ дифференциала ) \\ а\ т.к.\ d F(x,y)=0\ то\\ F(x,y)=const\ (2)\\ -неявное\ решение\ уравнения\ (1)\\ для\ нахождения\ F(x,y)\ нужно\ интегрировать\ одну\ из\ частных\ производных\\ по\ соответсвующему\ аргуметну.\ Например\ по\ x\ получим:\\ F(x,y)=\int F_{x}(x,y) d x = \int (x^2+y-y e^x) d x= \frac{x^3}{3}+x y-y e^x+C(y) \ \ \ (3)\\ здесь\ С(у)-постоянная\ интегрирования\ (не\ зависит\ от\ x\ но\ зависит\ от\ y).\\ Для\ нахождения\ C(y)\ надо\ продифференцировать\ (3)\ по\ y\ и\ прировнять\ к\\ F_{y}(x,y)=(x+2 y-e^x)\\ получим\ дифф\ уравнение:\\ C'(y)+x-e^x=(x+2 y-e^x)\\ C'(y)=2 y\\ \frac{d C(y)}{d y}=2 y\\ d C(y)=2 y d y\\ \int d C(y)=\int 2 y d y\\ C(y)=y^2+C\\ Подставляя\ C(y)\ в\ (3)\ найдем\ F(x,y)\ с\ точностью\ до\ константы:\\ F(x,y)=\frac{x^3}{3}+x y-y e^x+y^2+C\\ а\ затем,\ воспользовавшись\ (2)\ получим\ решение\ в\ неявной\ форме,\ с\ точностью\ до\ константы\ const:\\ \frac{x^3}{3}+x y-y e^x+y^2=const

Задача 4 (127)

В партии из 10 деталей 8 стандартных. Выбрали 2-е детали. Найти мат одидание (E(X)E(X)) и дисперсию (D(X)D(X)) числа стандартных деталий Решение:E(X)=∑(xip(X=xi))=0 p(0)+1 p(1)+2 p(2)p(1) можно найти как число способов выбрать 1 елемент из 8умножить на число способов выбрать 1 елемент из2деленное на число способов выбрать 2 елемента из 10:p(1)=(82)(21)(102)=281029p(2) можно найти как число способов выбрать 2 елемента из 8деленное на число способов выбрать 2 елемента из 10:p(2)=(82)(102)=8!(8−2)!2!10!(10−2)!2!=810∗79т.о.E(X)=1 p(1)+2 p(2)=1.6D(X)=E(X2)−E2(x)=(12p(1)+22p(2))−1,62=0.284Решение:\\ E(X)=\sum(x_{i}p(X=x_{i}))=0\ p(0)+1\ p(1)+2\ p(2)\\ p(1)\ можно\ найти\ как\ число\ способов\ выбрать\ 1\ елемент\ из\ 8\\ умножить\ на\ число\ способов\ выбрать\ 1\ елемент\ из 2\\ деленное\ на\ число\ способов\ выбрать\ 2\ елемента\ из\ 10:\\ p(1)=\frac{ \begin{pmatrix} 8\\ 2 \end{pmatrix} \begin{pmatrix} 2\\ 1 \end{pmatrix} } { \begin{pmatrix} 10\\ 2 \end{pmatrix} }=2 \frac{8}{10}\frac{2}{9}\\ p(2)\ можно\ найти\ как\ число\ способов\ выбрать\ 2\ елемента\ из\ 8\\ деленное\ на\ число\ способов\ выбрать\ 2\ елемента\ из\ 10:\\ p(2)=\frac{ \begin{pmatrix} 8\\ 2 \end{pmatrix}} { \begin{pmatrix} 10\\ 2 \end{pmatrix} }=\frac{\frac{8!}{(8-2)! 2!}}{\frac{10!}{(10-2)! 2!}}=\frac{8}{10}*\frac{7}{9}\\ т.о.\\ E(X)=1\ p(1)+2\ p(2)=1.6\\ D(X)=E(X^2)-E^2(x)=(1^2 p(1)+2^2 p(2))-1,6 ^2=0.284