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

勉強してて、
f:id:zf0rest:20171005185338j:plain
こんな綺麗なグラフ作ってみたいなぁと思い、やっぱり計算機でやってみた。

#include <stdio.h>
#include <math.h>
#define e 1.602*pow(10,-19)	//素電荷
#define ep 8.854*pow(10,-12)	//真空の誘電率
#define k 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:20171005190353p:plain
あれ、なんか思ってたのと違う。
これで合ってんのかも、怪しい…

まぁいいや。思い付きだし。