dsolve

    1热度

    2回答

    重新安装Python的位置,下面简单的代码 import sympy as sm x = sm.Symbol('x') f = sm.Function('f') y = sm.dsolve(sm.diff(f (x),x)-3*f(x)(1-0.5f(x)),f(x)) print(y) 给出了下面的输出后: Eq(x + 0.333333333333333*log(1.0*f

    0热度

    1回答

    如何绘制以下函数?我需要x的范围[0; 1]。 syms y(x) y(x) = dsolve(diff(y,x) == tan(x), y(0) == 1); plot(y, [0 1]);

    0热度

    1回答

    我用NDSolve对Mathematica进行了一些模拟,我需要介绍温度的影响。我定义了一个随机数表,然后做一个函数出来的,这样说: randomtablex = Table[RandomVariate[NormalDistribution[]], {i, 1, IntegerPart[3 tspacer/deltats] + 1}]; randomtabley

    0热度

    1回答

    我一直在玩这一段时间,任何人都可以解释为什么我从Code1和Code2得到不同的答案?关于'dsolve()'的实际脚本是什么使得2个代码的输出不同?可以输出,如果我只是用不同的语法相同(即”;”。“) %Code1: syms Qua t Area height rate_in = 3*Qua*(sin(t))^2; delta_Vol = dsolve('DAreaY = rate_

    1热度

    1回答

    我试图使用以下代码: ode1='D2y1=-1256.4*y1-5*Dy1+255.1*y2+182.781'; ode2='D2y2=-151.5*y2-5*Dy2+255.1*y1-14.0459'; CI='y1(0)=2,y2(0)=-2,Dy1(0)=0,Dy2(0)=0'; sol=dsolve(ode1,ode2,CI,'t'); sol.y1 sol.y2 和MA

    3热度

    1回答

    我有一个语法问题在Mathematica(第10版)中求解一个微分方程。 的公式,我需要解决的输入如下: solv = DSolve[{ a*u''[y] - b*u[y] == d, u'[0] == 0, u[1] == 0}, u, {y, -1, 1}] 其中使用ExpToTrig和FullSimplify后,我得到我要的答案: (d (-1 + Cosh[(Sqrt[b] y)/Sq

    0热度

    1回答

    我有我的程序解微分方程: p := dsolve({ic, sys}, numeric, method = rosenbrock); 解决,我们有以下后: print(p(10)); [t = 10., alpha(t) = HFloat(0.031724302221312055), beta(t) = HFloat(0.00223975915581258)] 我需要用这个α(

    0热度

    1回答

    基本目标:其中 a1=k*H*H/c; a0=-H*H*tau0/c; c=0.1; k=1; H=1; tau0=1; 我使用dsolve解决 & 情节的ODE D2f=a1*f+a0 的解决方案命令找到具有上述常量值的解决方案。 sol=dsolve('D2th=a1*th+a0','th(0)=0','th(1)=0','t'); 非但没有与常量的分配值的解决方案,我得到的

    -1热度

    1回答

    我想解微分方程。 MATLAB显示警告: clear all syms x f(x) theta eq = (-6*x+(-7+theta)*f*diff(f,x))*(1+diff(f,x)^2)+x*f*(diff(f,x,x))==0 cond = f(0)==1 dsolve(eq,cond) Warning: Explicit solution could not be fo

    5热度

    2回答

    我使用SymPy版本0.7.3,并在使用dsolve函数时遇到一些问题。似乎dsolve在输入方程有太多参数时有困难。 我试过的公式如下: from sympy import * p = Function('p') t, u1, u2, u3, u4, u5 = symbols('t u1 u2 u3 u4 u5') eq = Eq(Derivative(p(t),t), -(u3 + u4