euler.c

来自「【欧拉算法】 微分方程的本质特征是方程中含有导数项,数值解法的第一步就是...欧」· C语言 代码 · 共 42 行

#include "stdio.h"
#include "math.h"
#define N 100
main()
{
    int n,i;
    double x[N],y[N],ey[N],wc[N],a,b,h;
    
    printf("/*************************************\\\n");
    printf("       Use eular solove question!\n");
    printf("input y(1)=");
    scanf("%lf",&y[0]);
    printf("\ninput the rang a,b:");
    scanf("%lf,%lf",&a,&b);
    x[0]=a;
	ey[0]=y[0];
	wc[0]=y[0]-ey[0];
    printf("\ninput n(n<=N)=");
    scanf("%d",&n);
    h=(b-a)/n;
    for(i=1;i<n+1;i++)
		x[i]=x[0]+i*h;
    for(i=1;i<n+1;i++)
	{
    ey[i]=ey[i-1]+h*(0-(4*ey[i-1])/x[i-1]-2/(x[i-1]*x[i-1])-ey[i-1]*ey[i-1]);
	y[i]=(1-x[i])/(x[i]*x[i]-x[i]/2);
	wc[i]=y[i]-ey[i];
	}
	printf("/*************************************\\\n");
	printf("x[i]     Eular answer   sure answer      误差\n");
    for(i=0;i<n+1;i++)
    {
		printf("%lf   ",x[i]);
    	printf("%lf      ",ey[i]);
	    printf("%lf    ",y[i]);
	    printf("%lf",wc[i]);
    	printf("\n");
    }
}