縮退宇宙論Part2

1 ：NAS6 ◆n3AmnVhjwc ：2014/01/21(火) 22:55:23.33 ID:1azff/Tu.net

ttp://www5b.biglobe.ne.jp/~NAS6/secret/BlackHole.htm

シュワルツシルト時空を解いた

要旨抜粋
u=1/r
シュワルツシルト補正項3.0mu^2
u''+u=m/h^2+3.0mu^2
ケプラーの法則
u''+u=kM/H^2=GM/r^2=g
を参照
m/h^2でくくり
S=3.0mu^2/(m/h^2)=3.0(1/r^2)r^4ψ'^2=3.0r^2(dψ/dt)^2(1/c^2)
r(dψ/dt)/c=v、(v=rω)、(v=V/c)なので
S=3.0v^2
3.0mu^2をm/h^2でくくったので
F=mg(1.0+S)
g=GM/r^2

よってシュワルツシルト時空の重力による加速度は
a=(GM/r^2)(1+3(V/c)^2)
となりました

参照リンク
ttp://members3.jcom.home.ne.jp/nososnd/grel/peri.pdf

41 ：NAS6 ◆n3AmnVhjwc ：2014/01/23(木) 10:02:50.36 ID:T8E0hyf4.net

√3　log↓10　N
あたりが単純に正確に近似できるんじゃないの？

42 ：NAS6 ◆n3AmnVhjwc ：2014/01/23(木) 10:24:32.84 ID:T8E0hyf4.net

ttp://globe.asahi.com/feature/100201/04_2.html

リーマン予想

43 ：NAS6 ◆n3AmnVhjwc ：2014/01/23(木) 10:36:22.48 ID:T8E0hyf4.net

円周率テーラー展開
Σ↑∞↓n=0　((-1)^n/2ｎ+1)

(-1)^nはi^2nってことだ

44 ：NAS6 ◆n3AmnVhjwc ：2014/01/23(木) 10:52:44.29 ID:T8E0hyf4.net

ここまでのまとめ

π
=lim(n→∞)nsin(π/n)
=lim(n→∞)Σ↑∞↓n=1　((-1)^n/2ｎ+1)
=lim(n→∞)√(6　Σ↑∞↓n=1　(1/n^2))　
=lim(n→∞)√(6　Π↓p　(1/(1-1/p^2)))　

45 ：NAS6 ◆n3AmnVhjwc ：2014/01/23(木) 11:17:02.32 ID:T8E0hyf4.net

でようするに
π
=lim(n→∞)nsin(π/n)
なんだから

sinテーブル(波)を考えつつ
sin(π/n)×ｎ=π
ってことは

何がいいたいかは分かるだろ
0<x<π/4に注目して

半角公式
sin(θ/2)=√((1-cosθ)/2)
cos(θ/2)=√((1+cosθ)/2)

にほかならない

46 ：NAS6 ◆n3AmnVhjwc ：2014/01/23(木) 11:31:12.70 ID:T8E0hyf4.net

だから
cosπ/2=0→cos0=1
のcosテーブルは
cos(θ/2)=√((1+cosθ)/2)
で描けることになる
この方法で細かい(半角をたくさん計算した)とこまで計算すれば
荒いとこは細かいとこの鏡写しだよ

47 ：ご冗談でしょう？名無しさん：2014/01/23(木) 13:17:31.50 ID:Z4CNjSGc.net

2n=2log↓10(N)
数学的証明では合っているんだけど
誤りがあったら指摘してくれ

48 ：NAS6 ◆n3AmnVhjwc ：2014/01/23(木) 13:23:43.13 ID:T8E0hyf4.net

#include <stdio.h>
#include <stdlib.h>
#define _USE_MATH_DEFINES
#include <math.h>
#include <time.h>

int factorial(int n){
　int ret=n,i=n;
　if(n=0)return 1;
　for(i-=1;1<i;i-=1)ret=ret*i;
　return ret;
}
int Binomial(int n,int k){
　int a=factorial(n);
　int b=factorial(n-k);
　int c=factorial(k);
　if(b*c==0)return 1;
　return (a/(b*c));
}
double halftheta(double theta){
　return sqrt((1.0+theta)/2.0);
}
double ntheta(double theta,int n){
　double ret=0.0;
　int i,j;
　for(i=0;i<n/2;i++){
　　for(j=i;j<n/2;j++){
　　　ret+=pow(-1.0,(double)i)*(double)Binomial(n,2*j)*(double)Binomial(j,i)*pow(theta,(double)(n-2*i));
　　}
　}
　return ret;
}

49 ：NAS6 ◆n3AmnVhjwc ：2014/01/23(木) 13:24:36.73 ID:T8E0hyf4.net

void main(){
　double a,b,c,d;
　int i,MAX_N=30;
　FILE* fp;
　fp=fopen("data.txt","w");
　for(a=0.5,b=0.0,i=0;i<MAX_N;i++){
　　fprintf(fp,"%.15lf,%.15lf\n",a,b);
　　//if(i==26){c=a;d=b;}
　　if(i==4){c=a;d=b;}
　　a*=0.5;
　　b=halftheta(b);
　}
　fprintf(fp,"######################################\n");
　a=0.0;b=1.0;
　fprintf(fp,"%3d %.15lf,%.15lf\n",0,a,b);
　for(a=c,b=d,i=2;a<=1.0;a+=c,i++){
　　fprintf(fp,"%3d %.15lf,%.15lf\n",i-1,a,b);
　　b=ntheta(d,i);
　}

　fclose(fp);
}

どこがばぐってんだ？13倍角くらいからおかしい

50 ：NAS6 ◆n3AmnVhjwc ：2014/01/23(木) 13:31:25.59 ID:T8E0hyf4.net

　if(i==4){c=a;d=b;}
を
　if(i==6){c=a;d=b;}
　
にしても13番目の使徒に裏切られる