1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
|
.. Test the spline function.
function test
## Demonstrates the spline functions in RMAT.
"Click any number of points (x,y) with different x"
"The program will compute the interpolating spline"
"and the interpolating polynomial."
"Press any key, please!"
wait(180);
""
setplot(-1.1,1.1,-1.1,1.1); plot(-1.1,-1.1);
xgrid(-1:0.5:1); ygrid(-1:0.5:1); hold on;
title("Click with the mouse. Finish -> here <-");
t=zeros(1,0); s=t; n=0;
repeat;
m=mouse(); if abs(m(1))>1.1 || abs(m(2))>1.1; break; endif;
t=t|m(1); s=s|m(2);
n=n+1;
mark([m(1),m(1)],[m(2),m(2)]);
end;
{t,i}=sort(t); s=s(i);
sp=spline(t,s); x=linspace(min(t),max(t),n*10);
plot(x,splineval(t,s,sp,x));
color(2); plot(x,interpval(t,interp(t,s),x)); color(1);
wait(180);
hold off;
return 0;
endfunction
function test2
## Demonstrates the spline functions in RMAT.
"Click any number of points (x,y)"
"The program will compute the interpolating spline"
"and the interpolating polynomial in each variable."
"Press any key, please!"
wait(180);
""
setplot(-1.1,1.1,-1.1,1.1); plot(-1.1,-1.1);
xgrid(-1:0.5:1); ygrid(-1:0.5:1); hold on;
title("Click with the mouse. Finish -> here <-");
t=zeros(1,0); s=t; n=0;
repeat;
m=mouse(); if abs(m(1))>1.1 || abs(m(2))>1.1; break; endif;
t=t|m(1); s=s|m(2);
n=n+1;
mark([m(1),m(1)],[m(2),m(2)]);
end;
spt=spline(1:n,t); sps=spline(1:n,s);
x=linspace(1,n,600);
plot(splineval(1:n,t,spt,x),splineval(1:n,s,sps,x));
color(2);
plot(interpval(1:n,interp(1:n,t),x),interpval(1:n,interp(1:n,s),x));
color(1);
wait(180);
hold off;
return 0;
endfunction
test(); test2();
|
