🚀 go-pugleaf

hi

i want fitting a scattered data whit spline.
I have mathlab and i had used a nag library (E02DDF - E02DFF), but i have
problems.
I don't know, if the problem is in my source or in the selcted nag
library. 

Could you check my code, or propouse the new method?
Thanking you in advance.

Please use for reply my e-mail: qnjose@correo.usc.es
Jose A. Martinez Pontevedra
Universidad de Santiago de Compostela
Dpto de Quimica Analtica, Nutricion y Bromatologia

---------------------------------------------------
%Coordenada X
x=[0.0000 3.6742 5.5114 1.8371 1.8371 0.0000 3.6742 5.5114 1.8371 7.3485
...
   3.6742 0.0000 0.0000 1.8371 0.0000 4.9235 4.9235 2.4250 2.4250 0.0000
...
   0.0000 4.8500 1.2492 1.2492 2.425]';

%Coordenada Y
y=[0.5774 0.0000 -0.5774 -0.5774 -0.1925 -0.5774 0.5774 0.5774 0.1925
0.0000 ...
      -0.5774 0.2887 -0.2887 0.5774 0.0000 0.5774 -0.5774 0.5774 -0.5774
0.1963 ...
      -0.1963 0.0000 0.3408 -0.3408 -0.0086]';

%Coordenada Z
f=[1.75 1.85 1.91 11.31 2.71 39.97 1.64 1.67 1.77 1.67 ...
      1.96 2.16 25.24 1.96 3.83 1.67 1.95 1.68 5.87 2.54 ...
      19.24 1.67 1.76 9.47 1.88]';

Puntos=25;
s=1;
w=ones(25,1);
lamda=zeros(12,1);
mu=zeros(12,1);

%evaluation of a fitted bicubic splines with automatic knot placemnent,
scattered data
[lamda,mu,c,fp,rank,wrk,iwrk,nx,ny,ifail] = e02ddf(x,y,f,s,lamda,mu);

%personal data

     x1=[0.00 1.00 2.00 3.00 4.00 5.00 7.00]';
     y1=[-0.60 -0.30 -0.10 0.10 0.30 0.60]';

lamda1=lamda';
mu1=mu';
c1=c';

%Evaluation of a fitted bicubic spline at a mesh of points
	[ff,ifail] = e02dff(x1,y1,lamda1,mu1,c1)

------------------------------------------------------------