別の変数の関数を変数として取る関数の 3D プロットを作成したいと思います。
私の関数は部分関数でネストされているため、全体が複雑になっています
ここに私のコードがあります:
phi0=Function[u,1.21*10^-6/((u/10^44.25)^1.01 + (u/10^44.25)^2.38)][Lx]
zc=Function[v,Piecewise[{{2.49,v>=10^45.74},{2.49*(v/10^45.74)^0.2,v<10^45.74}}]][Lx]
e=Function[uu,Piecewise[{{(1+uu)^4.62,uu<=zc},{(1+zc)^4.62*((1+uu)/(1+zc))^-1.15,uu>zc}}]][z]
Plot3D[e[z,Lx],{z,0,7},{Lx,10^42,10^47}, PlotRange->Full]
しかし、それは何もプロットしておらず、何をすべきかわかりません
編集:
ヒントをありがとう、私はこの方法で解決したと思います。エラーは出ていませんが、1点でも結果を評価するのに時間がかかります...それは正常だと思いますか?
phizero[Lx_] := 1.21/10^6/((Lx/10^44.25)^1.01 + (Lx/10^44.25)^2.38)
zc[Lx_] :=
Piecewise[{{2.49, Lx >= 10^45.74}, {2.49*(Lx/10^45.74)^0.2,
Lx < 10^45.74}}]
e[z_, Lx_] :=
Piecewise[{{(1 + z)^4.62,
z <= zc[Lx]}, {(1 + zc[Lx])^4.62*((1 + z)/(1 + zc[Lx]))^-1.15,
z > zc[Lx]}}]
phi[z_, Lx_] := phizero[Lx]*e[z, Lx]
(*D[phi[z,Lx],Lx]:=Lx*phi[z,Lx]*)
p[z_, Lx_] = Integrate[Lx*phi[z, Lx], Lx]
p[4, 10^44]