四阶r k matlab程序,四阶R-K公式计算

#include

double fun(double x,double y)

{

double r;

r=y-2*x/y;

return r;

}

void jd_R_K(double xn,double yn,double h,int N)//经典公式

{

double K1,K2,K3,K4;

int i;

cout<

for(i=1;i<=N;i++)

{

K1=fun(xn,yn);

K2=fun(xn+h/2,yn+h/2*K1);

K3=fun(xn+h/2,yn+h/2*K2);

K4=fun(xn+h,yn+h*K3);

xn+=h;

yn+=h/6*(K1+2*K2+2*K3+K4);

cout<

}

cout<

}

void Gill_R_K(double xn,double yn,double h,int N)//吉尔(Gill)公式

{

double K1,K2,K3,K4;

int i;

cout<

for(i=1;i<=N;i++)

{

K1=fun(xn,yn);

K2=fun(xn+h/2,yn+h/2*K1);

K3=fun(xn+h/2,yn+(sqrt(2)-1)/2*h*K1+(1-sqrt(2)/2)*h*K2);

K4=fun(xn+h,yn-sqrt(2)/2*h*K2+(1+sqrt(2)/2)*h*K3);

xn+=h;

yn+=h/6*(K1+(2-sqrt(2))*K2+(2+sqrt(2))*K3+K4);

cout<

}

cout<

}

void main()

{

double a=0,b=1;//定义域；

double x0=0;

double y0=1;

double h=0.1;

int N=int((b-a)/h);

cout<

jd_R_K(x0,y0,h,N);

Gill_R_K(x0,y0,h,N);

}

四阶r k matlab程序,四阶R-K公式计算相关推荐

matlab单位阶跃序列,产生单位阶跃序列的MATLAB程序如下：k.PPT
产生单位阶跃序列的MATLAB程序如下:k 一.实验目的掌握典型连续信号和离散序列的MATLAB描述方法掌握序列的常用运算方法二.实验设备 * * 实验二基本信号的产生及其运算信息工程学院 ...
修正r s分析法matlab程序,关于R/S分析程序用法
看到有人用下面的源代码,但为何有时候求出的值是负数,请各位帮忙,我是新接触最好能稍微具体一点讲解,谢谢. function [logRS,logERS,V]=RSana(x,n,method,q) % ...
matlab 流程计算方法,吸波材料LLG公式计算复磁导率的过程及matlab程序
看到一篇paper,利用Landau-Lifshitz-Gilbert 公式计算片状颗粒的复磁导率.(JAP 107,033913, 2010) http://scitation.aip.org/co ...
四阶龙格库塔法c语言程序,四阶龙格_库塔算法的C语言实现_毋玉芝
2001年3月焦作大学学报 1 第1期JOURNALOFJIAOZUOUNIVERSITYMar.2001 四阶龙格库塔算法的C语言实现毋玉芝 (焦作财会学校) 摘要本文叙述了四阶龙格库塔算 ...
r语言 matlab 函数定义,R语言-编写自定义函数 ZZ
一.函数构造器每一个R函数都包括三个部分:函数名,程序主体以及参数集合,在编写自定义R函数时,需要将三个部分各自储存在一个R对象中.这里需要使用function函数,形如: my_function ...
matlab (k),MATLAB的sgolay（k，f）的Python等价物是什么？
如果您检查Matlab的^{}函数返回的矩阵b,您将看到中心行与SciPy的savgol_coeffs返回的一维数组相同.b的上半部分和下半部分各有(framelen - 1)/2行,是要应用到信号末 ...
Java 1.8（圆的面积和周长）编写程序，使用以下公式计算并显示半径为5.5的圆的面积和周长。
package Try; import java.util.Scanner;//在java.util中导入Scanner包public class Circle {public static void ...
C语言程序用英尺换算公式计算身高的小例子
#include <stdio.h> int main(){//转换公式为(英尺+英寸*12)*0.3048printf("请输入身高的英尺英寸:")float foo ...
有确定项微分方程的matlab程序,微分方程的数值解法matlab四阶龙格—库塔法课件...
<微分方程的数值解法matlab四阶龙格-库塔法课件>由会员分享,可在线阅读,更多相关<微分方程的数值解法matlab四阶龙格-库塔法课件(36页珍藏版)>请在人人文库网上搜索 ...

四阶r k matlab程序,四阶R-K公式计算

四阶r k matlab程序,四阶R-K公式计算相关推荐

最新文章

热门文章