1 files changed, 25 insertions, 18 deletions
diff --git a/Part 2/RLC_script.m b/Part 2/RLC_script.m
index 73977bc..881d941 100644
--- a/Part 2/RLC_script.m
+++ b/Part 2/RLC_script.m
@@ -1,38 +1,45 @@
-function RLC_script(f,h)
+function RLC_script (tf) 
+
 %initailise the circuits 
-R = 280;
-L = 600*10^(-3);
-C = 3*60^(-6);
-tf = 1/f;
-N = round(tf/h,0);
+R = 250;
+L = 650*10^(-3);
+C = 3*10^(-6);
+h = 0.00001; %step size
+
+%initailise the container 
 
-qc_dash = zeros(1,N);
-qc = zeros(1,N);
-t = zeros(1,N);
-Vout = zeros(1,N);
+N = round(tf/h); %number of iterations
+qc = zeros(1, N); %x
+qc_dash = zeros(1, N); %y  
+t = zeros(1, N);
+Vout = zeros(1, N);
 
 %input voltage
-Vin = @(t)5;
+% step function of 5 volt
+Vin = @(t)5*heaviside(t); 
+%Vin = @(t)
 
 %the coupled equation
-func1 = @(t,qc,qc_dash)qc_dash;
-func2 = @(t,qc,qc_dash)(Vin(t) - qc/C - R*qc_dash) / L;
+func1 = @(t, qc, qc_dash)qc_dash; 
+func2 = @(t, qc, qc_dash)(Vin(t) - qc/C - R*qc_dash)/L; 
 
 %the initial condition
+qc_dash(1) = 0; 
+qc(1) = 500*(10.^-9); 
 t(1) = 0;
-qc(1) = 500*10^(-9);
-qc_dash(1) = 0;
 
 %Rouge Kutta
-for i = 1 : N -1
-    [qc(i+1),qc_dash(i+1)] = RK4second(t(i),qc(i),qc_dash(i),h,func1,func2);
+for i = 1 : N - 1
     t(i + 1) = t(i) + h;
+    [qc(i + 1), qc_dash(i + 1)] = RK4second(t(i), qc(i), qc_dash(i), h, func1, func2); 
     Vout(i) = R*qc_dash(i);
 end
 
 %plot
+
 plot(Vout);
 xlabel('Time');
 ylabel('Amplitude');
-
 end
+
+