如何将以下用maple编写的代码转换为Python框架
ode:=diff(y(x),x)= #=An equation
sol:=dsolve([ode, y(1)=0.73}, numeric, output=listprocedure)
fy:=eval(y(x),sol);nn:=100
for x from 1 to nn do fy(5.1*x/nn) end do;
我找不到类似于上述程序的代码。谢谢你的帮助
我正在尝试使用SymPy将一些枫树代码转换为python。下面是我要转换的枫树代码:
下面是我尝试用python代码来做同样的事情,但这是我所得到的。
from sympy import *
c, k, x, t, tau, Y, xi = symbols('c k x t tau Y xi')
u = Function('u')
F = Function('F')
xi = k*(x-c*t)
U = Function('U')
pde = diff(u(x,t), t) + 6*u(x,t)*diff(u(x,t)
我对python完全陌生,并尝试集成以下ode:
$\dot{x} = -2x-y^2$
$\dot{y} = -y-x^2
这会导致数组中的所有内容都为0,但我做错了什么呢?它主要是复制代码,与另一个,而不是耦合的ode,它工作得很好。
import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import ode
def fun(t, z):
"""
Right hand side of the differential equations
dx/d
我有一个参数化的小图,我正试图从Mathematica转换成python。问题是我的python脚本与此不匹配,我正在努力使列表图工作,因为它将在Mathematica中使用。
您将如何使用matplotlib在python中编写这篇文章?
With[{r = 2.1, dt = 0.5, tfinal = 7.0},
x = 1;
y = 0;
xylist = {{x, y}};
Do[
x = x - r y dt;
y = y + r x dt;
AppendTo[xylist, {x, y}];
, {i, 0, tfinal, dt}];
p1 = Par
我在python的complex_ode解决程序上遇到了麻烦。
我试图解决以下问题:
dy/dt = -iAy - icos(Omegat)_B_y
其中A和B是NxN数组,未知y是Nx1数组,I是虚单位,欧米茄是参数。
这是我的密码:
import numpy as np
from scipy.integrate import ode,complex_ode
N = 3 #linear matrix dim
Omega = 1.0 #parameter
# define symmetric matrices A and B
A = np.random.ranf((N,N))
A =
我在试着解决一个汽车问题。首先,我有一个汽车问题的原始代码:
# Ampl Car Example
#
# Shows how to convert a minimize final time optimal control problem
# to a format pyomo.dae can handle by removing the time scaling from
# the ContinuousSet.
#
# min tf
# dxdt = 0
# dvdt = a-R*v^2
# x(0)=0; x(tf)=L
# v(0)=0; v(tf)=0
# -3<=a<
我正在学习sympy,并希望验证ODE的解决方案。我还不太明白同情的命名规则。
我不想执行在顶部加载所有包的标准方法,我只想使用import sympy,然后使用显式的长名称来引用所有其他名称。关于最新的conda python
Python 3.7.3 (default, Mar 27 2019, 22:11:17)
[GCC 7.3.0] :: Anaconda, Inc. on linux
打字时
import sympy
x = sympy.symbols('x')
y = sympy.Function('y')
ode = sympy.Eq(
我对ode45有个问题。我已经定义了一个函数,并试图用ode来解决它,但当我运行它时,它花费了这么长的时间。我试着在我的函数中显示"t“输入,它显示时间步长是10^-8!我没有从ode45得到任何错误
所以我在我的函数的末尾放了一个断点,在我一步之后,它转到sym.m文件并调用函数delet(h)
function dxr=Dynfun(t,x)
...
dxr=[A;B]
在步骤之后,它将转到
function delete(h)
if builtin('numel',h)==1 && inmem('-isloaded',
我正试图解决以下几个问题:
function [eta, sol] = compressible_similarity_wo
global Gamm Ma Pr omega;
Gamm = 1.4;
Ma = 2;
Pr = 0.7;
omega=0.76;
global eta_max_ode;
eta_max_ode = 20;
opt = optimset('Display','off','TolFun',1E-20);
F = fsolve(@(F) eva
根据我的语言知识,我的代码是正确的。但这并没有给我正确的解决方案(情节)。当我用mathematica解决了同一个ODE系统时,我得到了正确的解,两种解都是完全不同的。我正在编写一个研究项目,所以我需要在python中编写适当的代码。你能告诉我我密码的错误吗。
import numpy as np
import matplotlib.pyplot as plt
import scipy.integrate as si
##Three system
def func(state, T):
H = state[0]
P = state[1]
R = state[2]
我有一个代码(需要用户给出大量输入),它将给我一个n×n矩阵(例如A),我不得不用它来求解一个ODEs X'=AX系统。如何将此矩阵A包含在ode45的函数文件(ode45文件)中。如果我将代码包括在函数文件中,如下所示:
function xp=g_test(t,x);
k=input('something');
A=some manipulation of inputs;
xp=A*x;
end
Matlab要求在每个时间步骤输入(通常情况下,我的问题有30k的时间步骤)。那么,如何将矩阵A包含/传递给函数呢?
我正在为一个需要的项目从头开始编写我的第一个cmake文件。如果我正在做的其他坏习惯,对此的评论也很受欢迎。
使用cmake的makefile生成步骤正在工作,但是在使用make的构建步骤中,它会抱怨:
Scanning dependencies of target ode_oscillex
[ 50%] Building CXX object CMakeFiles/ode_oscillex.dir/ode_oscillex.cpp.o
[100%] Linking CXX executable ode_oscillex
ld: library not found for -lGLU
clan