On Mon, 14 Jun 2010, Leo Butler wrote:
>
>
> On Mon, 14 Jun 2010, Hugo Coolens wrote:
>
> < Here is my input-file:
> < m:matrix([100,16000],[200,8000],[500,3636],[1000,1818],[2000,870],[5000,323],[10000,160],[20000,83],[50000,33],[100000,20]);
> < lsquares_estimates (m, [x, y], y = sqrt(a^2 + (10e7/(2*%pi*x*b))^2),[a, b],
> < initial=[10, 10]);
> < /* => [[a = 19.2898942613692, b = 9.89381119455951]]*/
> < /* b = 10 is the exact value for b.*/
> < /* I'll construct the mean square error as a function of a alone*/
> < /* and then plot it. */
> < mse : lsquares_mse (m, [x, y], y = sqrt(a^2 + (10e7/(2*%pi*x*b))^2));
> < mse1 (a) := ''(subst (b=10, mse));
> < plot2d(mse1(a), [a,0,100]);
>
> Hugo, what if you try
>
> define(mse1(a), subst(b=10,mse));
>
> instead of the line above that uses ''.
I tried this suggestion but still I get:
plot [0.:100.] '-' notitle with lines 3
^
"/home/coolens/maxout.gnuplot", line 5: warning: Skipping data file with
no valid points
plot [0.:100.] '-' notitle with lines 3
^
"/home/coolens/maxout.gnuplot", line 5: all points y value undefined!
>
> Are you 'load'ing this input file, by any chance?
I called the input-file moremax.mac
I start Maxima as follows:
maxima <moremax.mac
hugo