69,371
社区成员
发帖
与我相关
我的任务
分享
#include "stdio.h"
#include<math.h>
#define F(y,t)(y+2*exp(4*t))
#define w(t)((2.0/3.0)*exp(4.0*t)-(11.0/3)*exp(t))
void main()
{
int i, n = 10000;
double h = 0.0001, y, t, y0 = -3, t0 = 0;
y = y0;
t = t0;
printf("t y(t) w(t) error\n");
printf("===========================================================================\n");
printf("%.21f %10.61f %10.61f %10.61f\n", t, y, w(t), fabs(y - w(t)));
for (i = 1; i <= n; i++)
{
y = y + h*F(y, t);
t = t + h;
if (i % 1000 == 0)
printf("%.21f %10.61f %10.61f %10.61f\n", t, y, w(t), fabs(y - w(t)));
}
return;
}