<jupyter_start><jupyter_text>Como as condições são mistas é necessário usar a série de senos de termos imparesAlém disso, como o lado não isolado não vale 0 existe uma parte permanente na solução, que simplesmente vale T1Assim temos que calcular os termos da série de senos de termos impares com f(x) - T1, com f(x) valendo T1 de 0 a 2L/3 e 0 de 2L/3 a LPara calcular a série de senos impares temos de substituir L por 2L e ajustar os c e d. Além de substituir os ns por 2n+1<jupyter_code>T1 = 10
L = 100
#f(x)
def calculacs1(ns):
    ks = (2*ns + 1)
    c = 0
    d = 1/3
    cn = (-4*(T1)/(np.pi*ks))*(np.cos(np.pi*ks*d)-np.cos(np.pi*ks*c))
    return cn
#Parte permanente
def calculacs2(ns):
    ks = (2*ns + 1)
    c = 0
    d = 1/2
    cn = (-4*(T1)/(np.pi*ks))*(np.cos(np.pi*ks*d)-np.cos(np.pi*ks*c))
    return cn
def soluchom(cs,ns,x,t):
    ks = (2*ns + 1)
    sine = np.sin(ks*np.pi*x/(2*L))*np.exp(-(ks**2*np.pi**2/(2*L)**2)*t)
    val = np.sum((sine*cs)) # parte transitóriarmanente
    val = val + T1 #parte permanente
    return val
n = 100
ns = np.arange(0,n)
ts = np.geomspace(1,20000, num= 9)
ts[0] = 0
xs = np.linspace(0,L,num=100)
fig, axes = plt.subplots(3,3, figsize = (18,9))
k = 0
for i in range(0,3):
    for j in range(0,3):
        axes[i,j].plot(xs,[soluchom(calculacs1(ns)-calculacs2(ns),ns,x,ts[k]) for x in xs])
        axes[i,j].set_title(repr(ts[k])[0:5])
        axes[i,j].set_ylim(0,12)
        k +=1<jupyter_output><empty_output>
{'fraction_non_alphanumeric': 0.16734972677595628, 'fraction_numerical': 0.05327868852459016, 'mean_word_length': 3.3601190476190474, 'pattern_counts': {'":': 0, '<': 5, '<?xml version=': 0, '>': 5, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'chain_length': '1', 'dir': 'jupyter-structured-clean-dedup', 'id': '289762', 'license': 'MIT', 'n_tokens_mistral': 708, 'n_tokens_neox': 666, 'n_words': 174, 'path': 'aula9/aula9.ipynb', 'repo_name': 'gbaraldi/MAP2320-PDEs', 'size': '1464'}