class Cap():
    def __init__(self,res1,cap1,source):
        self.res1 = res1
#        self.res2 = res2
        self.cap1 = cap1
#        self.cap2 = cap2
        self.vsrc = source
    def charge(self,vc_init):
        t=np.linspace(0.,1,1000,endpoint=False)
        vc = np.zeros(len(t))
        if vc_init is None:
            vc[1]=0
        else:
            vc[1]=vc_init
        for i in range(len(t)):
            vc_ = vc[i]
            ir1= (self.vsrc-vc_)/self.res1
            vc[i+1]=vc[i] + (ir1*i)/self.cap1
#            if sinput <=treshhold:
#                break
        print 'vc=', vc
        print len(vc)
        print  len(t)
#        print 'z=',z
        return vc