Mathematica 9 is now available
Student Support Forum
-----
Student Support Forum: 'Solve Error..?' topicStudent Support Forum > General > "Solve Error..?"

Next Comment >Help | Reply To Topic
Author Comment/Response
June E.S
10/31/07 6:15pm

I can't seem to find why I keep getting these errors when I try to do the following..

ClearAll;
<< Statistics`ContinuousDistributions`
<< Graphics`Graphics`
<< Statistics`DescriptiveStatistics`
Clear[a, b, s, t, v, alpha, mdist]
mdist = LogNormalDistribution[a, b];
lambdaNGL = 1/200974.
lambdaNGU = 2*lambdaNGL
s = Solve[{Quantile[mdist, 0.1] == lambdaNGL, Quantile[mdist, .50] == \
lambdaNGU}, {a, b}, Reals]
a = s[[1, 2, 2]]
b = s[[1, 1, 2]]
meanX = Exp[a + 0.5b^2]
varX = meanX^2(Exp[b^2] - 1)
tau = (1/b)^2
Plot[PDF[mdist, z], {z, 0, .0001}]
Export["LogNormPDFplot-LambdaKaw.eps", %]
Plot[CDF[mdist, z], {z, 0, 0.0001}, PlotRange -> {{0, .0001}, {0, 1}}]

The outputs are..

Out[49]=
\!\(4.975768009792311`*^-6\)

Out[50]=
\!\(9.951536019584623`*^-6\)

From In[43]:=
\!\(\*
RowBox[{\(Solve::"dinv"\), \(\(:\)\(\ \)\), "\<\"The expression \\!\\(
LogNormalDistribution[\\(\\(a, b\\)\\)]\\) involves unknowns in more \
than one argument, so inverse functions cannot be used. \\!\\(\\*ButtonBox[\\\
\"More…\\\", ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \
ButtonData:>\\\"Solve::dinv\\\"]\\)\"\>"}]\)

From In[43]:=
\!\(\*
RowBox[{\(Solve::"dinv"\), \(\(:\)\(\ \)\), "\<\"The expression \\!\\(
LogNormalDistribution[\\(\\(a, b\\)\\)]\\) involves unknowns in more \
than one argument, so inverse functions cannot be used. \\!\\(\\*ButtonBox[\\\
\"More…\\\", ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \
ButtonData:>\\\"Solve::dinv\\\"]\\)\"\>"}]\)

Out[51]=
\!\(Solve[{Quantile[
LogNormalDistribution[a, b], 0.1`] \[Equal] 4.975768009792311`*^-6, \
Quantile[LogNormalDistribution[a, b], 0.5`] \[Equal] 9.951536019584623`*^-6}, \
{a, b}, Reals]\)

Out[52]=
\!\(9.951536019584623`*^-6\)

From In[43]:=
\!\(\*
RowBox[{\(General::"ivar"\), \(\(:\)\(\ \
\)\), "\<\"\\!\\(9.951536019584623`*^-6\\) is not a valid
variable. \\!\\(\\*ButtonBox[\\\"More…\\\", \
ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \
ButtonData:>\\\"General::ivar\\\"]\\)\"\>"}]\)

From In[43]:=
\!\(\*
RowBox[{\(General::"ivar"\), \(\(:\)\(\ \
\)\), "\<\"\\!\\(9.951536019584623`*^-6\\) is not a valid
variable. \\!\\(\\*ButtonBox[\\\"More…\\\", \
ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \
ButtonData:>\\\"General::ivar\\\"]\\)\"\>"}]\)

Out[53]=
\!\(4.975768009792311`*^-6\)

Out[54]=
1.00001

Out[55]=
\!\(2.475868826341116`*^-11\)

Out[56]=
\!\(4.0390548676`*^10\)

From In[43]:=
\!\(\*
RowBox[{\(Plot::"plnr"\), \(\(:\)\(\ \)\), "\<\"\\!\\(PDF[\\(\\(mdist,
z\\)\\)]\\) is not a
machine-size real number at \\!\\(z\\) = \
\\!\\(4.166666666666667`*^-12\\). \\!\\(\\*ButtonBox[\\\"More…\\\", \
ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \
ButtonData:>\\\"Plot::plnr\\\"]\\)\"\>"}]\)

From In[43]:=
\!\(\*
RowBox[{\(Plot::"plnr"\), \(\(:\)\(\ \)\), "\<\"\\!\\(PDF[\\(\\(mdist,
z\\)\\)]\\) is not a
machine-size real number at \\!\\(z\\) = \
\\!\\(4.056699157291579`*^-6\\). \\!\\(\\*ButtonBox[\\\"More…\\\", \
ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \
ButtonData:>\\\"Plot::plnr\\\"]\\)\"\>"}]\)

From In[43]:=
\!\(\*
RowBox[{\(Plot::"plnr"\), \(\(:\)\(\ \)\), "\<\"\\!\\(PDF[\\(\\(mdist,
z\\)\\)]\\) is not a
machine-size real number at \\!\\(z\\) = \
\\!\\(8.480879985937368`*^-6\\). \\!\\(\\*ButtonBox[\\\"More…\\\", \
ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \
ButtonData:>\\\"Plot::plnr\\\"]\\)\"\>"}]\)

From In[43]:=
\!\(\*
RowBox[{\(General::"stop"\), \(\(:\)\(\ \)\), "\<\"Further output
of \\!\\(Plot :: \\\"plnr\\\"\\)
will be suppressed during this calculation. \
\\!\\(\\*ButtonBox[\\\"More…\\\", ButtonStyle->\\\"RefGuideLinkText\\\", \
ButtonFrame->None, ButtonData:>\\\"General::stop\\\"]\\)\"\>"}]\)

From In[43]:=
GraphicsData["PostScript", "\<\
%!
%%Creator: Mathematica
%%AspectRatio: .61803
MathPictureStart
/Mabs {
Mgmatrix idtransform
Mtmatrix dtransform
} bind def
/Mabsadd { Mabs
3 -1 roll add
3 1 roll add
exch } bind def
%% Graphics
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10 scalefont setfont
% Scaling calculations
0.0238095 0.952381 0.0147151 0.588604 [
[.21429 .00222 -9 -9 ]
[.21429 .00222 9 0 ]
[.40476 .00222 -9 -9 ]
[.40476 .00222 9 0 ]
[.59524 .00222 -9 -9 ]
[.59524 .00222 9 0 ]
[.78571 .00222 -9 -9 ]
[.78571 .00222 9 0 ]
[.97619 .00222 -3 -9 ]
[.97619 .00222 3 0 ]
[.01131 .13244 -18 -4.5 ]
[.01131 .13244 0 4.5 ]
[.01131 .25016 -18 -4.5 ]
[.01131 .25016 0 4.5 ]
[.01131 .36788 -18 -4.5 ]
[.01131 .36788 0 4.5 ]
[.01131 .4856 -18 -4.5 ]
[.01131 .4856 0 4.5 ]
[.01131 .60332 -6 -4.5 ]
[.01131 .60332 0 4.5 ]
[ 0 0 0 0 ]
[ 1 .61803 0 0 ]
] MathScale
% Start of Graphics
1 setlinecap
1 setlinejoin
newpath
0 g
.25 Mabswid
[ ] 0 setdash
.21429 .01472 m
.21429 .02097 L
s
[(0.2)] .21429 .00222 0 1 Mshowa
.40476 .01472 m
.40476 .02097 L
s
[(0.4)] .40476 .00222 0 1 Mshowa
.59524 .01472 m
.59524 .02097 L
s
[(0.6)] .59524 .00222 0 1 Mshowa
.78571 .01472 m
.78571 .02097 L
s
[(0.8)] .78571 .00222 0 1 Mshowa
.97619 .01472 m
.97619 .02097 L
s
[(1)] .97619 .00222 0 1 Mshowa
.125 Mabswid
.07143 .01472 m
.07143 .01847 L
s
.11905 .01472 m
.11905 .01847 L
s
.16667 .01472 m
.16667 .01847 L
s
.2619 .01472 m
.2619 .01847 L
s
.30952 .01472 m
.30952 .01847 L
s
.35714 .01472 m
.35714 .01847 L
s
.45238 .01472 m
.45238 .01847 L
s
.5 .01472 m
.5 .01847 L
s
.54762 .01472 m
.54762 .01847 L
s
.64286 .01472 m
.64286 .01847 L
s
.69048 .01472 m
.69048 .01847 L
s
.7381 .01472 m
.7381 .01847 L
s
.83333 .01472 m
.83333 .01847 L
s
.88095 .01472 m
.88095 .01847 L
s
.92857 .01472 m
.92857 .01847 L
s
.25 Mabswid
0 .01472 m
1 .01472 L
s
.02381 .13244 m
.03006 .13244 L
s
[(0.2)] .01131 .13244 1 0 Mshowa
.02381 .25016 m
.03006 .25016 L
s
[(0.4)] .01131 .25016 1 0 Mshowa
.02381 .36788 m
.03006 .36788 L
s
[(0.6)] .01131 .36788 1 0 Mshowa
.02381 .4856 m
.03006 .4856 L
s
[(0.8)] .01131 .4856 1 0 Mshowa
.02381 .60332 m
.03006 .60332 L
s
[(1)] .01131 .60332 1 0 Mshowa
.125 Mabswid
.02381 .04415 m
.02756 .04415 L
s
.02381 .07358 m
.02756 .07358 L
s
.02381 .10301 m
.02756 .10301 L
s
.02381 .16187 m
.02756 .16187 L
s
.02381 .1913 m
.02756 .1913 L
s
.02381 .22073 m
.02756 .22073 L
s
.02381 .27959 m
.02756 .27959 L
s
.02381 .30902 m
.02756 .30902 L
s
.02381 .33845 m
.02756 .33845 L
s
.02381 .39731 m
.02756 .39731 L
s
.02381 .42674 m
.02756 .42674 L
s
.02381 .45617 m
.02756 .45617 L
s
.02381 .51503 m
.02756 .51503 L
s
.02381 .54446 m
.02756 .54446 L
s
.02381 .57389 m
.02756 .57389 L
s
.25 Mabswid
.02381 0 m
.02381 .61803 L
s
0 0 m
1 0 L
1 .61803 L
0 .61803 L
closepath
clip
newpath
% End of Graphics
MathPictureEnd
\
\>"]

Out[57]=
⁃Graphics⁃

Out[58]=
LogNormPDFplot-LambdaKaw.eps

From In[43]:=
\!\(\*
RowBox[{\(Plot::"plnr"\), \(\(:\)\(\ \)\), "\<\"\\!\\(CDF[\\(\\(mdist,
z\\)\\)]\\) is not a
machine-size real number at \\!\\(z\\) = \
\\!\\(4.166666666666667`*^-12\\). \\!\\(\\*ButtonBox[\\\"More…\\\", \
ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \
ButtonData:>\\\"Plot::plnr\\\"]\\)\"\>"}]\)

From In[43]:=
\!\(\*
RowBox[{\(Plot::"plnr"\), \(\(:\)\(\ \)\), "\<\"\\!\\(CDF[\\(\\(mdist,
z\\)\\)]\\) is not a
machine-size real number at \\!\\(z\\) = \
\\!\\(4.056699157291579`*^-6\\). \\!\\(\\*ButtonBox[\\\"More…\\\", \
ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \
ButtonData:>\\\"Plot::plnr\\\"]\\)\"\>"}]\)

From In[43]:=
\!\(\*
RowBox[{\(Plot::"plnr"\), \(\(:\)\(\ \)\), "\<\"\\!\\(CDF[\\(\\(mdist,
z\\)\\)]\\) is not a
machine-size real number at \\!\\(z\\) = \
\\!\\(8.480879985937368`*^-6\\). \\!\\(\\*ButtonBox[\\\"More…\\\", \
ButtonStyle->\\\"RefGuideLinkText\\\", ButtonFrame->None, \
ButtonData:>\\\"Plot::plnr\\\"]\\)\"\>"}]\)

From In[43]:=
\!\(\*
RowBox[{\(General::"stop"\), \(\(:\)\(\ \)\), "\<\"Further output
of \\!\\(Plot :: \\\"plnr\\\"\\)
will be suppressed during this calculation. \
\\!\\(\\*ButtonBox[\\\"More…\\\", ButtonStyle->\\\"RefGuideLinkText\\\", \
ButtonFrame->None, ButtonData:>\\\"General::stop\\\"]\\)\"\>"}]\)

From In[43]:=
GraphicsData["PostScript", "\<\
%!
%%Creator: Mathematica
%%AspectRatio: .61803
MathPictureStart
/Mabs {
Mgmatrix idtransform
Mtmatrix dtransform
} bind def
/Mabsadd { Mabs
3 -1 roll add
3 1 roll add
exch } bind def
%% Graphics
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10 scalefont setfont
% Scaling calculations
-2.39609e-017 10000 0 0.618034 [
[.2 -0.0125 -21 -9 ]
[.2 -0.0125 21 0 ]
[.4 -0.0125 -21 -9 ]
[.4 -0.0125 21 0 ]
[.6 -0.0125 -21 -9 ]
[.6 -0.0125 21 0 ]
[.8 -0.0125 -21 -9 ]
[.8 -0.0125 21 0 ]
[1 -0.0125 -18 -9 ]
[1 -0.0125 18 0 ]
[-0.0125 .12361 -18 -4.5 ]
[-0.0125 .12361 0 4.5 ]
[-0.0125 .24721 -18 -4.5 ]
[-0.0125 .24721 0 4.5 ]
[-0.0125 .37082 -18 -4.5 ]
[-0.0125 .37082 0 4.5 ]
[-0.0125 .49443 -18 -4.5 ]
[-0.0125 .49443 0 4.5 ]
[-0.0125 .61803 -6 -4.5 ]
[-0.0125 .61803 0 4.5 ]
[ 0 0 0 0 ]
[ 1 .61803 0 0 ]
] MathScale
% Start of Graphics
1 setlinecap
1 setlinejoin
newpath
0 g
.25 Mabswid
[ ] 0 setdash
.2 0 m
.2 .00625 L
s
[(0.00002)] .2 -0.0125 0 1 Mshowa
.4 0 m
.4 .00625 L
s
[(0.00004)] .4 -0.0125 0 1 Mshowa
.6 0 m
.6 .00625 L
s
[(0.00006)] .6 -0.0125 0 1 Mshowa
.8 0 m
.8 .00625 L
s
[(0.00008)] .8 -0.0125 0 1 Mshowa
1 0 m
1 .00625 L
s
[(0.0001)] 1 -0.0125 0 1 Mshowa
.125 Mabswid
.05 0 m
.05 .00375 L
s
.1 0 m
.1 .00375 L
s
.15 0 m
.15 .00375 L
s
.25 0 m
.25 .00375 L
s
.3 0 m
.3 .00375 L
s
.35 0 m
.35 .00375 L
s
.45 0 m
.45 .00375 L
s
.5 0 m
.5 .00375 L
s
.55 0 m
.55 .00375 L
s
.65 0 m
.65 .00375 L
s
.7 0 m
.7 .00375 L
s
.75 0 m
.75 .00375 L
s
.85 0 m
.85 .00375 L
s
.9 0 m
.9 .00375 L
s
.95 0 m
.95 .00375 L
s
.25 Mabswid
0 0 m
1 0 L
s
0 .12361 m
.00625 .12361 L
s
[(0.2)] -0.0125 .12361 1 0 Mshowa
0 .24721 m
.00625 .24721 L
s
[(0.4)] -0.0125 .24721 1 0 Mshowa
0 .37082 m
.00625 .37082 L
s
[(0.6)] -0.0125 .37082 1 0 Mshowa
0 .49443 m
.00625 .49443 L
s
[(0.8)] -0.0125 .49443 1 0 Mshowa
0 .61803 m
.00625 .61803 L
s
[(1)] -0.0125 .61803 1 0 Mshowa
.125 Mabswid
0 .0309 m
.00375 .0309 L
s
0 .0618 m
.00375 .0618 L
s
0 .09271 m
.00375 .09271 L
s
0 .15451 m
.00375 .15451 L
s
0 .18541 m
.00375 .18541 L
s
0 .21631 m
.00375 .21631 L
s
0 .27812 m
.00375 .27812 L
s
0 .30902 m
.00375 .30902 L
s
0 .33992 m
.00375 .33992 L
s
0 .40172 m
.00375 .40172 L
s
0 .43262 m
.00375 .43262 L
s
0 .46353 m
.00375 .46353 L
s
0 .52533 m
.00375 .52533 L
s
0 .55623 m
.00375 .55623 L
s
0 .58713 m
.00375 .58713 L
s
.25 Mabswid
0 0 m
0 .61803 L
s
0 0 m
1 0 L
1 .61803 L
0 .61803 L
closepath
clip
newpath
% End of Graphics
MathPictureEnd
\
\>"]

Out[59]=
⁃Graphics⁃



URL: ,

Subject (listing for 'Solve Error..?')
Author Date Posted
Solve Error..? June E.S 10/31/07 6:15pm
Re: Solve Error..? yehuda ben-s... 11/04/07 01:09am
Next Comment >Help | Reply To Topic