summaryrefslogtreecommitdiff
path: root/26/CH6/EX6.4.1/6_4_1.sce
blob: beeb5b0e0b0dbabb919364bad38d50ae854a9612 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
disp('to orthogonalise the given vectors using Gram-Schimdt orthogonalisation')
x1=[3;0;-1]
x2=[8;5;-6]
disp(x2,'x2=',x1,'x1=')
disp('Let v1=x1')
v1=x1
disp('v2=x2-((x2.v1)/(v1.v1))*v1')
a1=x2'*v1
a2=v1'*v1
p=(a1/a2)*v1
v2=x2-p
disp(p,'-',x2,'=')
disp(v2,'=')
disp('Thus, an orthogonal basis is:')
disp(v2,v1)