matlab4阶龙格库塔法程序

上传人:沈*** 文档编号:109575267 上传时间:2022-06-17 格式:DOC 页数:5 大小:127KB
收藏 版权申诉 举报 下载
matlab4阶龙格库塔法程序_第1页
第1页 / 共5页
matlab4阶龙格库塔法程序_第2页
第2页 / 共5页
matlab4阶龙格库塔法程序_第3页
第3页 / 共5页
资源描述:

《matlab4阶龙格库塔法程序》由会员分享,可在线阅读,更多相关《matlab4阶龙格库塔法程序(5页珍藏版)》请在装配图网上搜索。

1、Four short words sum up what has lifted most successful individuals above the crowd: a little bit more.-author-datematlab4阶龙格库塔法程序matlab编的4阶龙格库塔法解微分方程的程序function varargout=saxplaxliu(varargin)clc,clearx0=0;xn=1.2;y0=1;h=0.1;y,x=lgkt4j(x0,xn,y0,h); n=length(x);fprintf( i x(i) y(i)n);for i=1:n fprintf

2、(%2d %4.4f %4.4fn,i,x(i),y(i);endfunction z=f(x,y) z=-2*x*y2;function y,x=lgkt4j(x0,xn,y0,h) x=x0:h:xn;n=length(x);y1=x;y1(1)=y0;for i=1:n-1 K1=f(x(i),y1(i); K2=f(x(i)+h/2,y1(i)+h/2*K1); K3=f(x(i)+h/2,y1(i)+h/2*K2); K4=f(x(i)+h,y1(i)+h*K3); y1(i+1)=y1(i)+h/6*(K1+2*K2+2*K3+K4);end y=y1;结果:i x(i) y(i)1 0.0000 1.00002 0.1000 0.99013 0.2000 0.96154 0.3000 0.91745 0.4000 0.86216 0.5000 0.80007 0.6000 0.73538 0.7000 0.67119 0.8000 0.609810 0.9000 0.552511 1.0000 0.500012 1.1000 0.452513 1.2000 0.4098-

展开阅读全文
温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

copyright@ 2023-2025  zhuangpeitu.com 装配图网版权所有   联系电话:18123376007

备案号:ICP2024067431-1 川公网安备51140202000466号


本站为文档C2C交易模式,即用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知装配图网,我们立即给予删除!