How to solve for gamma
5 次查看(过去 30 天)
显示 更早的评论
Hello, I'm trying to solve for gamma in this equation, I was hoping if someone can help out.
considering everything else is a constant that can be put as an input.
2 个评论
Torsten
2022-12-2
An analytical expression for gamma does not seem possible.
Use "fzero" or "fsolve" as numerical root finders for
left-hand side(gamma) - right-hand side(gamma) = 0
采纳的回答
Torsten
2022-12-2
编辑:Torsten
2022-12-2
format long
dbs=5.3e-10;
das=1.01e-9;
cb=4e3;
ca=1e3;
x0 =0;
options = optimset('TolX',1e-14,'TolFun',1e-14);
x = fsolve(@(x)myfun(x,dbs,das,cb,ca),x0,options)
myfun(x,dbs,das,cb,ca)
function F = myfun(x,dbs,das,cb,ca)
F = (1-erf(sqrt(x/dbs)))-(((cb/ca)*sqrt(dbs/das))*(erf(sqrt(x/das)))*(exp((x/das)-(x/dbs))));
end
5 个评论
Torsten
2022-12-5
编辑:Torsten
2022-12-5
Seems your equation has multiple solutions.
A good idea is always to plot the function to have a good guess for the solution.
dbs=5.3e-10;
das=1.01e-9;
cb=4e3;
ca=1e3;
F = @(x)(1-erf(sqrt(x/dbs)))-(((cb/ca)*sqrt(dbs/das))*(erf(sqrt(x/das))).*(exp((x/das)-(x/dbs))));
x = 0:1e-12:2e-10;
plot(x,F(x))
sol = fzero(F,[0 1e-9])
F(sol)
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Gamma Functions 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!