静電ポテンシャルのグラフを計算機で書いてみたい。[訂正]

2日前の話の間違いに気づいた。
$k=\Large \frac{1}{4 \pi \epsilon_0}$
にしなきゃならないところを、
$k=4 \pi \epsilon_0$
にしてた。
以下、修正コード

#include <stdio.h>
#include <math.h>
#define e 1.602*pow(10,-19)	//素電荷
#define ep 8.854*pow(10,-12)	//真空の誘電率
#define k 1/4.*M_PI*ep

float esp(double r);	//静電ポテンシャル,Electro Static Potential

int main(void){
	int i;
	double Us,Ua,Ub,r;
	for(i=-33;i<34;i++){
		r=i*pow(10,-12);
		Ua=esp(r-0.34*pow(10,-10));
		Ub=esp(r+0.34*pow(10,-10));
		Us=Ua+Ub;
		printf("%e %e %e %e\n",r,Ua,Ub,Us);
	}
	return 0;
}
float esp(double r){
	double U;
	U=-k*pow(e,2)/pow(r,2);
	return U;
}

結果
f:id:zf0rest:20171007194815p:plain
前回のと比較してみよう…