matlab程序设计常微分方程编程例题答案数学.doc

MATLAB程序设计常微分方程1、功能计算函数在指定点的值100,3SIN,XXFTD实验求，，，2F1FFFM文件FUNCTIONFFDXIFX0FINT SINT10 ,0,XELSEXFD2ANS14142FD1ANS1FD0ANS1ANS1FDPIANS63/256PI2、实验分别计算椭圆与的焦距与离心率2153XY21YXM文件FUNCTIONC,EFUNA,BCSQRTA2B2EC/AEND结果C,EFUN25,3C4E08000C,EFUN21,SQRT2C010000IE010000I3、实验分别绘制函数与的图像(要求图像光滑)SIN2YT10SIN56YTM文件FUNCTIONEAA,B,CK2PI/BT0000012KENDYASINBTCPLOTT,Y结果EA1,2,0COLUMNS62793THROUGH628000008000078000760007400072000700006800066COLUMNS62801THROUGH628080006400062000600005800056000540005200050COLUMNS62809THROUGH628160004800046000440004200040000380003600034COLUMNS62817THROUGH628240003200030000280002600024000220002000018COLUMNS62825THROUGH628320001600014000120001000008000060000400002EA10,500PI,PI/6YCOLUMNS1THROUGH85000062932743148386791355965939945299863COLUMNS9THROUGH169781593358866037771566913544644067425882COLUMNS17THROUGH241045305234207913583750000629327431483867COLUMNS25THROUGH324、实验求与180M文件FUNCTIONSQHNS1FORI2NSSIEND结果JC18ANS64024E015JC0ANS15、功能求最大的，使得N1123CN实验分别取，，，求15C012CNM文件FUNCTIONNMAXNCS1I2WHILESACOSBT,IZXAT09944I46A7,0,6,1,0,3,0,5T,IZXAT0I2572．用龙格－库塔方法求下列微分方程初值问题的数值解．1，201YX0YM文件FUNCTIONFLFX,YFY2XFFX,YODE45 LF ,0011,1X001000020000300004000050000600007000080000900010000Y1000011155126421449616755194622266426413307663578841548X,YODE45 LF ,0011,1X001000020000300004000050000600007000080000900010000Y10000111551264214496167551946222664264133076635788415482，，340DXYTT0X1YM文件FUNCTIONDXCZT,XDX3X14X24X13X2TS001PIX00,1OPTIONSODESET RELTOL ,1E6, ABSTOL ,1E9T,XODE45CZ,TS,X0,OPTIONST,XANS10E0040000001000000000100001000000000100001000000000200001000000000300000000010000400002000010000400004000010000300008000010000100011000010000700013000010001500013000010002600008000010003600003000010004400023000010004200052000010002500086000020001400121000020008100143000020017600135000020028900075000020039900059000020046500283000020043000596000020022100967000020023301319000030098401517000030202001369000030323200640000030435400903000030493903411000030434806838000030181110787