aboutsummaryrefslogtreecommitdiff
path: root/math/llsq/progs/prog4.f
diff options
context:
space:
mode:
authorJoseph Hunkeler <jhunkeler@gmail.com>2015-07-08 20:46:52 -0400
committerJoseph Hunkeler <jhunkeler@gmail.com>2015-07-08 20:46:52 -0400
commitfa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 (patch)
treebdda434976bc09c864f2e4fa6f16ba1952b1e555 /math/llsq/progs/prog4.f
downloadiraf-linux-fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4.tar.gz
Initial commit
Diffstat (limited to 'math/llsq/progs/prog4.f')
-rw-r--r--math/llsq/progs/prog4.f22
1 files changed, 22 insertions, 0 deletions
diff --git a/math/llsq/progs/prog4.f b/math/llsq/progs/prog4.f
new file mode 100644
index 00000000..090b1865
--- /dev/null
+++ b/math/llsq/progs/prog4.f
@@ -0,0 +1,22 @@
+c prog4
+c c.l.lawson and r.j.hanson, jet propulsion laboratory, 1973 jun 12
+c to appear in 'solving least squares problems', prentice-hall, 1974
+c demonstrate singular value analysis.
+c
+ dimension a(15,5),b(15),sing(15)
+ data names/1h /
+c
+ read (5,10) ((a(i,j),j=1,5),b(i),i=1,15)
+ write (6,20)
+ write (6,30) ((a(i,j),j=1,5),b(i),i=1,15)
+ write (6,40)
+c
+ call sva (a,15,15,5,15,b,sing,names,1,d)
+c
+ stop
+ 10 format (6f12.0)
+ 20 format (46h1prog4. demonstrate singular value analysis/53h list
+ 1ing of input matrix, a, and vector, b, follows..)
+ 30 format (1h /(5f12.8,f20.4))
+ 40 format (1h1)
+ end