数量を含む補間データを微分，積分する 

数量データのリストに対して三次元補間を実行する．

In[1]:=

X data = N@Flatten[ Outer[Quantity[Sin[#1 + #2^2], "Pounds"] &, Range[10], Range[10]], 1]; grid = Flatten[ Outer[List, Quantity[Range[10], "Seconds"], Quantity[Range[10], "Feet"]], 1]; data = Join[grid, Transpose[{data}], 2] /. {x_, y_, z_} :> {{x, y}, z};

In[2]:=

X if = Interpolation[data]

Out[2]=

解をプロットする．

完全なWolfram言語入力を表示する入力を隠す

In[3]:=

X fun[t_?NumericQ, p_?NumericQ] := QuantityMagnitude[ Quiet[if[Quantity[t, "Seconds"], Quantity[p, "Feet"]] ] ]; Plot3D[fun[tt, p], {tt, 1, 10}, {p, 1, 10}, AxesLabel -> {"Seconds", "Feet", "Pounds"}, PlotRange -> All, ImageSize -> Medium]

Out[3]=

補間された関数に関連付けられた入力および出力の数量は，特性"QuantityUnits"を使って抽出することができる．

In[4]:=

X if["QuantityUnits"]

Out[4]=

補間された関数を微分する．

In[5]:=

X Derivative[#, 0][if][Quantity[3, "Seconds"], Quantity[1, "Meters"]] & /@ Range[3]

Out[5]=

In[6]:=

X Derivative[#, #][if][Quantity[3, "Seconds"], Quantity[1, "Meters"]] & /@ Range[3]

Out[6]=

補間された関数を積分する．

In[7]:=

X Integrate[ if[x, y], {x, Quantity[1, "Seconds"], Quantity[10, "Seconds"]}, {y, Quantity[1, "Feet"], Quantity[10, "Feet"]}]

Out[7]=