circuit = Circuit('MOSFET AMPLIFIER')
circuit.model('mynmos','NMOS', level=1)

# Define the DC supply voltage value


circuit.V(1,'Vdd',circuit.gnd, 10@u_V)

input1=circuit.SinusoidalVoltageSource(2, 'in1', circuit.gnd, offset=0@u_V, amplitude=0.1@u_V, frequency=50@u_Hz)
circuit.C(1,'in1','in1c',1@u_uF)
circuit.R(2,'in1c','Vdd',220@u_kOhm)
circuit.R(3,'in1c',circuit.gnd,47@u_kOhm)


circuit.MOSFET(1, 'out', 'in1c', circuit.gnd, circuit.gnd, model='mynmos')
circuit.R(1,'out','Vdd',220@u_kOhm)
circuit.C(2,'out','out2',2200@u_uF)
circuit.R(4,'out2',circuit.gnd,220@u_kOhm)

simulator = circuit.simulator(temperature=25, nominal_temperature=25)
analysis = simulator.transient(step_time=input1.period/5000, end_time=input1.period*10)


#figure, ((ax1, ax2)) = plt.subplots(2, 1, figsize=(8, 5))

figure, ax = plt.subplots(figsize=(20, 10))
ax.grid()
#ax2.grid()
ax.set_xlabel('Time [s]')
ax.set_ylabel('Voltage [V]')

ax.plot(analysis['in1'])
ax.plot(analysis['in1c'])
ax.plot(analysis['out2'])
#ax1.plot(analysis['out2'])

#ax.legend(('in1','in2', 'out'), loc=(.05,.1))
#ax.set_ylim(-1,4)
plt.tight_layout()
plt.show()