Méthode euler circuit RLC
Bonjour à tous, voici mon problème :
je dois résoudre q''+w0/Qq'+w0²q=0 par la méthode d'euler, w0 et Q sont des paramètres connus et les conditions initiales sont telles que q(t=0)=q0 et q'(t=0)=0
lorsque j'execute mon code j'ai l'erreur suivante "TypeError: 'float' object is not callable" à la ligne 20
voici mon code help je dois le rendre demain soir
Code:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
| import matplotlib.pyplot as plt
import numpy as np
import scipy.integrate as spi
import math
R = 100
C = 0.00003
L = 1
w0 = (float(L)/float(C))**0.5
z = (float(L)/float(R))*(L/C)**0.5
q0 = 5
t=0
q=q0
dq=0
temps=[t]
Qt=[q]
h=0.01
while t<14:
q,dq=q+h*dq,dq-h*((-w0/z)(dq)-q*w0**2)
t=t+h
Qt.append(q)
temps.append(t)
plt.plot(temps,Qt,'r') |
