Быстрое введение в систему Mathematica. Ч.1. Израилевич Я.А - 24 стр.

UptoLike

Рубрика: 

Ke_^mxsZy kljhdZ ij_^klZ\ey_l j_ZebaZpbx f_lh^Z GvxlhgZ qbke_ggh]h
hlukdZgby dhjgy mjZ\g_gby/ \uiheg_ggmx \ klbe_ nmgdpbhgZevgh]h
ijh]jZffbjh\Zgby=
In[9]:=newton[f_,x0_]:=FixedPoint[(#-f[#]/f'[#])&,N[x0]]
A^_kv
&
ij_^klZ\ey_l kh[hc nhjfZevguc Zj]mf_gl \uibkZgghc \ kdh[dZo
Zghgbfghc +[_aufygghc, nmgdpbb/ hkms_kl\eyxs_c rZ] bl_jZpbb/
)
hlf_qZ_l
dhg_p aZibkb wlhc Zghgbfghc nmgdpbb/ nmgdpby
)L[HG3RLQW>K/H[SU@
hkms_kl\ey_l ihke_^h\Zl_evgh_ ijbf_g_gb_ +kmi_jihabpbx, nmgdpbb K/
gZqbgZy k \ujZ`_gby
H[SU
/ ihdZ ihemqZ_fu_ agZq_gby g_ i_j_klZgml
baf_gylvky/ nmgdpby
1
h[_ki_qb\Z_l \uiheg_gb_ qbke_gguo +g_ kbf\hevguo,
\uqbke_gbc1 >Zggh_ hij_^_e_gb_ iha\hey_l bkdZlv dhjgb mjZ\g_gbc=
In[10]:=g[x_]:=x^3-2x
In[11]:= newton[g,{0,2,-2}]
Out[11]= {0, 1.41421, -1.41421}
Fu \b^bf/ qlh \ klbe_ nmgdpbhgZevgh]h ijh]jZffbjh\Zgby \Z`g_crmx
jhev b]jZxl kihkh[u hj]ZgbaZpbb ihke_^h\Zl_evgh]h ijbf_g_gby
+kmi_jihabpbc, nmgdpbc/ ^ey q_]h \ kbkl_f_
³0DWKHPDWLFD´
bf__lky [h]Zluc
gZ[hj kj_^kl\1
Hq_gv \Z`gZ jZagbpZ f_`^m g_f_^e_gguf ijbk\Zb\Zgb_f
OHIW ULJKW
/
dhlhjh_ \uihegy_lky \ fhf_gl \uiheg_gby ijbk\Zb\Zgby/ b h leh`_gguf
ijbk\Zb\Zgb_f
OHIW= ULJKW/
dhlhjh_ \uihegy_lky dZ`^uc jZa/ dh]^Z
bkihevam_lky agZq_gb_
OHIW
1 Wlh ohjhrh \b^gh \ ke_^mxs_f ijbf_j_/ ]^_
5DQGRP>@
± dhfZg^Z ± ^Zlqbd kemqZcguo qbk_e1
In[12]:= x:=Random[]
In[13]:= y=Random[]
Out[13]= 0.199496
In[14]:= {x,x,y,y}
Out[14]= {0.784326, 0.453175, 0.199496, 0.199496}
;hevrb_ \hafh`ghklb ij_^hklZ\ey_l ijh]jZffbjh\Zgb_/ hkgh\Zggh_ gZ
ijZ\beZo ij_h[jZah\Zgbc1 Ijbf_g_gb_ ijZ\beZ ij_h[jZah\Zgby +beb ijZ\be,
UXOH
d \ujZ`_gbx
H[SU
h[hagZqZ_lky dZd
H[SU21UXOH
1
Ijbf_jhf ijZ\beZ ij_h[jZah\Zgby fh`_l kem`blv h[uqgZy ih^klZgh\dZ=
In[1]:= x^4-a^4
4 4
Out[1]= -a + x
In[2]:= %/.x->a
Out[2]= 0
?s_ h^gh ihe_agh_ kj_^kl\h ± mkeh\b_/ h[hagZqZ_fh_ dZd
H[SU2>FRQG
In[2]:= delta[i_Integer,j_Integer]:=0/;i!=j
delta[i_Integer,j_Integer]:=1/;i==j
In[3]:= {delta[2,3], delta[1,1], delta[2.5,3]}
     ��������� ������ ������������ ���������� ������ ������� ����������
��������� ����� ���������� ����������� � ����� ���������������
�����������������

In[9]:=newton[f_,x0_]:=FixedPoint[(#-f[#]/f'[#])&,N[x0]]

����� � ������������ ����� ���������� �������� ���������� � �������
��������� ������������ �������� �������������� ��� ��������� � ��������
����� ������ ���� ��������� �������� ������� ������������������
������������ ���������������� ���������� �������������� ������� ��
������� � ��������� ����� ���� ���������� �������� �� ����������
����������� ������� � ������������ ���������� ��������� ��� �����������
����������� ������ ����������� ��������� ������ ����� ����������
In[10]:=g[x_]:=x^3-2x
In[11]:= newton[g,{0,2,-2}]
Out[11]= {0, 1.41421, -1.41421}

      �� ������ ��� � ����� ��������������� ���������������� ���������
���� ������ ������� ����������� ����������������� ����������
�������������� �������� ��� ���� � ������� ������������� ������� �������
����� ��������
      ����� ����� ������� ����� ����������� ������������� �����������
������� ����������� � ������ ���������� ������������� � ����������
������������� ������������ ������� ����������� ������ ���� �����
������������ �������� ����� ��� ������ ����� � ��������� �������� ���
�������� � ������� � ������ ��������� ������
      In[12]:= x:=Random[]
      In[13]:= y=Random[]
      Out[13]= 0.199496
      In[14]:= {x,x,y,y}
      Out[14]= {0.784326, 0.453175, 0.199496, 0.199496}



  ������� ����������� ������������� ����������������� ���������� ��
�������� ��������������� ���������� ������� �������������� ���� �������
���� � ��������� ���� ������������ ��� �����������
�������� ������� �������������� ����� ������� ������� ������������
In[1]:= x^4-a^4
             4     4
Out[1]= -a + x
In[2]:= %/.x->a
Out[2]= 0
��� ���� �������� �������� � �������� ������������ ��� ����������
In[2]:= delta[i_Integer,j_Integer]:=0/;i!=j
           delta[i_Integer,j_Integer]:=1/;i==j
In[3]:= {delta[2,3], delta[1,1], delta[2.5,3]}