42/CH4/EX4.12/sadiku_4_12.sce


1
2
3
4
5
6
7
8
9
10

clear;
clc;
q=10*10^-6;
function[V]=pot(r,P,Q);
V=10*sin(P)*cos(Q)/r^2;
endfunction
Va=pot(1,%pi/6,2*%pi/3);
Vb=pot(4,%pi/2,%pi/3);
W=q*(Vb-Va);
disp(W*10^6,'Work done in uJoule');