%PROGRAM TO FIND OUT ELECTRIC FLUX DENSITY
%IN A COAXIAL CYLINDER
%BOTH INNER CYLINDER & OUTER CYLINDER

clear all;
clc;
  
a = input('Enter the radius of inner cylinder in meters = ');

b = input('Enter the radius of outer cylinder in meters = ');

Psinner = input('Enter the value of surface charge density in coulomb/square meter = ');

P = input('Enter the value of radial distance in meters = ');

if P<a
               
        disp('Electric Flux Density within inner conductor  = 0');
    
    elseif P>b
                 
        disp('Electric Flux Density outside the outer condcutor = 0');
        
        
    else       
          
        disp('Electric Flux Density between inner condcutor & outer conductor =');
        Dp = (Psinner*a)/P
        
        disp('Surface charge density in outer condcutor = ');
        Psouter = -Psinner*a/b
    
end