Практикум по предмету Математические методы и модели
Категория реферата: Рефераты по математике
Теги реферата: бесплатные конспекты, решебник по английскому
Добавил(а) на сайт: Shereshevskij.
Предыдущая страница реферата | 4 5 6 7 8 9 10 11 12 13 14 | Следующая страница реферата
Таблица П2. t-распределение Стьюдента
| |Значения t |
|( |При (=0,1 |При (=0,05 |
|10 |1,812 |2,228 |
|11 |1,796 |2,201 |
|12 |1,782 |2,179 |
|13 |1,771 |2,160 |
|14 |1,761 |2,145 |
Приложение 2
Текст программы численного решения системы семи дифференциальных уравнений
Sub DU()
x1=1 'начальные условия при t=0
x2=0 'начальные условия при t=0
x3=0 'начальные условия при t=0
x4=0 'начальные условия при t=0
x5=0 ' начальные условия при t=0
x6=0 ' начальные условия при t=0
x7=0 ' начальные условия при t=0
Sheets("1").Cells(k+2;2).Value=x1
Sheets("1").Cells(k+2;3).Value=x2
Sheets("1").Cells(k+2;4).Value=x3
Sheets("1").Cells(k+2;5).Value=x4
Sheets("1").Cells(k+2;6).Value=x5
Sheets("1").Cells(k+2;7).Value=x6
Sheets("1").Cells(k+2;8).Value=x7
dt=30/50 a12=Sheets("1").Cells(5;9).Value ' инт. потока a13=Sheets("1").Cells(5;10).Value ' инт. потока a21=Sheets("1").Cells(5;11).Value ' инт. потока a23=Sheets("1").Cells(5;12).Value ' инт. потока a34=Sheets("1").Cells(5;13).Value ' инт. потока a45=Sheets("1").Cells(5;14).Value ' инт. потока a52=Sheets("1").Cells(5;15).Value ' инт. потока a26=Sheets("1").Cells(5;16).Value ' инт. потока a62=Sheets("1").Cells(5;17).Value ' инт. потока a67=Sheets("1").Cells(5;18).Value ' инт. потока a72=Sheets("1").Cells(5;19).Value ' инт. потока
For k = 0 To 50 k1=One(x1;x2;x3;x4;x5;x6;x7;a12;a13;a21)*dt m1=Two(x1;x2;x3;x4;x5;x6;x7;a12;a26;a21;a23;a52;a62;a72)*dt n1=Three(x1;x2;x3;x4;x5;x6;x7;a13;a23;a34)*dt o1=Four(x1;x2;x3;x4;x5;x6;x7;a34;a45)*dt p1=Five(x1;x2;x3;x4;x5;x6;x7;a45;a52)*dt r1=Six(x1;x2;x3;x4;x5;x6;x7;a26;a67;a62)*dt s1=Seven(x1;x2;x3;x4;x5;x6;x7;a67;a72)*dt k2=One(x1+0,5*k1;x2+0,5*m1;x3+0,5*n1;x4+0,5*o1;x5+0,5*p1; x6+0,5*r1;x7+0,5*s1;a12;a13;a21)*dt m2=Two(x1+0,5*k1;x2+0,5*m1;x3+0,5*n1;x4+0,5*o1;x5+0,5*p1; x6+0,5*r1;x7+0,5*s1;a12;a26;a21;a23;a52;a62;a72)*dt n2=Three(x1+0,5*k1;x2+0,5*m1;x3+0,5*n1;x4+0,5*o1;x5+0,5*p1;x6+0,5*r1;x7
+0,5*s1;a13;a23;a34)*dt o1=Four(x1+0,5*k1;x2+0,5*m1;x3+0,5*n1;x4+0,5*o1;x5+0,5*p1;x6+0,5*r1;x7+
0,5*s1;a34;a45)*dt p1=Five(x1+0,5*k1;x2+0,5*m1;x3+0,5*n1;x4+0,5*o1;x5+0,5*p1;x6+0,5*r1;x7+
0,5*s1;a45;a52)*dt r1=Six(x1+0,5*k1;x2+0,5*m1;x3+0,5*n1;x4+0,5*o1;x5+0,5*p1;x6+0,5*r1;x7+0
,5*s1;a26;a67;a62)*dt s1=Seven(x1+0,5*k1;x2+0,5*m1;x3+0,5*n1;x4+0,5*o1;x5+0,5*p1;x6+0,5*r1;x7
+0,5*s1;a67;a72)*dt
k3=One(x1+0,5*k2;x2+0,5*m2;x3+0,5*n2;x4+0,5*o2;x5+0,5*p2;x6+0,5*r2;x7+0,5*s2
;a12;a13;a21)*dt
m3=Two(x1+0,5*k2;x2+0,5*m2;x3+0,5*n2;x4+0,5*o2;x5+0,5*p2;x6+0,5*r2;x7+0,5*s2
;a12;a26;a21;a23;a52;a62;a72)*dt
n3=Three(x1+0,5*k2;x2+0,5*m2;x3+0,5*n2;x4+0,5*o2;x5+0,5*p2;x6+0,5*r2;x7+0,5*
s2;a13;a23;a34)*dt
o3=Four(x1+0,5*k2;x2+0,5*m2;x3+0,5*n2;x4+0,5*o2;x5+0,5*p2;x6+0,5*r2;x7+0,5*s
2;a34;a45)*dt
p3=Five(x1+0,5*k2;x2+0,5*m2;x3+0,5*n2;x4+0,5*o2;x5+0,5*p2;x6+0,5*r2;x7+0,5*s
2;a45;a52)*dt
r3=Six(x1+0,5*k2;x2+0,5*m2;x3+0,5*n2;x4+0,5*o2;x5+0,5*p2;
x6+0,5*r2;x7+0,5*s2;a26;a67;a62)*dt
s3=Seven(x1+0,5*k2;x2+0,5*m2;x3+0,5*n2;x4+0,5*o2;x5+0,5*p2;x6+0,5*r2;x7+0,5*
s2;a67;a72)*dt k4=One(x1+k3;x2+m3;x3+n3;x4+o3;x5+p3;x6+r3;x7+s3;a12;a13;a21)*dt m4=Two(x1+k3;x2+m3;x3+n3;x4+o3;x5+p3;x6+r3;x7+s3;a12;a26;a21;a23;a52;a6
2;a72)*dt n4=Three(x1+k3;x2+m3;x3+n3;x4+o3;x5+p3;x6+r3;x7+s3;a13;a23;a34)*dt o4=Four(x1+k3;x2+m3;x3+n3;x4+o3;x5+p3;x6+r3;x7+s3;a34;a45)*dt p4=Five(x1+k3;x2+m3;x3+n3;x4+o3;x5+p3;x6+r3;x7+s3;a45;a52)*dt r4=Six(x1+k3;x2+m3;x3+n3;x4+o3;x5+p3;x6+r3;x7+s3;a26;a67;a62)*dt s4=Seven(x1+k3;x2+m3;x3+n3;x4+o3;x5+p3;x6+r3;x7+s3;a67;a72)*dt x1=x1+(k1+2*k2+2*k3+k4)/6 x2=x2+(m1+2*m2+2*m3+m4)/6 x3=x3+(n1+2*n2+2*n3+n4)/6 x4=x4+(o1+2*o2+2*o3+o4)/6 x5=x5+(p1+2*p2+2*p3+p4)/6 x6=x6+(r1+2*r2+2*r3+r4)/6 x7=x7+(s1+2*s2+2*s3+s4)/6
Sheets("1").Cells(k+3;2).Value=x1
Sheets("1").Cells(k+3;3).Value=x2
Sheets("1").Cells(k+3;4).Value=x3
Sheets("1").Cells(k+3;5).Value=x4
Sheets("1").Cells(k+3;6).Value=x5
Sheets("1").Cells(k+3;7).Value=x6
Sheets("1").Cells(k+3;8).Value=x7
Next
End Sub
Function One(x1;x2;x3;x4;x5;x6;x7;a12;a13;a21)'Вер.P1
One=-(a12+a13)*x1+a21*x2
End Function
FunctionTwo(x1;x2;x3;x4;x5;x6;x7;a12;a26;a21;a23;a52;a62;a72)'Вер.P4
Two=a12*x1-(a26+a21+a23)*x2+a52*x5+a62*x6+a72*x7
End Function
Function Three(x1;x2;x3;x4;x5;x6;x7;a13;a23;a34)'Вер.P3
Three=a13*x1+a23*x2-a34*x3
End Function
Function Four(x1;x2;x3;x4;x5;x6;x7;a34;a45)'Вер.Р4
Four=a34*x3-a45*x4
End Function
Function Five(x1;x2;x3;x4;x5;x6;x7;a45;a52)'Вер.Р5
Five=a45*x4-a52*x5
End Function
Function Six(x1;x2;x3;x4;x5;x6;x7;a26;a67;a62)'Вер.Р6
Six=a26*x2-(a67+a62)*x6
End Function
Function Seven(x1;x2;x3;x4;x5;x6;x7;a67;a72)'Вер.Р7
Seven=a67*x6-a72*x7
End Function
Приложение 3
Рекомендуем скачать другие рефераты по теме: диплом, оформление доклада титульный лист.
Предыдущая страница реферата | 4 5 6 7 8 9 10 11 12 13 14 | Следующая страница реферата