x=@(t) 2*cos(t); y=@(t) sin(t); z=@(t) t;
Dx=@(t) -2*sin(t); Dy=@(t) cos(t); Dz=@(t) 1+0*t;
T=@(x,y,z) x.^2+y.^2+z.^2;
N=50; tv=linspace(0,2*pi,N+1);
xv=x(tv); yv=y(tv); zv=z(tv);
Dxv=Dx(tv); Dyv=Dy(tv); Dzv=Dz(tv);
ntv=sqrt(Dxv.^2+Dyv.^2+Dzv.^2);
mitjanaaprox=intaprox/longaprox;
fprintf('mitjanaaprox = %.6f\n',mitjanaaprox)
x=@(t) 3*t./(1+t.^3); y=@(t) 3*t.^2./(1+t.^3);
Dx=@(t) 3*(1-2*t.^3)./(1+t.^3).^2; Dy=@(t) 3*(2*t-t.^4)./(1+t.^3).^2;
nt=@(t) sqrt(Dx(t).^2+Dy(t).^2);
g=@(t) rho(x(t),y(t)).*nt(t);
fprintf('massa = %.6f\n',massa)
x=@(t) 3*cos(t); y=@(t) 2*sin(t);
Dx=@(t) -3*sin(t); Dy=@(t) 2*cos(t);
N=20; tv=linspace(0,2*pi,N+1);
fprintf('intaprox = %.6f\n',intaprox)
x=@(t) cos(t.^2); y=@(t) sin(t.^2); z=@(t) t.^2;
Dx=@(t) -2*t.*sin(t.^2); Dy=@(t) 2*t.*cos(t.^2); Dz=@(t) 2*t;
R=@(x,y,z) 3*x-cos(2*z.^2);
xv=x(tv); yv=y(tv); zv=z(tv);
Dxv=Dx(tv); Dyv=Dy(tv); Dzv=Dz(tv);
gv=Pv.*Dxv+Qv.*Dyv+Rv.*Dzv;
intaprox(i)=trapz(tv,gv);
Delta=intaprox(2)-intaprox(1);
fprintf('Delta = %.6e\n',Delta)
x=@(t) t.^2; y=@(t) 1./(1+t); z=@(t) cos(t);
Dx=@(t) 2*t; Dy=@(t) -1./(1+t).^2; Dz=@(t) -sin(t);
P=@(t,a) 1+a*x(t).*y(t).*z(t);
Q=@(t,a) 1./(a+x(t).*y(t)+x(t).*z(t)+y(t).*z(t));
R=@(t,a) exp(a*x(t).*y(t).*z(t));
g1=@(t) P(t,a).*Dx(t)+Q(t,a).*Dy(t)+R(t,a).*Dz(t);
g2=@(t) P(t,a).*Dx(t)+Q(t,a).*Dy(t)+R(t,a).*Dz(t);
fprintf('Delta = %.6e\n',Delta)