uiopen('/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/BScurve0.m', true);
knots=[0,0,0,1,2,3,3,3];
t=linspace(0,1,10);
k=3;
B = BSplineEval(k,knots,t);
B
B =
Columns 1 through 3
[1x10 double] [1x10 double] [1x10 double]
Columns 4 through 5
[1x10 double] [1x10 double]
B{:}
ans =
Columns 1 through 5
1.0000 0.7901 0.6049 0.4444 0.3086
Columns 6 through 10
0.1975 0.1111 0.0494 0.0123 0
ans =
Columns 1 through 5
0 0.2037 0.3704 0.5000 0.5926
Columns 6 through 10
0.6481 0.6667 0.6481 0.5926 0.5000
ans =
Columns 1 through 5
0 0.0062 0.0247 0.0556 0.0988
Columns 6 through 10
0.1543 0.2222 0.3025 0.3951 0.5000
ans =
Columns 1 through 8
0 0 0 0 0 0 0 0
Columns 9 through 10
0 0
ans =
Columns 1 through 8
0 0 0 0 0 0 0 0
Columns 9 through 10
0 0
cat(1,B{:})
ans =
Columns 1 through 5
1.0000 0.7901 0.6049 0.4444 0.3086
0 0.2037 0.3704 0.5000 0.5926
0 0.0062 0.0247 0.0556 0.0988
0 0 0 0 0
0 0 0 0 0
Columns 6 through 10
0.1975 0.1111 0.0494 0.0123 0
0.6481 0.6667 0.6481 0.5926 0.5000
0.1543 0.2222 0.3025 0.3951 0.5000
0 0 0 0 0
0 0 0 0 0
cat(2,B{:})
ans =
Columns 1 through 5
1.0000 0.7901 0.6049 0.4444 0.3086
Columns 6 through 10
0.1975 0.1111 0.0494 0.0123 0
Columns 11 through 15
0 0.2037 0.3704 0.5000 0.5926
Columns 16 through 20
0.6481 0.6667 0.6481 0.5926 0.5000
Columns 21 through 25
0 0.0062 0.0247 0.0556 0.0988
Columns 26 through 30
0.1543 0.2222 0.3025 0.3951 0.5000
Columns 31 through 35
0 0 0 0 0
Columns 36 through 40
0 0 0 0 0
Columns 41 through 45
0 0 0 0 0
Columns 46 through 50
0 0 0 0 0
cat(3,B{:})
ans(:,:,1) =
Columns 1 through 5
1.0000 0.7901 0.6049 0.4444 0.3086
Columns 6 through 10
0.1975 0.1111 0.0494 0.0123 0
ans(:,:,2) =
Columns 1 through 5
0 0.2037 0.3704 0.5000 0.5926
Columns 6 through 10
0.6481 0.6667 0.6481 0.5926 0.5000
ans(:,:,3) =
Columns 1 through 5
0 0.0062 0.0247 0.0556 0.0988
Columns 6 through 10
0.1543 0.2222 0.3025 0.3951 0.5000
ans(:,:,4) =
Columns 1 through 8
0 0 0 0 0 0 0 0
Columns 9 through 10
0 0
ans(:,:,5) =
Columns 1 through 8
0 0 0 0 0 0 0 0
Columns 9 through 10
0 0
cat(1,B{:})
ans =
Columns 1 through 5
1.0000 0.7901 0.6049 0.4444 0.3086
0 0.2037 0.3704 0.5000 0.5926
0 0.0062 0.0247 0.0556 0.0988
0 0 0 0 0
0 0 0 0 0
Columns 6 through 10
0.1975 0.1111 0.0494 0.0123 0
0.6481 0.6667 0.6481 0.5926 0.5000
0.1543 0.2222 0.3025 0.3951 0.5000
0 0 0 0 0
0 0 0 0 0
B
B =
Columns 1 through 3
[1x10 double] [1x10 double] [1x10 double]
Columns 4 through 5
[1x10 double] [1x10 double]
B{1:3}
ans =
Columns 1 through 5
1.0000 0.7901 0.6049 0.4444 0.3086
Columns 6 through 10
0.1975 0.1111 0.0494 0.0123 0
ans =
Columns 1 through 5
0 0.2037 0.3704 0.5000 0.5926
Columns 6 through 10
0.6481 0.6667 0.6481 0.5926 0.5000
ans =
Columns 1 through 5
0 0.0062 0.0247 0.0556 0.0988
Columns 6 through 10
0.1543 0.2222 0.3025 0.3951 0.5000
B
B =
Columns 1 through 3
[1x10 double] [1x10 double] [1x10 double]
Columns 4 through 5
[1x10 double] [1x10 double]
B(1:3)
ans =
[1x10 double] [1x10 double] [1x10 double]
B=B(1:3)
B =
[1x10 double] [1x10 double] [1x10 double]
b
Warning: Could not find an exact (case-sensitive)
match for 'b'.
/home/accounts/personale/clrmrc90/aa1213/advanced_numerical_analysis/B.m
is a case-insensitive match and will be used
instead.
You can improve the performance of your code by
using exact
name matches and we therefore recommend that you
update your
usage accordingly. Alternatively, you can disable
this warning using
warning('off','MATLAB:dispatcher:InexactCaseMatch').
This warning will become an error in future
releases.
??? Input argument "n" is undefined.
Error in ==> B at 8
knots = linspace(0,1,n+1-k+2);
b=rand(2,6);
t=linspace(0,1);
S=BScurve3(b,k,t);
S1=BScurve4(b,k,t);
plot(S(1,:),S(2,:),S1(1,:),S1(2,:),'o')
b=rand(2,1000);
t=linspace(0,1,200);
tic,BScurve3(b,k,t);,toc
Elapsed time is 0.041484 seconds.
tic,BScurve3(b,k,t);,toc
Elapsed time is 0.036293 seconds.
tic,BScurve3(b,k,t);,toc
Elapsed time is 0.037249 seconds.
tic,BScurve3(b,k,t);,toc
Elapsed time is 0.036905 seconds.
tic,BScurve4(b,k,t);,toc
Elapsed time is 0.005303 seconds.
tic,BScurve4(b,k,t);,toc
Elapsed time is 0.001099 seconds.
tic,BScurve4(b,k,t);,toc
Elapsed time is 0.001053 seconds.
tic,BScurve3(b,k,t);,toc
Elapsed time is 0.037516 seconds.
tic,BScurve3(b,k,t);,toc
Elapsed time is 0.038075 seconds.
tic,BScurve3(b,k,t);,toc
Elapsed time is 0.037411 seconds.
edit spval
b=rand(2,10000);
t=linspace(0,1,1000);
tic,BScurve3(b,k,t);,toc
Elapsed time is 1.182929 seconds.
tic,BScurve3(b,k,t);,toc
Elapsed time is 1.118311 seconds.
tic,BScurve4(b,k,t);,toc
Elapsed time is 0.002148 seconds.
tic,BScurve4(b,k,t);,toc
Elapsed time is 0.002197 seconds.
tic,BScurve4(b,k,t);,toc
Elapsed time is 0.002197 seconds.
edit spmak
b=rand(2,5);
t=linspace(0,1);
[S,knots]=BScurve3(b,3,t);
knots
knots =
Columns 1 through 5
0 0 0 0.3333 0.6667
Columns 6 through 8
1.0000 1.0000 1.0000
tbar=rand
tbar =
0.1963
[b,knots]=knotinsertion(b,knots,tbar);
S1=BScurve3(b,3,t,knots);
plot(S(1,:),S(2,:),S1(1,:),S1(2,:),'o')
close all
plot(S(1,:),S(2,:),S1(1,:),S1(2,:),'o')
b=rand(2,10000);
t=linspace(0,1,1000);
BScurve4(b,3,t);
Elapsed time is 0.000973 seconds.
Elapsed time is 0.001295 seconds.
BScurve4(b,3,t);
Elapsed time is 0.000573 seconds.
Elapsed time is 0.001376 seconds.
BScurve3(b,3,t);
Elapsed time is 0.881757 seconds.
Elapsed time is 0.244868 seconds.
B{1:10}=zeros(2)
??? The right hand side of this assignment has too
few values to satisfy
the left hand side.
help deal
DEAL Deal inputs to outputs.
[A,B,C,...] = DEAL(X,Y,Z,...) simply matches up the input and
output lists. It is the same as A=X, B=Y, C=Z, ...
[A,B,C,...] = DEAL(X) copies the single input to all
the requested outputs. It is the same as A=X, B=X, C=X, ...
DEAL is most useful when used with cell arrays and structures
via comma separated list expansion. Here are some useful
constructions:
[S.FIELD] = DEAL(X) sets all the fields with the name FIELD
in the structure array S to the value X. If S doesn't
exist, use [S(1:M).FIELD] = DEAL(X);
[X{:}] = DEAL(A.FIELD) copies the values of the field with
name FIELD to the cell array X. If X doesn't exist,
use [X{1:M}] = DEAL(A.FIELD).
[A,B,C,...] = DEAL(X{:}) copies the contents of the cell
array X to the separate variables A,B,C,...
[A,B,C,...] = DEAL(S.FIELD) copies the contents of the fields
with the name FIELD to separate variables A,B,C,...
Examples:
sys = {rand(3) ones(3,1) eye(3) zeros(3,1)};
[a,b,c,d] = deal(sys{:});
direc = dir; filenames = {};
[filenames{1:length(direc),1}] = deal(direc.name);
See also lists, paren.
Reference page in Help browser
doc deal
B{1:10}=deal(zeros(2))
??? The right hand side of this assignment has too
few values to satisfy
the left hand side.
clear B
B{1:10}=deal(zeros(2))
??? The right hand side of this assignment has too
few values to satisfy
the left hand side.
[B{1:10}]=deal(zeros(2))
B =
Columns 1 through 3
[2x2 double] [2x2 double] [2x2 double]
Columns 4 through 6
[2x2 double] [2x2 double] [2x2 double]
Columns 7 through 9
[2x2 double] [2x2 double] [2x2 double]
Column 10
[2x2 double]
BScurve3(b,3,t);
??? Error using ==> zeros
Too many output arguments.
Error in ==> BSplineEval at 27
[b{1:n-1}] = zeros(size(xx));
Error in ==> BScurve3 at 13
B = BSplineEval(k,knots,t);
BScurve3(b,3,t);
Elapsed time is 1.311980 seconds.
Elapsed time is 0.240732 seconds.
BScurve3(b,3,t);
Elapsed time is 1.249103 seconds.
Elapsed time is 0.239237 seconds.
BScurve3(b,3,t);
Elapsed time is 0.870847 seconds.
Elapsed time is 0.242312 seconds.
BScurve5(b,3,t);
??? Error: File: BSplineEvalnew.m Line: 40 Column:
15
Unbalanced or unexpected parenthesis or bracket.
Error in ==> BScurve5 at 13
B = BSplineEvalnew(k,knots,t);
BScurve5(b,3,t);
??? Cell contents assignment to a non-cell array
object.
Error in ==> BSplineEvalnew at 29
b{i}(t(i) <= xx & xx < t(i+1)) = 1;
Error in ==> BScurve5 at 13
B = BSplineEvalnew(k,knots,t);
BScurve5(b,3,t);
Elapsed time is 1.675147 seconds.
Elapsed time is 0.024137 seconds.
BScurve5(b,3,t);
Elapsed time is 1.182775 seconds.
Elapsed time is 0.034674 seconds.
BScurve5(b,3,t);
Elapsed time is 1.187192 seconds.
Elapsed time is 0.020875 seconds.
BScurve5(b,3,t);
Elapsed time is 0.018988 seconds.
BScurve5(b,3,t);
Elapsed time is 1.187635 seconds.
Elapsed time is 0.047493 seconds.
BScurve4(b,3,t);
BScurve4(b,3,t);
ans =
2000 3
size(b)
ans =
2 10000
BScurve4(b,3,t);
ans =
2000 3
S=BScurve4(b,3,t);
ans =
2000 3
size(S)
ans =
2 1000
ver
-------------------------------------------------------------------------------------
MATLAB Version 7.6.0.324 (R2008a)
MATLAB License Number: 629186
Operating System: Linux 3.5.0-18-generic #29-Ubuntu SMP Fri Oct 19 10:26:51 UTC 2012 x86_64
Java VM Version: Java 1.6.0_01 with Sun Microsystems Inc. Java HotSpot(TM) 64-Bit Server VM mixed mode
-------------------------------------------------------------------------------------
MATLAB Version 7.6 (R2008a)
Simulink Version 7.1 (R2008a)
Bioinformatics Toolbox Version 3.1 (R2008a)
Control System Toolbox Version 8.1 (R2008a)
Image Processing Toolbox Version 6.1 (R2008a)
Neural Network Toolbox Version 6.0 (R2008a)
Optimization Toolbox Version 4.0 (R2008a)
Signal Processing Toolbox Version 6.9 (R2008a)
Simulink Control Design Version 2.3 (R2008a)
Spline Toolbox Version 3.3.4 (R2008a)
Statistics Toolbox Version 6.2 (R2008a)
Symbolic Math Toolbox Version 3.2.3 (R2008a)
b=rand(2,10);
t=linspace(0,1,5);
BScurve5
??? Input argument "b" is undefined.
Error in ==> BScurve5 at 7
n = size(b,2)-1;
BScurve5(b,3,t)
Elapsed time is 0.000946 seconds.
B =
1.0000 0 0 0 0
0 0 0 0 0
0 0.5000 0 0 0
0 0.5000 0 0 0
0 0 0.5000 0 0
0 0 0.5000 0 0
0 0 0 0.5000 0
0 0 0 0.5000 0
0 0 0 0 0
0 0 0 0 1.0000
Elapsed time is 0.000026 seconds.
ans =
0.8592 0.1532 0.7542 0.7387 0.7264
0.8721 0.1775 0.2715 0.8818 0.7949
t=linspace(0,1,10);
BScurve5(b,3,t)
Elapsed time is 0.000047 seconds.
B =
Columns 1 through 5
1.0000 0.0123 0 0 0
0 0.5926 0.0247 0 0
0 0.3951 0.6728 0.0556 0
0 0 0.3025 0.7222 0.0988
0 0 0 0.2222 0.7469
0 0 0 0 0.1543
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
Columns 6 through 10
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0.1543 0 0 0 0
0.7469 0.2222 0 0 0
0.0988 0.7222 0.3025 0 0
0 0.0556 0.6728 0.3951 0
0 0 0.0247 0.5926 0
0 0 0 0.0123 1.0000
Elapsed time is 0.000022 seconds.
ans =
Columns 1 through 5
0.8592 0.6414 0.1901 0.2466 0.6678
0.8721 0.4712 0.2192 0.1324 0.2275
Columns 6 through 10
0.7839 0.8103 0.7023 0.6276 0.7264
0.3571 0.7487 0.8656 0.3995 0.7949
BScurve5(b,3,t)
Elapsed time is 0.000048 seconds.
Elapsed time is 0.000019 seconds.
ans =
Columns 1 through 5
0.8592 0.6414 0.1901 0.2466 0.6678
0.8721 0.4712 0.2192 0.1324 0.2275
Columns 6 through 10
0.7839 0.8103 0.7023 0.6276 0.7264
0.3571 0.7487 0.8656 0.3995 0.7949
BScurve5(b,3,t)
Elapsed time is 0.000104 seconds.
ans =
26
Elapsed time is 0.000019 seconds.
ans =
Columns 1 through 5
0.8592 0.6414 0.1901 0.2466 0.6678
0.8721 0.4712 0.2192 0.1324 0.2275
Columns 6 through 10
0.7839 0.8103 0.7023 0.6276 0.7264
0.3571 0.7487 0.8656 0.3995 0.7949
b=rand(2,10000);
t=linspace(0,1,1000);
BScurve5(b,3,t)
Elapsed time is 1.183601 seconds.
ans =
2996
Elapsed time is 0.019988 seconds.
ans =
Columns 1 through 5
0.2600 0.6019 0.9531 0.4816 0.5863
0.4132 0.2405 0.7384 0.5623 0.4740
Columns 6 through 10
0.6667 0.1218 0.3230 0.7529 0.4785
0.6475 0.5752 0.8707 0.3972 0.0439
Columns 11 through 15
0.6813 0.3320 0.3124 0.4527 0.3199
0.4828 0.1780 0.4030 0.4912 0.8463
Columns 16 through 20
0.6393 0.3201 0.5325 0.3752 0.5384
0.1685 0.7007 0.4150 0.5873 0.6635
Columns 21 through 25
0.5708 0.3329 0.8243 0.3765 0.2796
0.5686 0.6628 0.3971 0.3449 0.7221
Columns 26 through 30
0.6045 0.6809 0.8759 0.6485 0.2848
0.8763 0.5824 0.6009 0.3245 0.6655
Columns 31 through 35
0.5574 0.7093 0.5375 0.5066 0.3236
0.5496 0.5513 0.4082 0.3514 0.1003
Columns 36 through 40
0.5488 0.4122 0.6803 0.4286 0.5695
0.2417 0.7349 0.2875 0.8939 0.6291
Columns 41 through 45
0.3667 0.5875 0.6765 0.2060 0.7490
0.4318 0.5023 0.3411 0.8155 0.5451
Columns 46 through 50
0.5611 0.2046 0.2750 0.5379 0.4285
0.8428 0.3458 0.6272 0.6708 0.7211
Columns 51 through 55
0.5162 0.2524 0.6686 0.5267 0.3539
0.2570 0.7107 0.5233 0.5769 0.2951
Columns 56 through 60
0.6461 0.7429 0.5738 0.5371 0.1824
0.1841 0.2721 0.4484 0.5099 0.2926
Columns 61 through 65
0.2310 0.7869 0.3929 0.3997 0.2975
0.2446 0.0841 0.5589 0.8374 0.3215
Columns 66 through 70
0.2523 0.6241 0.4616 0.4894 0.2006
0.2820 0.3901 0.6783 0.7439 0.8027
Columns 71 through 75
0.5552 0.6141 0.6942 0.4250 0.7982
0.3012 0.4258 0.3709 0.6591 0.7241
Columns 76 through 80
0.8915 0.3115 0.5201 0.5289 0.2648
0.9451 0.6191 0.1223 0.3283 0.4400
Columns 81 through 85
0.7230 0.4605 0.6937 0.7485 0.3614
0.2896 0.7195 0.4564 0.5263 0.5000
Columns 86 through 90
0.6082 0.1146 0.5522 0.1368 0.5808
0.2703 0.4550 0.6029 0.5166 0.4617
Columns 91 through 95
0.6791 0.3330 0.5058 0.8603 0.2351
0.3354 0.7343 0.0857 0.5498 0.4399
Columns 96 through 100
0.6093 0.8426 0.3150 0.4637 0.5934
0.3174 0.4532 0.3237 0.6303 0.7530
Columns 101 through 105
0.4836 0.9219 0.4478 0.6394 0.3699
0.4578 0.2948 0.5026 0.8390 0.3229
Columns 106 through 110
0.7712 0.4553 0.8519 0.4620 0.5077
0.3333 0.4922 0.3741 0.2114 0.2897
Columns 111 through 115
0.5952 0.3921 0.1242 0.5421 0.4279
0.4763 0.1932 0.7385 0.0387 0.6305
Columns 116 through 120
0.7205 0.1308 0.3371 0.8382 0.5564
0.2098 0.1415 0.7045 0.3593 0.3189
Columns 121 through 125
0.4285 0.6651 0.2712 0.7618 0.3293
0.3586 0.4809 0.8692 0.3273 0.4058
Columns 126 through 130
0.5966 0.7192 0.8078 0.7164 0.7330
0.7489 0.5108 0.5629 0.8364 0.5051
Columns 131 through 135
0.0976 0.4501 0.6478 0.2393 0.4875
0.6180 0.3927 0.6118 0.5820 0.0178
Columns 136 through 140
0.4473 0.1500 0.6719 0.6411 0.4287
0.3198 0.6510 0.5185 0.7443 0.7156
Columns 141 through 145
0.2511 0.2532 0.4251 0.4367 0.1008
0.3595 0.4538 0.4138 0.5146 0.2684
Columns 146 through 150
0.7660 0.6487 0.8631 0.0956 0.6468
0.5532 0.4401 0.1798 0.2110 0.6078
Columns 151 through 155
0.5694 0.3170 0.3500 0.1684 0.6068
0.5402 0.8404 0.6473 0.2268 0.2249
Columns 156 through 160
0.6802 0.4874 0.2369 0.6020 0.8361
0.4411 0.4424 0.6398 0.2622 0.3036
Columns 161 through 165
0.3762 0.8364 0.4750 0.3064 0.6058
0.5391 0.2464 0.1559 0.3465 0.4384
Columns 166 through 170
0.7455 0.7155 0.3864 0.6695 0.1897
0.7386 0.7099 0.4719 0.3075 0.1152
Columns 171 through 175
0.6848 0.2383 0.3492 0.5816 0.5471
0.7304 0.4684 0.6786 0.2133 0.3541
Columns 176 through 180
0.4158 0.4121 0.3003 0.2804 0.7341
0.8165 0.2568 0.9184 0.5049 0.2152
Columns 181 through 185
0.7019 0.7197 0.0548 0.6290 0.5590
0.5330 0.3248 0.1589 0.6750 0.4552
Columns 186 through 190
0.4469 0.2535 0.7272 0.4791 0.5489
0.1770 0.4069 0.6019 0.8903 0.4764
Columns 191 through 195
0.1828 0.4705 0.8581 0.5686 0.3458
0.6964 0.0793 0.2955 0.4909 0.9029
Columns 196 through 200
0.3944 0.3540 0.2088 0.4792 0.8062
0.7381 0.1831 0.6384 0.5720 0.8725
Columns 201 through 205
0.8148 0.6203 0.6815 0.8138 0.3680
0.6325 0.7350 0.2723 0.0693 0.6562
Columns 206 through 210
0.2552 0.3970 0.4202 0.6151 0.6993
0.0973 0.1080 0.8048 0.1442 0.5757
Columns 211 through 215
0.7820 0.7483 0.5608 0.7233 0.3858
0.4828 0.3062 0.8221 0.7687 0.2511
Columns 216 through 220
0.3276 0.5104 0.4929 0.8188 0.3415
0.9475 0.5523 0.6204 0.6650 0.6122
Columns 221 through 225
0.8105 0.4248 0.5192 0.4444 0.5562
0.2787 0.7208 0.3551 0.8118 0.8446
Columns 226 through 230
0.6554 0.4955 0.0854 0.3097 0.3397
0.2122 0.7118 0.0719 0.9535 0.8518
Columns 231 through 235
0.8151 0.5271 0.5816 0.3008 0.4298
0.5964 0.6231 0.4176 0.8744 0.8123
Columns 236 through 240
0.7324 0.7618 0.2103 0.4778 0.5779
0.7247 0.4063 0.3535 0.4490 0.3802
Columns 241 through 245
0.3605 0.2803 0.8803 0.2907 0.4903
0.5661 0.7568 0.4837 0.6479 0.5379
Columns 246 through 250
0.1065 0.2901 0.2087 0.5597 0.6460
0.3049 0.8089 0.7095 0.6934 0.5313
Columns 251 through 255
0.7591 0.5087 0.3595 0.9435 0.6754
0.4909 0.2219 0.4720 0.0637 0.8429
Columns 256 through 260
0.5184 0.6017 0.5003 0.4367 0.8282
0.4368 0.7066 0.2848 0.0896 0.6790
Columns 261 through 265
0.6910 0.2943 0.4467 0.4954 0.3341
0.3187 0.6755 0.3058 0.7179 0.7647
Columns 266 through 270
0.6967 0.4490 0.3695 0.5186 0.4841
0.7608 0.7879 0.7883 0.5261 0.3353
Columns 271 through 275
0.8643 0.3145 0.7677 0.8458 0.7370
0.4806 0.1839 0.7795 0.6301 0.7655
Columns 276 through 280
0.5499 0.1369 0.6057 0.4250 0.7514
0.6195 0.6790 0.9567 0.6501 0.2970
Columns 281 through 285
0.4572 0.8989 0.7594 0.3698 0.7996
0.5464 0.7274 0.6842 0.3549 0.0956
Columns 286 through 290
0.6544 0.4887 0.4960 0.6496 0.1393
0.4650 0.8369 0.2121 0.2631 0.5211
Columns 291 through 295
0.6821 0.6840 0.4608 0.4585 0.4722
0.7539 0.3668 0.2124 0.6183 0.1062
Columns 296 through 300
0.1860 0.4862 0.4462 0.6817 0.4113
0.8211 0.1729 0.4429 0.3016 0.2517
Columns 301 through 305
0.7324 0.6752 0.3911 0.5633 0.4928
0.8706 0.5757 0.4505 0.4648 0.7953
Columns 306 through 310
0.6059 0.5822 0.8184 0.3102 0.3788
0.4465 0.3281 0.3687 0.7185 0.3024
Columns 311 through 315
0.7909 0.6368 0.6586 0.4159 0.6999
0.1506 0.1927 0.2476 0.6492 0.4443
Columns 316 through 320
0.6712 0.5760 0.9605 0.5333 0.8412
0.2682 0.4652 0.8716 0.5113 0.1231
Columns 321 through 325
0.7367 0.5717 0.7737 0.2479 0.7275
0.0760 0.3735 0.7905 0.6650 0.5984
Columns 326 through 330
0.9171 0.1896 0.2203 0.1644 0.3927
0.1922 0.5572 0.5429 0.5324 0.6929
Columns 331 through 335
0.2648 0.2392 0.8332 0.3168 0.7587
0.4204 0.5781 0.8077 0.1888 0.7523
Columns 336 through 340
0.6622 0.1639 0.3598 0.6484 0.5605
0.2810 0.4972 0.3019 0.4031 0.9312
Columns 341 through 345
0.4860 0.3544 0.2103 0.5920 0.6133
0.5820 0.6292 0.5767 0.7289 0.4593
Columns 346 through 350
0.9537 0.0557 0.9914 0.0906 0.7013
0.3230 0.1080 0.3503 0.5312 0.3780
Columns 351 through 355
0.4691 0.1948 0.1507 0.5861 0.4895
0.3104 0.5612 0.2828 0.5638 0.3253
Columns 356 through 360
0.0648 0.6829 0.5363 0.7437 0.4401
0.2787 0.2367 0.6557 0.7222 0.5685
Columns 361 through 365
0.4914 0.1146 0.3174 0.2757 0.3874
0.5455 0.4940 0.6272 0.6599 0.3657
Columns 366 through 370
0.8695 0.2773 0.7271 0.6274 0.3148
0.8474 0.3999 0.4255 0.7266 0.4129
Columns 371 through 375
0.7360 0.4703 0.2094 0.7441 0.8319
0.3755 0.3126 0.3191 0.3284 0.7661
Columns 376 through 380
0.1616 0.5053 0.0508 0.6397 0.2043
0.5872 0.2020 0.5628 0.0622 0.5375
Columns 381 through 385
0.7513 0.3833 0.4325 0.4236 0.4737
0.8418 0.7909 0.2352 0.0895 0.3356
Columns 386 through 390
0.1839 0.2987 0.6764 0.8329 0.1574
0.8494 0.7551 0.7938 0.4433 0.2572
Columns 391 through 395
0.4891 0.4985 0.0542 0.4011 0.3274
0.4387 0.2010 0.2077 0.3901 0.7350
Columns 396 through 400
0.7928 0.8435 0.2579 0.4457 0.2735
0.4607 0.8416 0.4400 0.2108 0.4611
Columns 401 through 405
0.3098 0.2667 0.2679 0.4584 0.5086
0.4717 0.4402 0.6358 0.4794 0.7340
Columns 406 through 410
0.2402 0.7522 0.2448 0.3971 0.6253
0.6307 0.3841 0.6296 0.6263 0.8358
Columns 411 through 415
0.6694 0.3918 0.8161 0.3589 0.0791
0.6500 0.2885 0.8432 0.5288 0.6551
Columns 416 through 420
0.6309 0.6614 0.4583 0.5755 0.2716
0.8267 0.5318 0.2825 0.7720 0.6594
Columns 421 through 425
0.8514 0.7843 0.7870 0.8604 0.5997
0.2992 0.1280 0.8594 0.7335 0.6570
Columns 426 through 430
0.1581 0.4479 0.5977 0.1617 0.5544
0.9087 0.4569 0.8471 0.7002 0.5827
Columns 431 through 435
0.6198 0.5993 0.7282 0.2901 0.3148
0.5465 0.4334 0.4954 0.4550 0.5698
Columns 436 through 440
0.2040 0.3351 0.7042 0.9659 0.7357
0.2787 0.8716 0.6906 0.5862 0.3519
Columns 441 through 445
0.7986 0.7150 0.6752 0.8033 0.3962
0.4942 0.1626 0.6697 0.1668 0.8251
Columns 446 through 450
0.4301 0.2273 0.4326 0.6080 0.1836
0.1982 0.5278 0.6970 0.3051 0.6820
Columns 451 through 455
0.1226 0.7800 0.8835 0.3804 0.4992
0.2041 0.3223 0.1427 0.3401 0.5972
Columns 456 through 460
0.8105 0.1439 0.5491 0.4122 0.4150
0.7311 0.5985 0.6428 0.7076 0.2952
Columns 461 through 465
0.7184 0.6101 0.8108 0.4267 0.2824
0.3958 0.2041 0.2406 0.7544 0.5977
Columns 466 through 470
0.3738 0.6470 0.5964 0.7017 0.7271
0.6176 0.1245 0.6841 0.5670 0.2477
Columns 471 through 475
0.3035 0.3280 0.4126 0.2248 0.2898
0.3344 0.6879 0.0814 0.0724 0.5641
Columns 476 through 480
0.1777 0.2764 0.6031 0.3124 0.4046
0.4900 0.6071 0.3149 0.1525 0.4436
Columns 481 through 485
0.2971 0.6359 0.7196 0.8218 0.7823
0.7652 0.4587 0.3724 0.7845 0.1940
Columns 486 through 490
0.6383 0.8047 0.1819 0.7858 0.3866
0.2356 0.6502 0.8038 0.7431 0.2309
Columns 491 through 495
0.3341 0.4467 0.6912 0.5985 0.3857
0.7875 0.4060 0.6059 0.2867 0.5151
Columns 496 through 500
0.2718 0.2916 0.4308 0.2022 0.4991
0.5619 0.5912 0.5692 0.6258 0.1067
Columns 501 through 505
0.8775 0.5480 0.5536 0.2199 0.7077
0.6199 0.0577 0.4066 0.7976 0.3173
Columns 506 through 510
0.8035 0.4353 0.1380 0.7736 0.4802
0.3512 0.5683 0.3812 0.5476 0.4643
Columns 511 through 515
0.1203 0.8978 0.7144 0.6129 0.4346
0.6675 0.5859 0.7113 0.5549 0.4554
Columns 516 through 520
0.5701 0.3185 0.4403 0.2931 0.4988
0.4461 0.5662 0.8167 0.8816 0.5271
Columns 521 through 525
0.4429 0.3617 0.6830 0.7242 0.6453
0.1953 0.9057 0.1947 0.3597 0.4173
Columns 526 through 530
0.3144 0.6017 0.7120 0.3670 0.3599
0.4261 0.5272 0.4587 0.7710 0.4877
Columns 531 through 535
0.4153 0.2964 0.5801 0.3580 0.4265
0.6958 0.7512 0.6121 0.5095 0.6552
Columns 536 through 540
0.7301 0.2585 0.1891 0.3221 0.4669
0.4656 0.2024 0.0778 0.5766 0.8046
Columns 541 through 545
0.7899 0.4669 0.3779 0.3370 0.2232
0.2616 0.8879 0.1585 0.3186 0.6052
Columns 546 through 550
0.7165 0.3713 0.0636 0.7529 0.3048
0.4933 0.3192 0.3008 0.3930 0.6114
Columns 551 through 555
0.3899 0.6978 0.4620 0.7255 0.6610
0.7526 0.7263 0.6293 0.1868 0.4787
Columns 556 through 560
0.0535 0.3866 0.4840 0.4033 0.7136
0.8997 0.6985 0.0260 0.7923 0.6444
Columns 561 through 565
0.7093 0.2963 0.6369 0.4782 0.5932
0.7712 0.7481 0.4207 0.3922 0.7162
Columns 566 through 570
0.2320 0.2342 0.8383 0.4432 0.1037
0.7466 0.5091 0.8598 0.8633 0.3385
Columns 571 through 575
0.6612 0.1681 0.5953 0.4558 0.7737
0.8698 0.4644 0.6742 0.4511 0.1407
Columns 576 through 580
0.1151 0.8526 0.6434 0.7267 0.1843
0.5376 0.2195 0.5345 0.3573 0.7900
Columns 581 through 585
0.6259 0.1082 0.2307 0.5554 0.3206
0.7716 0.1920 0.4979 0.7072 0.2894
Columns 586 through 590
0.2149 0.5837 0.5399 0.5589 0.5789
0.1932 0.5326 0.3444 0.7377 0.1620
Columns 591 through 595
0.1445 0.5968 0.8500 0.4591 0.3436
0.7016 0.3915 0.7076 0.3648 0.4862
Columns 596 through 600
0.1654 0.4919 0.6756 0.7639 0.3224
0.4488 0.3389 0.4518 0.3458 0.4623
Columns 601 through 605
0.5851 0.5414 0.3685 0.1744 0.6167
0.8704 0.3127 0.3327 0.4967 0.4049
Columns 606 through 610
0.7986 0.1536 0.3511 0.2237 0.4367
0.6456 0.3241 0.7017 0.6777 0.1727
Columns 611 through 615
0.5932 0.4941 0.7913 0.6981 0.6629
0.4743 0.5246 0.3138 0.7435 0.5320
Columns 616 through 620
0.5593 0.8486 0.6547 0.1178 0.1843
0.1655 0.2044 0.8552 0.4917 0.7672
Columns 621 through 625
0.5692 0.0247 0.8493 0.4528 0.4676
0.7142 0.3582 0.8295 0.6874 0.8154
Columns 626 through 630
0.4653 0.2984 0.4113 0.5712 0.3279
0.6570 0.6752 0.5792 0.3830 0.6490
Columns 631 through 635
0.4178 0.4315 0.3723 0.5496 0.7802
0.5058 0.4160 0.8040 0.3385 0.4809
Columns 636 through 640
0.4985 0.1743 0.5805 0.6948 0.5225
0.6021 0.3791 0.5777 0.4667 0.5705
Columns 641 through 645
0.3129 0.4838 0.7670 0.3355 0.6413
0.8189 0.8538 0.7639 0.3436 0.5600
Columns 646 through 650
0.5249 0.3112 0.9358 0.6617 0.8534
0.4295 0.7625 0.7454 0.3413 0.3757
Columns 651 through 655
0.3658 0.1233 0.5728 0.1834 0.1827
0.7178 0.8738 0.5244 0.6930 0.6020
Columns 656 through 660
0.1168 0.4941 0.3020 0.5523 0.6382
0.6663 0.3562 0.4633 0.1807 0.7358
Columns 661 through 665
0.1896 0.6596 0.1707 0.3566 0.0485
0.4446 0.2977 0.9003 0.2528 0.4812
Columns 666 through 670
0.1432 0.5265 0.5205 0.1394 0.1443
0.1400 0.1389 0.1454 0.5736 0.8471
Columns 671 through 675
0.5767 0.5236 0.9326 0.6752 0.8596
0.2970 0.6564 0.7635 0.1263 0.6123
Columns 676 through 680
0.1949 0.4023 0.1753 0.4524 0.1331
0.5481 0.3332 0.2890 0.6087 0.1017
Columns 681 through 685
0.8027 0.2251 0.6579 0.7290 0.4883
0.5156 0.5948 0.5508 0.6734 0.8319
Columns 686 through 690
0.4143 0.5323 0.6114 0.8141 0.2205
0.3149 0.4367 0.2729 0.4342 0.6377
Columns 691 through 695
0.7731 0.4151 0.3011 0.4946 0.4825
0.8304 0.4280 0.3437 0.4999 0.5430
Columns 696 through 700
0.3813 0.0746 0.3994 0.6466 0.3045
0.7645 0.4262 0.5049 0.4951 0.4954
Columns 701 through 705
0.5401 0.2582 0.2319 0.3680 0.2775
0.4858 0.2415 0.2025 0.3431 0.7094
Columns 706 through 710
0.4658 0.5178 0.7778 0.3433 0.7821
0.7372 0.5496 0.8461 0.7161 0.4072
Columns 711 through 715
0.3814 0.6678 0.6325 0.4752 0.2795
0.7180 0.1282 0.4024 0.8392 0.3255
Columns 716 through 720
0.4369 0.1509 0.5858 0.3429 0.2212
0.3445 0.4942 0.9328 0.7161 0.8775
Columns 721 through 725
0.5218 0.2718 0.8409 0.1878 0.1769
0.5310 0.5649 0.2568 0.7351 0.1723
Columns 726 through 730
0.5251 0.5201 0.3858 0.2517 0.6931
0.2127 0.7260 0.3830 0.7443 0.4243
Columns 731 through 735
0.6795 0.1675 0.3783 0.8938 0.8479
0.3846 0.6038 0.2005 0.1419 0.2107
Columns 736 through 740
0.3915 0.4077 0.2363 0.1239 0.5380
0.2600 0.5617 0.7302 0.4761 0.2913
Columns 741 through 745
0.6981 0.7681 0.8923 0.6495 0.3616
0.3246 0.3115 0.1457 0.4903 0.5943
Columns 746 through 750
0.4943 0.3206 0.4267 0.2114 0.7099
0.8548 0.7934 0.4770 0.0744 0.4449
Columns 751 through 755
0.6722 0.6833 0.4368 0.6841 0.3264
0.4812 0.1995 0.7434 0.4582 0.1745
Columns 756 through 760
0.6413 0.2322 0.6799 0.4183 0.2781
0.3815 0.3041 0.0954 0.2498 0.7133
Columns 761 through 765
0.7104 0.5398 0.2912 0.3841 0.2032
0.3635 0.5693 0.1892 0.7926 0.6736
Columns 766 through 770
0.3022 0.0983 0.7841 0.6640 0.1754
0.1243 0.5713 0.8084 0.9700 0.5617
Columns 771 through 775
0.8234 0.6500 0.7166 0.5418 0.6048
0.4772 0.7576 0.1156 0.6121 0.1458
Columns 776 through 780
0.3700 0.2920 0.7916 0.2778 0.2076
0.7351 0.5873 0.6436 0.1472 0.4579
Columns 781 through 785
0.7858 0.4844 0.6357 0.9427 0.7915
0.8487 0.3282 0.5009 0.6387 0.3906
Columns 786 through 790
0.1519 0.6212 0.3144 0.0841 0.6433
0.5745 0.6996 0.6213 0.4777 0.5000
Columns 791 through 795
0.8218 0.2469 0.7945 0.5891 0.5612
0.3441 0.6439 0.2443 0.7634 0.9110
Columns 796 through 800
0.4343 0.7348 0.8386 0.2239 0.2264
0.5508 0.1695 0.6783 0.1811 0.8020
Columns 801 through 805
0.4230 0.1972 0.6629 0.6044 0.1960
0.1048 0.1148 0.5667 0.5182 0.5039
Columns 806 through 810
0.2295 0.3145 0.7677 0.5318 0.3328
0.6847 0.3756 0.7635 0.3193 0.7255
Columns 811 through 815
0.2733 0.3031 0.7043 0.4605 0.3088
0.2367 0.9225 0.1963 0.6632 0.6224
Columns 816 through 820
0.5447 0.4891 0.1597 0.7263 0.8942
0.5938 0.8122 0.3303 0.7698 0.8113
Columns 821 through 825
0.7342 0.4298 0.5436 0.6633 0.4334
0.1019 0.2322 0.8354 0.9084 0.3955
Columns 826 through 830
0.4736 0.1217 0.7783 0.1990 0.6995
0.7120 0.7976 0.4064 0.9203 0.6577
Columns 831 through 835
0.6823 0.5417 0.6016 0.6938 0.2249
0.7754 0.6983 0.4968 0.4928 0.7700
Columns 836 through 840
0.2981 0.6670 0.2229 0.0954 0.1762
0.5173 0.5546 0.1131 0.5578 0.5731
Columns 841 through 845
0.5291 0.6236 0.2198 0.6275 0.9588
0.3629 0.5678 0.2846 0.7525 0.5821
Columns 846 through 850
0.4247 0.0986 0.6633 0.2652 0.9224
0.6269 0.4207 0.6095 0.6632 0.3609
Columns 851 through 855
0.3108 0.7143 0.3269 0.6552 0.5442
0.8312 0.7484 0.3860 0.5655 0.6276
Columns 856 through 860
0.4111 0.5343 0.2115 0.3148 0.4653
0.1057 0.3320 0.5199 0.4799 0.6410
Columns 861 through 865
0.3757 0.6253 0.5973 0.1590 0.3438
0.3429 0.7516 0.2425 0.1407 0.3410
Columns 866 through 870
0.7040 0.5153 0.4344 0.4271 0.5703
0.5050 0.6962 0.5691 0.1145 0.2794
Columns 871 through 875
0.3310 0.3064 0.4744 0.7880 0.4672
0.2742 0.2486 0.6261 0.8983 0.5702
Columns 876 through 880
0.5824 0.3804 0.5682 0.4670 0.6115
0.6523 0.6363 0.2759 0.4898 0.3189
Columns 881 through 885
0.8819 0.4377 0.5148 0.5657 0.4814
0.3516 0.6314 0.3933 0.3726 0.6196
Columns 886 through 890
0.6332 0.6617 0.1075 0.2913 0.4653
0.5127 0.6737 0.8641 0.5291 0.4369
Columns 891 through 895
0.2954 0.7567 0.7697 0.3572 0.4000
0.1064 0.3857 0.3356 0.5043 0.2421
Columns 896 through 900
0.8930 0.5833 0.8304 0.8143 0.4788
0.4766 0.6799 0.6139 0.1985 0.4208
Columns 901 through 905
0.6084 0.3901 0.5416 0.5193 0.0651
0.3543 0.1481 0.6747 0.3865 0.7525
Columns 906 through 910
0.3630 0.3375 0.4274 0.4122 0.2279
0.5480 0.4014 0.6240 0.5814 0.8423
Columns 911 through 915
0.8177 0.4452 0.7498 0.2261 0.4054
0.2449 0.2521 0.5970 0.3474 0.6270
Columns 916 through 920
0.4903 0.2173 0.8216 0.7182 0.5862
0.5524 0.3411 0.3500 0.7275 0.4353
Columns 921 through 925
0.4380 0.5331 0.5605 0.2603 0.3498
0.8154 0.8956 0.2057 0.1347 0.2416
Columns 926 through 930
0.1102 0.5829 0.6839 0.3747 0.8361
0.7240 0.3339 0.3886 0.7175 0.1478
Columns 931 through 935
0.8715 0.2694 0.5158 0.8499 0.7445
0.6570 0.6335 0.7055 0.2765 0.7741
Columns 936 through 940
0.4881 0.6817 0.7236 0.4331 0.1186
0.6634 0.8444 0.9275 0.6527 0.3933
Columns 941 through 945
0.4473 0.6161 0.6106 0.5283 0.6207
0.5182 0.4842 0.4450 0.8084 0.3218
Columns 946 through 950
0.8492 0.7702 0.8613 0.3274 0.5851
0.1653 0.2425 0.4486 0.2963 0.3659
Columns 951 through 955
0.8359 0.8664 0.1128 0.6381 0.6829
0.6464 0.5308 0.2527 0.5682 0.5138
Columns 956 through 960
0.5528 0.8021 0.3645 0.5003 0.6390
0.4375 0.5552 0.4194 0.2630 0.8296
Columns 961 through 965
0.4225 0.3755 0.9810 0.9015 0.4628
0.5847 0.3148 0.4262 0.7340 0.4772
Columns 966 through 970
0.1962 0.4153 0.5379 0.4699 0.4218
0.6407 0.3685 0.6219 0.4272 0.2935
Columns 971 through 975
0.6881 0.9468 0.6021 0.9146 0.3316
0.7233 0.5224 0.8638 0.7510 0.3128
Columns 976 through 980
0.5347 0.6869 0.6987 0.9372 0.7292
0.5122 0.8979 0.7669 0.5111 0.0440
Columns 981 through 985
0.3305 0.2251 0.2373 0.6218 0.6149
0.3432 0.6770 0.5224 0.2715 0.4119
Columns 986 through 990
0.2986 0.5948 0.1061 0.7921 0.5757
0.2627 0.2338 0.4309 0.6582 0.3230
Columns 991 through 995
0.4448 0.4116 0.2003 0.4020 0.6271
0.5988 0.3530 0.5566 0.4122 0.7376
Columns 996 through 1000
0.4904 0.4341 0.7593 0.6748 0.8192
0.5020 0.6067 0.5920 0.8614 0.5463
BScurve5(b,3,t);
Elapsed time is 1.187500 seconds.
ans =
2996
Elapsed time is 0.025372 seconds.
t=[0,0,0,1,2,3,4,5,6,6,6]
t =
Columns 1 through 8
0 0 0 1 2 3 4 5
Columns 9 through 11
6 6 6
tbar = 3.5
tbar =
3.5000
find(knots > tbar)
ans =
Empty matrix: 1-by-0
find(t > tbar)
ans =
7 8 9 10 11
t
t =
Columns 1 through 8
0 0 0 1 2 3 4 5
Columns 9 through 11
6 6 6
find(t > tbar,1)
ans =
7
find(t > tbar,1)-1
ans =
6
tbar=6
tbar =
6
find(t > tbar,1)-1
ans =
Empty matrix: 1-by-0
xy=[0,1,2;10,11,12]
xy =
0 1 2
10 11 12
xx=[0,0.5,1,1.5,2];
yy=[10,10.5,11,11.5,12];
[X,Y]=ndgrid(xx,yy);
X
X =
0 0 0 0 0
0.5000 0.5000 0.5000 0.5000 0.5000
1.0000 1.0000 1.0000 1.0000 1.0000
1.5000 1.5000 1.5000 1.5000 1.5000
2.0000 2.0000 2.0000 2.0000 2.0000
Y
Y =
10.0000 10.5000 11.0000 11.5000 12.0000
10.0000 10.5000 11.0000 11.5000 12.0000
10.0000 10.5000 11.0000 11.5000 12.0000
10.0000 10.5000 11.0000 11.5000 12.0000
10.0000 10.5000 11.0000 11.5000 12.0000
XY=[X(:)';Y(:)']
XY =
Columns 1 through 5
0 0.5000 1.0000 1.5000 2.0000
10.0000 10.0000 10.0000 10.0000 10.0000
Columns 6 through 10
0 0.5000 1.0000 1.5000 2.0000
10.5000 10.5000 10.5000 10.5000 10.5000
Columns 11 through 15
0 0.5000 1.0000 1.5000 2.0000
11.0000 11.0000 11.0000 11.0000 11.0000
Columns 16 through 20
0 0.5000 1.0000 1.5000 2.0000
11.5000 11.5000 11.5000 11.5000 11.5000
Columns 21 through 25
0 0.5000 1.0000 1.5000 2.0000
12.0000 12.0000 12.0000 12.0000 12.0000
[xyd,XYd]=ndgrid(xy(1,:),XY(1,:))
xyd =
Columns 1 through 8
0 0 0 0 0 0 0 0
1 1 1 1 1 1 1 1
2 2 2 2 2 2 2 2
Columns 9 through 16
0 0 0 0 0 0 0 0
1 1 1 1 1 1 1 1
2 2 2 2 2 2 2 2
Columns 17 through 24
0 0 0 0 0 0 0 0
1 1 1 1 1 1 1 1
2 2 2 2 2 2 2 2
Column 25
0
1
2
XYd =
Columns 1 through 5
0 0.5000 1.0000 1.5000 2.0000
0 0.5000 1.0000 1.5000 2.0000
0 0.5000 1.0000 1.5000 2.0000
Columns 6 through 10
0 0.5000 1.0000 1.5000 2.0000
0 0.5000 1.0000 1.5000 2.0000
0 0.5000 1.0000 1.5000 2.0000
Columns 11 through 15
0 0.5000 1.0000 1.5000 2.0000
0 0.5000 1.0000 1.5000 2.0000
0 0.5000 1.0000 1.5000 2.0000
Columns 16 through 20
0 0.5000 1.0000 1.5000 2.0000
0 0.5000 1.0000 1.5000 2.0000
0 0.5000 1.0000 1.5000 2.0000
Columns 21 through 25
0 0.5000 1.0000 1.5000 2.0000
0 0.5000 1.0000 1.5000 2.0000
0 0.5000 1.0000 1.5000 2.0000
xyd-XYd
ans =
Columns 1 through 5
0 -0.5000 -1.0000 -1.5000 -2.0000
1.0000 0.5000 0 -0.5000 -1.0000
2.0000 1.5000 1.0000 0.5000 0
Columns 6 through 10
0 -0.5000 -1.0000 -1.5000 -2.0000
1.0000 0.5000 0 -0.5000 -1.0000
2.0000 1.5000 1.0000 0.5000 0
Columns 11 through 15
0 -0.5000 -1.0000 -1.5000 -2.0000
1.0000 0.5000 0 -0.5000 -1.0000
2.0000 1.5000 1.0000 0.5000 0
Columns 16 through 20
0 -0.5000 -1.0000 -1.5000 -2.0000
1.0000 0.5000 0 -0.5000 -1.0000
2.0000 1.5000 1.0000 0.5000 0
Columns 21 through 25
0 -0.5000 -1.0000 -1.5000 -2.0000
1.0000 0.5000 0 -0.5000 -1.0000
2.0000 1.5000 1.0000 0.5000 0
doc tpaps
tpss
diary off