fork download
copy 
  1. format longG;
  2. pkg load optim;
  3. vx=[1,4,7,10,12,15,16,18,22,23];
  4. vy=[9.3,23.4,34.7,45.2,52,61.9,65.2,71.8,84.7,87.9];
  5. vi=[1.1,1.2,1.3];
  6. fa=@(c,x)c(1)*x+c(2)*log10(x)+c(3)*nthroot(x,5);
  7. c=nlinfit(vx,vy,fa,vi)
  8. f=@(x)fa(c,x);
  9.  x=1:0.01:23;
  10. plot(vx,vy,".r",x,f(x),"-b");
  11. round(rref([1,4,6,2,f(2); 2,3,1,5,f(3); 3,1,4,7,f(4); 4,2,5,6,f(5)])(:,5)'*1e2)/1e2
Success #stdin #stdout 0.4s 65176KB
comments (?)
stdin
copy 
1
2
88
42
10
stdout
copy 
c =

  3.003159895798863
  5.176803672003957
  6.291133162616028

ans =

  3.4  1.61  0.33  1.47
https://ideone.com/3C6Epy
language:
Octave (octave 4.4.1)
created:
14 days ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!