Rascal

Preview onto version 0.3.4

SourceForge Logo Here you can do some experiments with Rascal. The computation time is restricted to about 10 secs to prevent high load on server. Have a look at the documentation or at the examples. If you see wrong plots reload may help. You job is cancelled as soon someone else accesses this page, so consider a local installation.
Rascal-toolbox-0.3.4 (C) 2001,2002 Sebastian Ritterbusch

Enter 'help' for help or 'licence' for information about licence and warranty.

 "rascal.rc loaded, predefined e, i, pi"
>"Middle value and higher order interval evaluation";
>
>M(f,Y,z)=f([z,z])+f'(Y)*(Y-z);
>H(f,Y,z)=f([z,z])+f'([z,z])*(Y-z)+f''(Y)*sqr(Y-z)/2;
>
>p=(X-17)^9
 X^9-153*X^8+10404*X^7-412692*X^6+10523646*X^5-178901982*X^4+2027555796*X^3-14772192228*X^2+62781816969*X-118587876497
>
>p(17+1/10)*1.
 1.0000000000E-009
>p(17.1)
   -0.0005035400
>p([17.1,17.1])
 [  -0.0078430176,   0.0070800782]
>
>w=1e-4;Y=[17.1,17.1]+w*[-1,1]/2
 [  17.0999499999,  17.1000500001]
>
>p(Y)
 [-9.1144963283E+007,9.1145229790E+007]
>M(p,Y,mid(Y))
 [-1.2027985001E+003,1.2028008806E+003]
>H(p,Y,mid(Y))
 [  -0.0092670592,   0.0116474509]
>
>y(w)=[17.1,17.1]+w*[-1,1]/2;
>order(t)=[log(diam(p(y(10^t)))) log(diam(M(p,y(10^t),17.1))) log(diam(H(p,y(10^t),17.1)))]/log(10);
>
>order(-4)
 [   8.2607633049    3.3812247688   -1.5371185426]
>
>plotGrid(1);plot(order,-12,0);


>"log-log diagram in width of incoming interval to resulting interval of the three methods, ";
>"whose order can be read from the diagram.";
>


Questions, Problems, Bugs? Please use this form.