function u2 = RCdif (t, u)
u1=220;
RC=10;
u2= (u1 - u)/RC;