/*
x_{n+1}=y_n+a*x_n+2/(1+xn^2)
y_{n+1}=-x_n
*/
#include
#include
#define a 1.96
main()
{
/*initial condition*/
double x=4.0,y=0.0,nx=0.0,ny=0.0, u=0.0,v=0.0;
int i,n=120000;
for(i=1; i<=n;i++)
{
nx=y+a*x+2.0/(1.0+x*x);
ny=-x;
u=(x+y)/sqrt(2.0-a);
v=(x-y)/sqrt(2.0+a);
if(i>=60000){
printf("%f %f\n",u,v);
}
x=nx;
y=ny;
}
return 0;
}
出力結果
Gif Animation
【How to create a Gif Animation】
Back to C Language
