Classic Fourier Series

[salvomic]

(02-01-2015 05:12 PM)Snorre Wrote:  another approach:

Code:

#cas
...
  f:=argv(1);
  k:=argv(argc);
  IF argc=3 THEN f:=zip('unapply',f,argv(2)); END;
 ...
#end

This doesn't check types, but number of args, so usage is either fourier(expr,var,k) or fourier(func,k).

I like this approach!

please, Snorre, help with this code (see thread with Han in http://www.hpmuseum.org/forum/thread-2992.html...
I would like calculate curvilinear integral with input a function (x,y,z), a parametric curve [r(t), r(t), r(t)], low bound, high bound; then I would write also a program for line integral (vectorial function, so: [x,y,z] fun, [t,t,t] param, l, h)...
I need some controls with the expression or function (for now if I put f and not f(x) Prime reset, as we now), for the parametric vector [t,...,...] e mostly for the variables in those two: they could be 2 or 3...
Any help?

With expression in x,y,z, vector [t...],l,h the function works well, the only control now is that there be 4 parameters...

Code:


#cas
intcur(args):=
BEGIN
local a, b, f, r, dr, ft;

argv:=[args];
argc:=size(argv);
IF argc !=4 THEN
return "Input: scalar func, [param curve t] ,low, high"; 
 ELSE
f:=argv(1);
r:=argv(2);
a:=argv(3);
b:=argv(4);
dr:=diff(r,t);
ft:=subst(f,[x,y,z]=r);
return int(dot(ft,l2norm(dr)),t,a,b);
END;

END;
#end

I hope don't bore you long time

P.S. we could also continue in the other thread, if you want, as here we are a bit off topic