数値計算法レポート2-オイラー - 阪大理学部科学計算・開発サークル -Circle for Science Caluculation and Development-

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "report2.h"
 
int main(){
  double x,y,vx,vy;
  int i;
  FILE *fpdata;
  if((fpdata=fopen("euler.dat","w")) == NULL){
    printf("error:cannot open file\n");
    exit(1);
  }
 
  for(i=-10;i<=10;i++){
    x=min;
    y=i*20;
    vx=sqrt(Energy/mass)*c;
    vy=0;
    while(x<=max && y<=max && y>=min && x>=min){
      fprintf(fpdata,"%lf %lf ",x,y);
      x=x+vx*dt;
      y=y+vy*dt;
      vx=vx+Z*197/137/mass*x/(pow(x*x+y*y,1.5))*c*c*dt;
      vy=vy+Z*197/137/mass*y/(pow(x*x+y*y,1.5))*c*c*dt;
    }
  }
  fclose(fpdata);
}