基于TI-83ODE系统求解实验.pdfVIP

LAB:SolvingSystemsofOrdinaryODE

Scenario:

SolvethefollowingsystemofODEs

dx/dt=2x-1.2xy

dy/dt=-y+1.2xy

I.C.y(0)=1,x(0)=1.75

Leth=0.1(StepSize)

TI-83Plus

ClearallcolumnsinL-L

16

1.cedx/dtfunctionintoY

1

\Y=2*X-1.2*X*Y

1

2.cedy/dtfunctionintoY

2

\Y2=-Y+1.2*X*Y

2.EXECtheEULERSYSProgram

DispSystems

DispStepSize

PromptH

DispMAXITER

PromptN

DispICforY,F(0)=C

PromptY

DispICforX,G(0)=D

PromptX

XZ

0T

YW

0I

Lbl2

IfIN

Then

Goto3

：普通ODE求解系统

场景：求解以下ODE系统

dx/dt=2x‑1.2xy

dy/dt=‑y+1.2xy

y(0)=1,x(0)=1.75

设h=0.1（步长）

TI‑83加号

清除L‑L中的所有列

16

1.将dx/dt函数放入Y\Y=2

11

*X‑1.2*X*Y2.将dy/dt函数

放入Y\Y=‑Y+1.2*X*Y

22

2.执行EULERSYS程序

显示“系统”显示

“步长”提示H显示

“MAXITER”提示

N显示“ICfor

Y,F(0)C”提示Y显

示“ICforX,G(0)

D”提示XXZ0

TYW0ILbl2如

果IN那么转到3

Else

Y1J

YK

2

W+H*JW

Z+H*KZ

(T+H)T

IL1(I+1)

XL2(I+1)

YL3(I+1)

WL4(I+1)

ZL(I+1)

5

WY

ZX

I+1I

Goto2

End

Lbl3

DispEND

3.YoucanviewtheoutputinthelistcolumnsofSTAT.

L2isXandL3isY.

4.YoucanplotXvsY(ifyouremovethefunctionsinYandY.)

12