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"
>" Verification of the solution of linear (point) systems ";
>
>A=[1e+7 1;1e+7 1.0001];
>b=[2;7];
>R=1/A
 [   0.0010001000   -0.0010000000;-1.0000000000E+004 1.0000000000E+004]
>
>x=R*b
 [  -0.0049998000;5.0000000000E+004]
>A*x-b
 [7.2759576142E-012;7.2759576142E-012]
>
>" Iteration xk=(1-R*A)*xk+R*b starting with inflated x0:";
>xk=x*[0.8,1.2];
>i=0;
>C=(1-R*[1,1]*A);
>z=R*[1,1]*b;
>"[# Enclosure xk]";
>xn=C*xk+z;[i++ less(xn(1),xk(1))&less(xn(2),xk(2)) xk=xn]
 [1 1 [[  -0.0049998001,  -0.0049997999];[4.9999999999E+004,5.0000000001E+004]]]
>xn=C*xk+z;[i++ less(xn(1),xk(1))&less(xn(2),xk(2)) xk=xn]
 [2 1 [[  -0.0049998001,  -0.0049997999];[4.9999999999E+004,5.0000000001E+004]]]
>xn=C*xk+z;[i++ less(xn(1),xk(1))&less(xn(2),xk(2)) xk=xn]
 [3 0 [[  -0.0049998001,  -0.0049997999];[4.9999999999E+004,5.0000000001E+004]]]
>xn=C*xk+z;[i++ less(xn(1),xk(1))&less(xn(2),xk(2)) xk=xn]
 [4 0 [[  -0.0049998001,  -0.0049997999];[4.9999999999E+004,5.0000000001E+004]]]
>
>yk=(z=R*(b-A*[1,1]*x))*[0.8,1.2];
>C=(1-R*[1,1]*A);
>
>" Krawczyk Iteration yk=(1-R*A)*yk+R*(b-A*x) starting with inflated y0:";
>yk
 [[-1.7463171389E-014,2.6193447412E-014];[-2.6193447412E-007,1.7462298275E-007]]
>i=0;
>"[# Enclosure yk+x]";
>yn=C*yk+z;[i++ less(yn(1),yk(1))&less(yn(2),yk(2)) (yk=yn)+x]
 [1 1 [[  -0.0049998001,  -0.0049997999];[4.9999999999E+004,5.0000000001E+004]]]
>yn=C*yk+z;[i++ less(yn(1),yk(1))&less(yn(2),yk(2)) (yk=yn)+x]
 [2 1 [[  -0.0049998001,  -0.0049997999];[4.9999999999E+004,5.0000000001E+004]]]
>yn=C*yk+z;[i++ less(yn(1),yk(1))&less(yn(2),yk(2)) (yk=yn)+x]
 [3 0 [[  -0.0049998001,  -0.0049997999];[4.9999999999E+004,5.0000000001E+004]]]
>yn=C*yk+z;[i++ less(yn(1),yk(1))&less(yn(2),yk(2)) (yk=yn)+x]
 [4 0 [[  -0.0049998001,  -0.0049997999];[4.9999999999E+004,5.0000000001E+004]]]
>


Questions, Problems, Bugs? Please use this form.