// §4 常微分方程的差分方法
// 程序4.1 — 欧拉方法
#include <stdlib.h>
#include <stdio.h>
#include <conio.h>
#include <string.h>
#include <graphics.h>
#include <math.h>
#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<n;++i)
{
yi=yi+h*f(x0+(i-1)*h,yi);
printf("\ny%d=%f",i,yi);
}
}
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",FxyString,&x0,&y0,&h,&n);
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.508966
y7=1.580338
y8=1.649784
y9=1.717780
y10=1.784771
*/
