问题如图 下面是我计算的代码： DiffusionEquation={D[\[Rho][r,t],t]-d*1/r*D[D[\[Rho][r,t],r]*r,r]==0,\[Rho][r,0]==c,Piecewise[{{D[\[Rho][r,t],r]==c,r<R0},{D[\[Rho][r,t],r]==0,r>R0}}]}; DSolve[DiffusionEquation,\[Rho][r,t],{r,0,+\[Infinity]},{t,0,+\[Infinity]}] 这样的边界条件应该怎么处理？求助各位！

症状： 顶楼中的 [ 和 ] 会被自动转换为 [lbk] 和 [rbk]。此现象在回复帖中不存在。 影响范围： 此bug至少在去年10月16日即已出现，至今（安卓12.58.1.0）未修。IPhone受影响情况不明，欢迎实验补充。 对策： 对于提问者，最简单的对策就是改用电脑访问贴吧网页版，不要使用使用手机客户端发帖。如果确有相关需求（所以说你们到底为什么会有这种需求？电脑没连网？那也可以开手机热点啊，访问网页又花不了什么流量），避免在顶楼张贴代码也是一个

新手绘图时遇到个不知如何下手的问题，理论上曲线中间还有一小段圈将左右两段勾连在一起，但绘图出来后总是缺这么一小段关键部分断开，求助各位吧友，实在不知道问题出在哪里了不清楚是参数太小了精细度不足还是什么其他原因，以下是代码： b0 = 0.5; b1 = 4; b2 = 0.05; b3 = 0.003; x0 = 0.11; sigma = 0.002; k = 0.2; gamma = 0.2; lambda = 0.3; alpha = 0.1; beta[lbk]x_[rbk] := (b0 + b1*x/(b2 + x))*(1 - b3*E^(-(x - x0)^2/sigma)); bb[lbk]x_[rbk] := alpha*k + beta[lbk]x[rbk]*(x + k + gamma); cc[lbk]x_[rbk] :

这个关于编译（Compile）的帖子我也算构思了满久了，不过还是不敢说已经定了稿，为了便于修订并避免链接的迅速报废，这次尝试着放在了有道云笔记上，大家多点一下鼠标吧： http://note.youdao.com/share/?id=bf2a031f7d0d69064941032c7857dfd

代码如下： Plot[{(3600.` 2.5`[250 - 4 (10 + e)])/( 1800 - 22.5` (1 + e)^2 (5 + 3 e) - 0.75` 4[-e^2 + (2.5` + e)^2][4800 - 9.` (1 + e)^2])}, {e, 0, 7}, PlotLegends -> "Expressions", AxesLabel -> {"e", "p1"}, LabelStyle -> Directive[15, FontFamily -> "Times New Roman", FontColor -> Black], ImageSize -> {300}]

In[72]:= f = 5 - r^2 - r; ff[r_] := 5 - r^2; cf1 = Compile[{{r, _Real}}, f] cf2 = Compile[{{r, _Real}}, ff[r]] cf1[2] cf2[2] Out[74]= CompiledFunction[{11, 14., 5470}, { Blank[Real]}, {{3, 0, 0}, {3, 0, 1}}, {{5, {2, 0, 0}}}, {0, 1, 4, 0, 0}, {{46, Function[{r}, f], 3, 0, 0, 3, 0, 1}, {1}}, Function[{r}, f], Evaluate] Out[75]= CompiledFunction[{11, 14., 5470}, { Blank[Real]}, {{3, 0, 0}, {3, 0, 1}}, {}, {0, 0, 2, 0, 0}, {{47, ff, 3, 0, 0, 3, 0, 1}, {1}}, Function[{r}, ff[r]], Evaluate] \:6B63\:5728\:8BA1\:7B97In[72]:= CompiledFunction::cfse: 被编译表达式 5-r-r^2 应该是一个 machine-s

贴吧新弄了个贴条，贴条只能同时存在10个并且只能实际只能显示26字，于是我就在想着是不是该总结些Mathematica新手最常见的问题贴到那上面滚动播出（虽然很多人或许不会看，但有总比没有好嘛……）。试着弄了个草稿（二十戒），大家可以看看，哪些需要修改，哪些需要删节或补充： 二十戒 1　不要为了节约硬盘选择过老的版本。至少要用版本7 2 　自带帮助是最好的教材，提问前先打开软件按下F1仔细查查 3　将光标停在不认识的函数前/中/后再

我现在的动画干涉过程和效果都不太明显，想问下可以怎么去调节呢？ Clear[lbk]"Global`*"[rbk]A1 = 1; w1 = 1; \[lbk]CurlyPhi[rbk]1 = 0; k1 = 1; k2 = 1;Manipulate[lbk]Plot3D[lbk]A1*Cos[lbk]w1*t + \[lbk]CurlyPhi[rbk]1 - k1*Sqrt[lbk]x^2 + y^2[rbk][rbk] + A2*Cos[lbk]w2*t + \[lbk]CurlyPhi[rbk]2 - k2*Sqrt[lbk]x^2 + (y - 5)^2[rbk][rbk], {x, -20, 20}, {y, -20, 20}, AxesLabel -> {"x", "y", "wave1+wave2"}, Mesh -> None, ColorFunction -> "Rainbow", ColorFunctionScaling -> False, PlotPoints -> 150[rbk], {t, 0, 10}, {

也算是对这帖（http://tieba.baidu.com/p/2964416898）的2，3，4条的解说吧。有什么意见或建议欢迎提出： http://note.youdao.com/share/?id=abd51087f44c0b6a41ff6022d549dc4

FullSimplify[Solve[\[Pi]mYR - \[Pi]mNR == 0, \[Lambda]]] {{\[Lambda] -> Root[-17916796875000000000 + 31250000000000000000 \[Omega] - 136643898437500000000 #1 + (-218550898945312500000 + 2187500000000000000 \[Omega]) #1^2 - 14729351792968750000 #1^3 + (121134832362753906250 - 22050000000000000000 \[Omega]) #1^4 + 351782505267375000000 #1^5 + (348916432147793125000 - 12691000000000000000 \[Omega]) #1^6 - 125878931825762500000 #1^7 + (-85747006371034312500 - 2689120000000000000 \[Omega]) #1^8 - 5676317980870600000 #1^9 + (-4448713765162793000 - 201684000000000000 \[Omega]) #1^10 + 130974452817

我运行之后他一直显示正在运行，但是不出图，想问一下代码哪里错了吗 α = 0.6; γ = 0.5; t[lbk]lbk[rbk]x_[lbk]rbk[rbk] := (1/π) Integrate[lbk]lbk[rbk](-1)^(Gamma[lbk]lbk[rbk]1 + α + 0.5 t[lbk]rbk[rbk])/(Gamma[lbk]lbk[rbk]1 + t[lbk]rbk[rbk]) Cos[lbk]lbk[rbk](1 + α + 0.5 t) π/2[lbk]rbk[rbk] x^(1 + α + 0.5 t), {t, 0, ∞}[lbk]rbk[rbk] Plot[lbk]lbk[rbk]t[lbk]lbk[rbk]x[lbk]rbk[rbk], {x, 0, 1}[lbk]rbk[rbk] 原公式是图二

吞私信标准并不明确。 网址，包括不限于qq号手机号的各类号码，过长的mathematica代码，似乎都有可能被吞。 如果确有发私信的需求，建议在发完之后补发一条“我刚发了X条私信，你是否收到了。” 另，自从百度贴吧封杀pc端私信（是的，封杀了，尽管你在pc端看到的提示是“对方不接受私信”）后我个人就再也没回答过私信中的问题；我个人也一直不喜欢非公开的mathematica讨论，希望大家尽量在吧内提问以便知识的继承。

如果要将动态的图像放在ppt里，要怎么做，需要保留互动按键

各位大佬，，大家知道如何把带控件的交互式动态图放到PPT中吗？

Clear["Global`*"] mu = SetPrecision[5/1000, 20]; \[CapitalLambda] = SetPrecision[65/100, 20]; Gs = SetPrecision[493/100, 20]; r = SetPrecision[1/10, 20]; a0 = SetPrecision[675/100, 20]; a1 = SetPrecision[-195/100, 20]; a2 = SetPrecision[2625/1000, 20]; a3 = SetPrecision[-744/100, 20]; b3 = SetPrecision[75/100, 20]; b4 = SetPrecision[75/10, 20]; T0 = SetPrecision[27/100, 20]; b2[T_] := a0 + a1*(T0/T) + a2*(T0/T)^2 + a3*(T0/T)^3; U[\[CapitalPhi]_, \[CapitalPhi]\[CapitalPhi]_, T_] := (-(1/2)* b2[T]*\[CapitalPhi]*\[CapitalPhi]\[CapitalPhi] - (b3/ 6)*(\[CapitalPhi]^3 + \[CapitalPhi]\[Capi

像如图的常微分方程，想对图中的x(t)与复合函数i(x(t))以t为变量进行绘图，需要怎么处理呢其中i(x）与D(x)的形式都比较复杂，小白刚上手mathmatica没几天实在搞不懂，提前谢谢各位佬们阿里嘎多思密达

怎么提取网页https://www.2bt0.com/movie/1.html内的所有文本信息

含有参数α和Q的二阶常微分方程：\left[1 - (1 + Q^2)(1 - x) + Q^2(1 - x)^2\right] f''(x) + \left[(1 + Q^2) - 2Q^2(1 - x)\right] f'(x) = \alpha Q^2 f(x) 有初始条件：f[0] == 1, f'[0] == α Q^2/(Q^2 - 1) 和边界条件f(1)=0 要怎么得到可以得到参数α和Q的关系啊，应该没有解析解 用gpt写一直报错

我把软件里面将生成的图像保存为png的形式放到word里面，然后我学长说不够清晰不能放进论文里面...

关键词：Export Import Read ReadList BinaryRead BinaryReadList DumpSave 导入 导出 Excel 用户友好 CSV XLS XLSX TXT 图片 PDF EPS …… 想想还是开个帖子整理下常见的导入导出问题及技巧。为方便修订照例放在云笔记上，请大家多点几下鼠标： note.youdao.com/noteshare?id=cca8f44b56fcab28a6b9bfb394397e21 有什么意见或已经议欢迎提出，不过回帖不看帖的照例直接删。

如何算一列点的下凸包！

一起交流

提问时请附上引发问题的代码。对标题论断有疑虑的同学请阅读下面这篇对 FindRoot::nlnum 警告的简要解说： https://zhidao.baidu.com/question/810397720290271692.html

其实本来想用的标题是《【交流】“这些问题你换软件也没用！”——受限于数学发展而无法或很难求解却经常有人想用Mathematica解的问题》。Mathematica作为一个软件，必然无法超越现有的数学理论，也就是说，对于尚未在理论上得到很好解决的问题，Mathematica是不太可能直接求解出来的——这话说出来大家都明白，但是，一般用户常常很难分辨哪些问题目前在理论上是很难或无法求解的，故开此交流帖，大家可以来谈谈自己所知道的相关问题。为避免

三次方怎么解出有效的解析解呀，难道是要简化参数了吗？好难呀，想做一下两种模式的利润比较，结果比不出来。 Manipulate[ Module[{result}, result = Simplify[-((\[Lambda] ((-2 + \[Theta]) \[Theta] + \[Beta] \ (-1 + \[Theta]) \[Xi])^2)/( 4 (\[Theta]^2 + (-1 + \[Theta]) (\[Alpha]s + \[Theta] (-1 + \ \[Alpha]d + \[Theta])) \[Lambda]))) + ( k ((-2 + \[Theta]) \[Theta] + \[Beta] (-1 + \[Theta]) \[Lambda] \ \[Xi])^2)/( 4 k (\[Alpha]s (-1 + \[Theta]) + \[Theta] - \[Alpha]d \[Theta] + \ \[Theta]^2 (-1 + \[Alpha]d + \[Theta])) - 2 \[Beta]^2 (-1 + \[Theta])^2

为方便修订，文章放在了云笔记上，请大家多点一下鼠标。有什么意见或建议欢迎提出： http://note.youdao.com/noteshare?id=a9c14381c11b44ad47980c44c200

请问这种可以显示子单元的符号是怎么跳出来的

也差不多该针对语义式匹配（Semantic Pattern Matching）专开个帖子了。为方便修订，照例开个云笔记，请各位多点一下鼠标： http://note.youdao.com/noteshare?id=77b86cb08aaf29ce8747990c9e71 有什么意见或建议欢迎提出。

如图 定义了一个函数H（数学上复杂，不必理会），对其关于两个变量l,p画柱状图。由于柱子高度差异比较大，想用对数刻度。但是用对数刻度后，发现柱子高度从1往下增长，感觉应该从0往上增长才对。 求助各位热心吧友 找一下问题出现在哪里，非常感激，谢谢~ 代码单独粘一楼，请吧主勿删，感谢~

{a, b, K, A, S, P, \!\(TraditionalForm\` \*SubscriptBox[\(c\), \(I\)]\), \!\(TraditionalForm\` \*SubscriptBox[\(c\), \(u\)]\), Subscript[c, S], Subscript[c, T], Subscript[c, 1], Subscript[c, E], h, \!\(TraditionalForm\` \*SuperscriptBox[\(h\), \(\[Prime]\)]\), Subscript[h, R], Subscript[h, E], X, DD, R, Subscript[t, T], m, \!\(TraditionalForm\`\[Rho]\)} = {50000, 5, 100, 200, 100, 50, 0.5, 25, 20, 2, 5, 40, 5, 4, 6, 8, 175200, 50000, 50000, 0.01, 0.2, 0.02} Subscript[Z, 1] = P - Subscript[c, u] - Subscript[c, I] - (1 - m) (Subscript[c, 1] + 2 Subscript[c, T] + Derivative[1][h] Subscript[t, T])

各位大神，想画出U0,1[Y]关于Y的函数图像，程序没报错但图出不来，请各位指教，代码及截图如下： \[CapitalGamma] = 2; \!\(\*OverscriptBox[\(\[Kappa]\), \(_\)]\) = 40; \[Beta] = 0.01; \ \[CapitalLambda] = 0.8; Pe = 0.05; \[Gamma] = 0.02; \[Chi] = 0.5; Subscript[\[Kappa], 1] = 1 - Tanh[ \!\(\*OverscriptBox[\(\[Kappa]\), \(_\)]\)]/ \!\(\*OverscriptBox[\(\[Kappa]\), \(_\)]\) + \!\(\* SubsuperscriptBox[ StyleBox[ RowBox[{"\[Epsilon]", StyleBox["De", "TI"]}]], "\[Kappa]", "2"]*\(\((Cosh[3* \*OverscriptBox[\(\[Kappa]

r = 5*10^-5(*半径：m*); \[Alpha] = 5.5*10^-7(*线膨胀系数：1/℃*); k = 1.3(*导热系数：W/(m*℃)*); \[Rho] = 2650(*石英密度：KG/m^3*); L = 20(*长度：m*); m = \[Pi]*r^2*L*\[Rho](*质量：KG*); c = 966(*石英比热容J/KG*℃*); S = Integrate[Integrate[r, {\[Theta], 0, 2 \[Pi]}], {l, 0, 20}]; sol = NDSolve[{ W'[t] == (1/2 - (W[t]/(m*c)))*k*S/r, f[t] == t/2, g[t] == Integrate[W[t]/(m*c), t]; g[0] == 0, f[0] == 0, W[0] == 0 }, {W[t], f[t], g[t]}, {t, 0, 160}]; Plot[Evaluate[{f[t], g[t], W[t]} /. sol], {t, 0, 160}, PlotLegends -> {"f[t]", "g[t]", &q

Integrate[\[Sigma]/(4 Pi \[CurlyEpsilon]) R (x - R Cos[\[Theta]])/((x - R Cos[\[Theta]])^2 + (y - R Sin[\[Theta]])^2 + z^2)^(3/2), {\[Theta], 0, 2 Pi}] 本科毕业论文里面的，本垃圾不知道为啥算出来是这个

这个也是讲烂了的内容，但吧里好像没有标题足够显眼的帖，姑且专开一帖，方便索引。 比较常用的“超级函数”基本都有与之对应的数值求解器。如标题所说，如果符号解没算出来并且符号解对你来说不是必要的，那么不妨去试试数值解。 只写这么几行的话内容好像太少，那就再稍微谈一下比较常见的成对的符号/数值求解器及其特性。 对于代数方程：Solve（虽然现在已经发展成了超级求解器，但它主要是个多项式方程求解器）/NSolve（引入时间非

通过远程桌面打开另外一个电脑的mma 页面空白显示，输入框是空白的？啥也看不到，有没有哪位前辈遇到过？

各位大神，想画出θf0，0（x）关于x的函数图像，程序没报错但图出不来，请各位指教，代码如下： \[CapitalGamma] = 1.5; \!\(\*OverscriptBox[\(\[Kappa]\), \(_\)]\) = 40; \[Beta] = 0.01; \ Subscript[\[Epsilon]De, \!\(\*OverscriptBox[\(\[Kappa]\), \(_\)]\)] = 1; \[CapitalLambda] = \ 0.8; Pe = 0.05; \[Gamma] = 0.01; Subscript[\[Kappa], 1] = 1 - Tanh[ \!\(\*OverscriptBox[\(\[Kappa]\), \(_\)]\)]/ \!\(\*OverscriptBox[\(\[Kappa]\), \(_\)]\) + Subscript[\[Epsilon]De, \!\(\*OverscriptBox[\(\[Kappa]\), \(_\)]\)]*(Cosh[3* \!\(\*OverscriptBox[\(\[Kappa]\), \(_\)]\)] - Sinh

近日在写毕业设计，有一个学术上的问题，使得StreamPlot函数显示的结果不够严谨。 参考文献上的源程序： StreamPlot[ {NIntegrate[(1 - y*Cos[\[Alpha]])/( 2 \[Pi] (1 + y^2 + x^2 - 2*y*Cos[\[Alpha]])^1.5), {\[Alpha], 0, 2 \[Pi]}], NIntegrate[(x*Cos[\[Alpha]])/( 2 \[Pi] (1 + y^2 + x^2 - 2*y*Cos[\[Alpha]])^1.5), {\[Alpha], 0, 2 \[Pi]}]}, {x, -2, 2}, {y, -2, 2}, StreamScale -> 0.05] 绘制的是两条通电直导线产生的磁感线，绘制的平面是导线的横截面。 然而，磁感线是闭合的，我们一般认为磁感线应该是一整条曲线，

Clear[lbk]"`*"[rbk]a[lbk]x_[rbk] := Cos[lbk]Pi*x[rbk]; A = E^Sqrt[lbk]1/2[rbk]; B = 0; j = 1; M = 7; (*u=\[lbk]ExponentialE[rbk]^Sqrt[lbk]1/2-x[rbk]*)uuu[lbk]x[rbk] = E^Sqrt[lbk]1/2[rbk] + b0;f[lbk]x_[rbk] := -(Pi*Sin[lbk]Pi*x[rbk] + Cos[lbk]Pi*x[rbk]*1/(2 Sqrt[lbk](1/2 - x)[rbk]))*uuu[lbk]x[rbk];step = 1;(*ai*)fs = Normal[lbk]Series[lbk](f[lbk]x[rbk]) /. x -> (1/2 - x), {x, 0, 2}[rbk][rbk];gamma = Exponent[lbk]fs, x, List[rbk];b = gamma[lbk][lbk]j[rbk][rbk];c = Coefficient[lbk]fs, x, b[rbk];imax = IntegerPart[lbk]M - 2 - b[rbk];alpj = Table[lbk] 2 + b + 2*i, {i, 0, imax}[rbk];(*\

各位好，本人在用NDSolve求解PDEs时，将结果画图得到如图结果，一半没画出来，画出来的一半极为平坦。但程序没有任何报错，想请问各位有无办法排查。源码如下： {Subscript[\[Gamma], S], Subscript[\[Gamma], L]} = {0, 0}; {Subscript[v, S], Subscript[v, L]} = {0.8, 0.1}; {Subscript[\[Omega], S], Subscript[\[Omega], L]} = {0.7, 0.3}; {Subscript[k, i], Subscript[k, L]} = {1.9, 1.68}; {Subscript[\[Beta], 1], Subscript[\[Beta], 2]} = {0.2, -0.04}; \[Alpha] = 1; eqns = {\!\( \*SubscriptBox[\(\[PartialD]\), \(t\)]\(a[x, t]\)\) + Subscript[\[Gamm

现在我对这个微分方程的解猜出了一个形式解，如图所示，现在我想要级数展开然后求解前面的常系数，但是使用Series函数只能得到幂级数的展开解，带log z的这些项体现不出来，我要怎么做才能得到我想要的形式展开呢？ 以下是代码： In[4]:= SetAttributes[{u1, u0, zh, v4}, Constant] In[2]:= eq1 = y''[z] + (-3/z + 2*u1^2*z + 4*z^3/(z^4 - zh^4)) y'[z] + 1/z^2*zh^4/(zh^4 - z^4) (3*y[z] - 4*v4*y[z]^3) == 0 Out[2]= (zh^4 (3 y[z] - 4 v4 y[z]^3))/( z^2 (-z^4 + zh^4)) + (-(3/z) + 2 u1^2 z + (4 z^3)/( z^4 - zh^4)

大佬们，这个我算不出来，能帮我看看吗？eqn = D[lbk]u[lbk]t, x, y[rbk], t[rbk] == D[lbk]u[lbk]t, x, y[rbk], x, x[rbk] + D[lbk]u[lbk]t, x, y[rbk], y, y[rbk] - u[lbk]t, x, y[rbk] + u[lbk]t, x, y[rbk]^3;ic = {u[lbk]0, x, y[rbk] == Exp[lbk]-(x^2 + y^2)[rbk], Derivative[lbk]1, 0, 0[rbk][lbk]u[rbk][lbk]0, x, y[rbk] == 0};sol = NDSolve[lbk]{eqn, ic}, u, {t, 0, 1}, {x, -5, 5}, {y, -5, 5}[rbk];Plot3D[lbk]Evaluate[lbk]u[lbk]t, x, y[rbk] /. sol[rbk], {t, 0, 1}, {x, -5, 5}, {y, -5, 5}, PlotRange -> All[rbk]

um=4; list=Table[ParametricPlot3D[{Cos[u]+n/2,Sin[u]+n/2,v}, {v,0,1},{u,0,um},Axes->False,Boxed->False, PlotStyle->Directive[FaceForm[Green,Texture["ABCDEFG"]]], TextureCoordinateFunction->({2#5,#4}&),Lighting->"Neutral", Mesh->None,PlotRange->All,ViewPoint->{-30,10,20}, ViewAngle->.025],{n,0,2}] Show[list[[1 ;; 3]]]

看了帮助文档和网上的经验贴都没解决问题 无奈之下向贴吧各位大佬们求助 谢谢各位 报错截图如下 具体代码文本如下 v0 = 100; Rp = 0.1; R = 0.025; \[Nu] = 0.3; \[Mu] = 8.5; EE = 50*10^9; mp = 0.03; h = 0.005; G = 2*10^10; \[Rho]s = 1200; k1 = 1; k2 = 1; \[Sigma]ys = 50*10^6; n1 = 1; n2 = 1; \[Omega][r_, t_] := \[Omega]0[ t]*((2*r^3)/\[Eta][t]^3 - (3*r^2)/\[Eta][t]^2 + 1); \[Gamma]rz[r_, t_] := \[Gamma]rz0[t]*Sin[(Pi*r)/\[Eta][t]]^2; \[Rho][Z] = 0.8*Z + 0.1; \[Rho]Z0 = Integrate[\[Rho][Z], {Z, 0, 1}]; \[Sigma]0[\[Rho]] = k1*\[Sigma]ys*\[Rho]Z0^n1; CC[\[

请问有大佬知道从哪里可以找到现成的高阶隐式龙格库塔方法的butcher tableau吗? 之前学计算物理时候好像RK方法的系数是用Taylor展开推出来的，但IRK的表感觉有点难推

Clear["Global`*"]; delta[a_] := Z^2*(1 - 2 a)/(a*(Z + 1)); a11[a_] := (3*Pi/32)*(288 + 604*Sqrt[2]*delta[a] + 217*delta[a]^2)/ 4/(72 + 61*Sqrt[2]*delta[a] + 16*delta[a]^2); dc[a_] := Piecewise[{{1, a <= 0}, {a11[a]/(1 - a*(2 - mi/mI*(Z + 1)))/(a*(Z - 1) + 1), a > 0 && a <= 1/2}, {dc[1/2], a > 1/2}}]; dp[a_] := dc[a]*a*(1 - 2 a)*(Z - 1)^2/(Z - a*(Z - 1)); D0 = 278*4/100*Sqrt[2]*5^(5/2)/(2*1.81)*(Z + 1)/Z^2; (*cm^2/s*) Z = 6; mi = 2; mI = 36; a0[z_] := Piecewise[{{1/2, z < 0}, {0, z >= 0}}]; sol = NDSolve[{D[a[z, t], t] == D[D0*(dc[a[z, t]] + dp[a[z, t]])*D[a[