Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #include <stdio.h>
- #include <math.h>
- #define nx 100
- #define ny 100
- #define x0 0.0
- #define dx 0.1
- #define y0 0.0
- #define dy 0.1
- #define q 1.0
- #define xp 5.01
- #define yp 5.01
- // x0 e y0 é posição do canto da matriz
- // xp e yp é posiçao da particula
- int main(void)
- {
- double M[nx][ny],x,y,V;
- int i,j;
- M[i][j]=0.0;
- x=x0;
- y=y0;
- for (i=0;i<nx;i++)
- {
- for(j=0;j<ny;j++)
- {
- x=x0+dx*i;
- y=y0+dy*j;
- V=q/(sqrt(pow((x-xp),2)+pow((y-yp),2)));
- M[i][j]=V;
- printf("%lf \t %lf \t %lf \n",x,y,V);
- }
- }
- return 0;
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement