ラザフォードの散乱実験のソースコードを書いたが、関数に関するコンパイルエラーが出ました。

前提・実現したいこと

ラザフォードの散乱実験をテーマにプログラミングを書いています。
コンパイルしたところ、16、36行目の”｛”にエラーが出ました。
プログラミングについては初心者ですが、アドバイスよろしくお願いします。

発生している問題・エラーメッセージ

Main.c:16:1: error: function definition is not allowed here

該当のソースコード

C言語

#include <stdio.h>
#include<math.h>

#define N 10
int main(void)
{ 
  double x=-2000, b=2000,y[N]={10,100,150,200,300};
  double n,h,t,vx,vy,E=5.30,ma=6.664*pow(10,-27),M;
  
  M=2*E/ma;
  vx=pow(M,0.5);
  n=40000;
  h=(b-x)/n;

int main()
{
 for(i = 0 ; i <=n ; i++);
  {
     double ki,k2,k3,k4;

        i+=1;
        t+=h;
        k1 = g(t,vy);
        k2 = g(t+h*0.5,vy+k1*0.5);
        k3 = g(t+h*0.5,vy+k2*0.5);
        k4 = g(t+h, y+k3);
        vy = vy + (k1+2.0*k2+2.0*k3+k4)/6.0;
        y=vy*t:
        x=vx*t;
        printf("%2f %lf %lf\n",t,x,y);

    return 0;
   }
}  
int main(void)
  {
  double g(double t, double vy);
  int z=79,m=3.141592,e=8.854*pow(10.-12),r=sqrt(pow(x,2)+pow(y,2));
  
  g=2*z*z*y/(4*m*e*ma*pow(r,3);
     return(g)
  }
}

試したこと

”｛”のインデントを変えたり、指定する関数を変えたりしましたがあまり成果はなかったです。

補足情報（FW/ツールのバージョンなど）

ここにより詳細な情報を記載してください。

行動規範の内容に同意します

回答2件

ベストアンサー

main関数が複数あることも問題ですが、中途半端なインデントは誤読の元です。
きれいに書く習慣をつけましょう。

インデントしていて気づきましたが、２つ目のmain関数内のfor文は自己完結しますよ。

K&R風インデント例

C
1#include <stdio.h>
2#include<math.h>
3
4#define N 10
5int main(void) {
6    double x=-2000, b=2000,y[N]={10,100,150,200,300};
7    double n,h,t,vx,vy,E=5.30,ma=6.664*pow(10,-27),M;
8
9    M=2*E/ma;
10    vx=pow(M,0.5);
11    n=40000;
12    h=(b-x)/n;
13
14    int main() {
15        for(i = 0 ; i <=n ; i++);   /* ←セミコロンがついてるのでループは自己完結しますよ */
16        {
17            double ki,k2,k3,k4;
18
19            i+=1;
20            t+=h;
21            k1 = g(t,vy);
22            k2 = g(t+h*0.5,vy+k1*0.5);
23            k3 = g(t+h*0.5,vy+k2*0.5);
24            k4 = g(t+h, y+k3);
25            vy = vy + (k1+2.0*k2+2.0*k3+k4)/6.0;
26            y=vy*t:
27            x=vx*t;
28            printf("%2f %lf %lf\n",t,x,y);
29
30            return 0;
31        }
32    }  
33    int main(void) {
34        double g(double t, double vy);
35        int z=79,m=3.141592,e=8.854*pow(10.-12),r=sqrt(pow(x,2)+pow(y,2));
36
37        g=2*z*z*y/(4*m*e*ma*pow(r,3);
38        return(g)
39    }
40}
41