07/06/27 09:15:34
>>344
y' = z とう条件は不要なんじゃないかと思うが
#include <stdio.h>
#include <math.h>
#define N 100
double f(double y, double z)
{
return -z - 2.0 * y + cos(y);
}
int main(void)
{
double y[N], z[N], dx = 0.2;
int i;
y[0] = z[0] = 0.0;
for (i = 0; i < N; i++) {
y[i] = y[0] + i * dx;
z[i + 1] = z[i] + dx * f(y[i], z[i]);
printf("y = %f, z = %f\n", y[i], z[i]);
}
return 0;
}