:
cat > .$$.c <<'End of Source'
/*
하나씨께서 옳다고 믿어 의심치 않는 계산 법으로 정지길이 Lo인 점이
임의의 속도로 이동할때 관측계에서 측정된 변화된 길이 L과
점의 위치와 이동계의 시간을 계산하는 code입니다.
입력은 광속 c, 길이 Lo, 계산하고자 하는 최대 시간 T Max, 시간간격 dt
그리고 매 시간의 속도 V입니다.
출력은 매 시간의 속도로 부터 sqrt(1-(v/c)**2)의 관계식을 이용하여 계산된
변화된 길이 L, 위치 X 그리고 이동계의 시간 T'입니다.
계산을 좋아하시니 지겹도록 계산해 보시고 논문써서 PRL에도 계제 하세요.
물론 뒷일은 책임 못집니다. 전적으로 하나씨 아이디어니까요.
아침에 짬을 내서 만들었습니다. 하나씨의 학문적 성취에 도움이 되어
뛰어난 업적에 조금이라도 믿거름이 될수 있다면 기쁘기 그지없겠습니다.
버그가 발견될 수 있습니다. 일일이 확인해보지는 않아서.
있으면 고쳐서 사용하세요. 간단한 계산이니까.
*/
#include <math.h>
#define G(a) pow(1.-pow(a/c,2),0.5)
void main() {double c,Lo,tmax,dt,t,x,v,tt; for(printf("c = "),scanf("%lf",&c),
printf("Lo = "),scanf("%lf",&Lo),printf("T Max = "),scanf("%lf",&tmax),
printf("dt = "),scanf("%lf",&dt),x=0.,t=0.,tt=0.;t<=tmax;
printf("T = %E, V = ",t),scanf("%lf",&v),printf("L = %E, X = %E, T' = %E\n\n",
G(v)*Lo,x+G(v)*Lo,tt),x+=v*dt,tt+=dt*G(v),t+=dt);}
End of Source
cc -o hana .$$.c -lm; \rm .$$.c; hana