Помогите найти критическую толщину

[drWho]

Всем привет!
Помогите, пожалуйста, получить график зависимости критической толщины пленки от состава, если
$a_0=(6.05-0.4x)$ - параметр решетки тройного раствора,
$a_s=5.65$ - параметр решетки подложки;
$v=1/3$ коэф.Пуассона;
$\varphi =60^\circ$;
$\Theta=54.74^\circ$;
$b={a_0/(\sqrt{2})}$;
$f=(a_0-a_s)/a_0$;

$h_c=b(1-v\cos^2(\varphi))\frac {1}{ 8\pi f(1+v)\sin(\varphi)\cos(\Theta)} $$(\ln(\frac {h_c} {b\sin(\Theta)})+1)$ - критическая толщина

-- 22.02.2020, 12:44 --

clc;
clear;
x=linspace(0, 0.5, 100);
a0=(6.0583-0.405.*x);
as=5.65325;
b=a0./2;
v=1/3;
fi=60*pi/180;
theta=54.74*pi/180;
f=(a0-as)./a0;
m=b.*(1-v*cos(fi)^2)./(8*pi.*f*(1+v)*cos(theta)*sin(fi));
W=0.971846*b-(0.11.*b.*lambertw(8.835*exp(8.835.*f)).*f)./f;
hc=m+W;
plot(x,hc);

Скорее всего я неправильно считаю функцию Ламберта.

[Pphantom]

Тема перемещена из форума «Помогите решить / разобраться (Ф)» в форум «Карантин»
по следующим причинам:

- неправильно набраны формулы (краткие инструкции: «Краткий FAQ по тегу [math]» и видеоролик Как записывать формулы).

Исправьте все Ваши ошибки и сообщите об этом в теме Сообщение в карантине исправлено.
Настоятельно рекомендуется ознакомиться с темами Что такое карантин и что нужно делать, чтобы там оказаться и Правила научного форума.

[Pphantom]

Тема перемещена из форума «Карантин» в форум «Помогите решить / разобраться (Ф)»