解不了图片中要求得值，请教大神，非常感谢！

busahuang

function [TD] = menxian2(L)
5 c. {' t- v$ w6 S" ~9 k%UNTITLED4 此处显示有关此函数的摘要
%   此处显示详细说明
& E! e3 m. R' e9 N3 v$ d1 O9 }  p! Jsyms s y x t
  O" `' |' _. z0 fs=int(exp(-(t^2/2)),t,-inf,x)*1/sqrt(2*pi);
9 x3 j* K" X; B) U) @/ ry=matlabFunction(s);
f1=finverse(s,t);
f2=matlabFunction(f1);
TD=f1(1-L);
/ i4 }( D! H3 I+ V1 C2 f3 I" q# `end
. q) y& V5 a! E# [4 i9 x6 v; O
代码如上，解不了图片中要求得值，请教大神，非常感谢！
5 a. U* H- H) G! d
关于这个函数得解，Pfa已知，求Td
5 @0 m1 ^1 t$ i9 N( q8 X

: e3 m; T3 K. {9 l5 s4 V, Y
: S7 r* c1 o* {( W4 q1 ?1 j
3 {( }& _6 u, @" b

Ferrya

5 w% v& n) R( ?2 z5 T直接用TD = norminv( 1- L );就行了
非要解方程的话用vpasolve求数值解
; V& O! x( S/ F& Z: Ksyms x t real
; y) E9 i- G: q0 a! Z0 I7 ^  Bs = int(exp(-(t^2/2)),t,-inf,x)*1/sqrt(2*pi);
+ c% F* X$ p' O7 U$ N3 b# X% d7 ^y = matlabFunction(s);
TD = vpasolve( y(x)==1-L );

busahuang

Ferrya 发表于 2020-5-28 15:03
1 n  s7 [5 J* S5 G! d# x/ K直接用TD = norminv( 1- L );就行了
6 W" v. O  Z: |" N8 M2 q9 y+ A. d  z非要解方程的话用vpasolve求数值解
syms x t real

大佬看看求这个Td用的代码跟前面这个一样吗，n已知。