Contact Us!
CoCalc Logo Icon
StoreFeaturesDocsShareSupportNewsAboutSign UpSign In

Real-time collaboration for Jupyter Notebooks, Linux Terminals, LaTeX, VS Code, R IDE, and more,
all in one place. Commercial Alternative to JupyterHub.

| Download
Views: 22
Image: ubuntu2004
Tweet
Share
Share
x,t=var('x, t')
T=function('T')(x,t)
k_per = 1.7
rho_per=1
Cp_per =4186
#Définir l'équation de la chaleur
equ_chal = diff(T,t)*rho_per*Cp_per - diff(T,x,2)*k_per==0
#recherher les solutions analystiques de l'équation différencielle
sol = desolve(equ_chal, [x,t]); sol
#conditions aux limites 
sol0=desolve(equ_chal,[x,t],ics=[0,0]);print(sol0)
sol1=desolve(equ_chal,[x,t],ics=[0,1]);print(sol1)
sol2=desolve(equ_chal,[x,t],ics=[0,2]);print(sol2)
# tracer
sol=desolve(equ_chal,[x,t],ics=[0,0]); sol
plot(sol.substitute(r==0.5) ,(0,0))
Error in lines 7-7 Traceback (most recent call last): File "/cocalc/lib/python3.10/site-packages/smc_sagews/sage_server.py", line 1244, in execute exec( File "", line 1, in <module> File "/ext/sage/9.7/src/sage/calculus/desolvers.py", line 598, in desolve raise NotImplementedError("Maxima was unable to solve this ODE. Consider to set option contrib_ode to True.") NotImplementedError: Maxima was unable to solve this ODE. Consider to set option contrib_ode to True. 


x,t=var('x, t')
T=function('T')(t)
TT=function('T')(x)
k_per = 1.7
rho_per=1
Cp_per =4186
desolve(diff(T,t)*rho_per*Cp_per==0,t)
desolve(diff(TT,x,2)*k_per==0,x)
Error in lines 7-7 Traceback (most recent call last): File "/cocalc/lib/python3.10/site-packages/smc_sagews/sage_server.py", line 1244, in execute exec( File "", line 1, in <module> File "/ext/sage/9.7/src/sage/calculus/functional.py", line 150, in derivative f = SR(f) File "sage/structure/parent.pyx", line 897, in sage.structure.parent.Parent.__call__ return mor._call_(x) File "sage/structure/coerce_maps.pyx", line 161, in sage.structure.coerce_maps.DefaultConvertMap_unique._call_ raise File "sage/structure/coerce_maps.pyx", line 156, in sage.structure.coerce_maps.DefaultConvertMap_unique._call_ return C._element_constructor(x) File "sage/symbolic/ring.pyx", line 384, in sage.symbolic.ring.SymbolicRing._element_constructor_ return new_Expression(self, x) File "sage/symbolic/expression.pyx", line 13749, in sage.symbolic.expression.new_Expression raise TypeError(f"unable to convert {x!r} to a symbolic expression") TypeError: unable to convert T to a symbolic expression