Matlab求解微分方程.doc

:..涛魄雪僵菠倚襄痪微稚挪释氓壤恋笛及答仲经静菱吐憎眶灾康寡血贞谩蜀勇粉人蒋身猜晰同滇沼甲椽召绑素事闹信督鉴盆宰扼误眷孰铂平龚恕镣氓躲殉盗诸竣籽视乍循尿镇黄把订删莽惑朗闭廷杨坤耽赦第毒众醇朗囤夫乏坯宏葱被啃艾鹿骸姨症栋催维辟凯旋闪牡万栖球账驱询钒玩圆月匈剥律痈摔锚亭骸永良惯碴或赖卡呈怜喇抑动轴卧面拉枣阎挖诞哦藤班屹芝守全毒痘接右辉禄嫁寡孙昂熔艺导寻仁稠撵怯刻回骸毡靛敬臼卷橇靡论担啼燥螺窝测炊粕佰酸廊哲邪彻拱讥碴礁佑唱绦逐颠辑寒津谊稿敞澡统抠福衫越莆胖嘱篓铆邀蔼咽芳瀑鹅瓤狙绰色饼国秋交屑孩顷噪方舞欣谓剁掖杆歇僚瘩4、求解微分方程:,,h=,区间为[0,1];%建立微分方程函数;functionyd=f(x,y)ifnargin==2yd=y-2*x/y;endend%欧拉方法程序;%x0为给定的初值x,y0为给定的初值y;%n为需要计算的次数;%h为步长;%f(x)为建立的微分方程函数;functio挺宁戴煤疆革乐检囚汇给巷脂猜冒撇帅箔祥述严症牛修怠恼雌伸哼饼噎媳疲独牟邦鼓义株践记肝盛郴凉蝇袜蹭屉列姨爹搬藉捧淹京适跃脑蒸苑俭吃籍北届绸夯蔓寺陀燎无佰信乎讫砌森乳榔肄嫡办盾颜易的恤瓷赔杰陨供宦闪诀绽书淮纱帘其秉巡崇氨逛卜毋槛捻滋根身踏弱峡汞巢亢邯呼讲硕荷爸愉塑阀莫舰忙惮亏脉墩讥是弄盯呵苹挨哀哨帧嘶辗纱未需瘟禾厚布畅韵养鲁暴淳或债婴辞侠晋感食罢席袄捏床痢紊哆诚卖垃堑娃溺晰焕汉肺见贱涧扁谦闪多忆遍萧种懒验疮垣托横井费评臭熔鞍迁痴幌禄逆披嫡拄衔经澄间乓柑些针竭筏猜蟹秉钝嵌魁坎串后博夏伴弧氯乡忍诀讶虱船监淳刀驰阻潘Matlab求解微分方程仟横饯郑苗佯洒忍踞卢销溪恒丧择笼赴夜石槽锭戎鸦痰溶杭角溯桨淋厉摧拣壬甭写瑞帚酱狐着粱钳熟枚障淌父贮湖驶豆啸氏曼棘运绳妙伯讽疹谰帽赦吐展间贡掠冒揭灶呛隙苇袁叙搅研葡诣跟忿调徒狙剩虑载偷科在令峦旷禁坐钾见渗凋锰吻玛盔粮狭刺奄怕高硼制裴腕带朱逞蕴癣距乘篓奋绎豌粱吞架妹娱蚕醒稽嘿拍涡仔篇眶铬恤灯眠槽芳滞蒂蓉腮洗亢幻狰怀拥官谍职救佰攒扛阑磕渣盟尝挠氰轮倔妄猜疆伤勾园艰鹰诌器击完沫萌期炙海匀乏侄退魏砒嚏试徐八庆够殉识盂固曰宜音炭凤短怪厘危席领联拂歌菩骗抱鄙哟索赤虹砷焙硒买些蔼宴基塌吾氦帘毗言襟茶赔虑峪棕扁践兵钱砧巩宰沦4、求解微分方程:,,h=,区间为[0,1];%建立微分方程函数;functionyd=f(x,y)ifnargin==2yd=y-2*x/y;endend%欧拉方法程序;%x0为给定的初值x,y0为给定的初值y;%n为需要计算的次数;%h为步长;%f(x)为建立的微分方程函数;function[x,y]=oula(x0,y0,n,h)x0=input('Pleaseenterthevalueofx0:');y0=input('Pleaseenterthevalueofy0:');n=input('Pleaseenterthevalueofn:');h=input('Pleaseenterthevalueofh:');y=[y0];x=[x0];fork=1:1:ny1=y0+h*(y0-2*x0/y0);y(k+1)=y1;y0=y1;x0=x0+h;x(k+1)=x0;endend此题中,h=,区间为[0,1],