|
|
1.1 ! root 1: real u(3), d(3,3), xyz(3,5), work(25) ! 2: data u / .5, .7, 0 / ! 3: data d / 0., 0., 0., ! 4: $ 1., 0., 1., ! 5: $ 0., 1., 10. / ! 6: data xyz / 0., 0., 0., ! 7: $ 1., 0., 0., ! 8: $ 0., 1., 10., ! 9: $ 0., 2., 20.5, ! 10: $ 0., 3., 29.4 / ! 11: call interp(u,d(1,1), d(1,2),d(1,3)) ! 12: write(6,1001)u(3) ! 13: 1001 format(' interp:',f8.3) ! 14: c call l2fit(5,xyz,u,work) ! 15: c write(6,1002)u ! 16: 1002 format(' l2fit:',3f8.3) ! 17: end
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.