2つ目のエラーは配列の添字が範囲外になったときでますね。π二乗は重力定数(9.80665)とだいたい同じで、それを天井関数で切り上げたので10になったのでしょう。 3つ目は半径0を指定したためなので、初期値を1にすると改善されます。