地下倾斜电偶极子电流场正演

今天地磁学试验，做了地下倾斜电偶极子的试验，实验要求为

1、推导并计算地下电偶极子的空间分布；

2、熟悉使用相关软件绘出地下不同倾角电偶极子地面电场强度、电位场分；

3、熟悉不同参数偶极子地面场分布规律及特征。

我选了15度的倾角做了一下，公式就不列出来了可以从代码中看出。
#include<stdio.h>
#include<math.h>
#include int main()
{
float p=200,I=1.2,pai=3.141592653,L=200,qingjiao=15,h=300,E;
int i,j,x,y;
FILE*fp;
fp=fopen("xxx.txt","r+");
for(x=-800;x<=800;x=x+10)
{
for(y=-800;y<=800;y=y+10)
{
E=-10000*I*p/(2*pai)*pow(pow((x+L/2*cos(15/180.0*pai)),2.0)+y*y+pow((h-L/2*sin(15/180.0*pai)),2),(-1.5))
+10000*I*p/(2*pai)*pow(pow((x-L/2*cos(15/180.0*pai)),2.0)+y*y+pow((h+L/2*sin(15/180.0*pai)),2),(-1.5)) ;
printf("%d\t%d\t%f\n",x,y,E);
fprintf(fp,"%d\t%d\t%f\n",x,y,E);
}
}
fclose(fp);
return 0;
}

数据还是在文件中，把数据可以用sufer成图，很简单的