欧拉法解常微分方程代码

来源：爱go旅游网

// §4 常微分方程的差分方法

// 程序4.1 — 欧拉方法

#include

#include \"expressi.cpp\"

void Euler(char FxyString[],float x0,float y0,float h,int n)

{

int i;

float yi;

if(CreateFxy(FxyString)) return;

yi=y0;

for(i=1;i{

yi=yi+h*f(x0+(i-1)*h,yi);

printf(\"\\ny%d=%f\

}

void main()

{

float x0,y0,h; //x0:自变量初值,y0:y(x0),h:步长

int n; //步数

char FxyString[200]; //存放C表达式格式的f(x,y)

printf(\"\\nInput function,x0,y0,h,n: \");

scanf(\"%s %f %f %f %d\

Euler(FxyString,x0,y0,h,n);

getch();

}

运行实例

Input function,x0,y0,h,n: y-2*x/y 0 1 0.1 11

y1=1.100000

y2=1.191818

y3=1.277438

y4=1.358213

y5=1.435133

y6=1.5066

y7=1.580338

y8=1.9784

y9=1.717780

y10=1.784771

因篇幅问题不能全部显示，请点此查看更多更全内容

查看全文