blob: 32b9bb618d8be352a285017e7c9fa2f48719e52a (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
//clear//
//Caption: Program to Determine the electric field of a two infinite radial planes with an interior angle alpha
//Example 7.3
//page 180
clc;
phi = sym('phi');
A = sym('A');
B = sym('B');
Vo = sym('Vo');
alpha = sym('alpha');
aphi = sym('aphi');
r = sym('r');
V = integ(A,phi)+B;
disp(V,'V =');
V = limit(V,B,0);
V = limit(V,A,Vo/alpha);
disp(V,'Potential V after applying boundary conditions =')
E = -(1/r)*diff(V,phi)*aphi;
disp(E,'E =')
//Result
// V = B+phi*A
// Potential V after applying boundary conditions = phi*Vo/alpha
// E = -aphi*Vo/(alpha*r)
|