How do i do this in scilab code?

In the following figure, the switch is closed at $t=0.$ The current, $i,$ flowing through the circuit shown, can be described by a second$-$order differential equation:

$L\frac{d^{2}i}{d{t}^2}+R\frac{di}{dt}+\frac{1}{C}i=V$

which can be rearranged to give:

$\frac{d^{2}i}{d{t}^2}=V-\frac{R}{L}\frac{di}{dt}-\frac{1}{LC}i$

The behavior of this system depends upon the relative values of $L,C,$ and $R($the inductance, capacitance, and resistance$).$ When

$R^2>\frac{4L}{C}$

$R^2<mfrac><mrow><mn>4</mn><mi>L</mi></mrow><mrow><mi>C</mi></mrow></mfrac>$

$R^2=\frac{4L}{C}$

The system is over$-$damped'.

The system is 'over$-$damped'.

The system is 'under$-$damped'.

The system is 'critically damped'.

Use the graphical programming method in Scilab to model the system response, assuming that, V $=5V,R=100,000,$ and $C=1F.$ Assume initial current $i(0)=0A.$

Select values of $L$ to meet each of the damping conditions described above, $(3$ pts $)$

Plot $i$ vs $t$ for all three values of $L.$ There will be three plots altogether. The plots need to be properly labeled. Choose a suitable time scale range and take integration steps to be $0\mathrm{.}001$ s by changing the period of the clock. $(10$ pts $)$

Comment on the differences between the plots in terms of damping. $(7$ pts$)$