1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
|
// =============================================================================
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) ????-2008 - INRIA
//
// This file is distributed under the same license as the Scilab package.
// =============================================================================
// Test 1
deff("[z]=f(x,y,t)","z=x^2+y^2+t");
res = feval(1:10,1:5,list(f,10));
res1 = [];
for j=1:5
res1=[res1, ((1:10).*(1:10))'+(j**2+10)*ones(10,1)];
end
if res<>res1 then bugmes();quit;end
// Test 2
deff("[z]=g(x,y)","z=x+%i*y");
res =feval(1:10,1:5,g);
res1 =[];
for j=1:5
res1=[res1, (1:10)'+%i*j*ones(10,1)];
end
if res<>res1 then bugmes();quit;end
// An other example from :
// http://numerics.tinabargs.com/2008/07/scilab-functions-as-input-for-other-scilab-functions/
// Thanks to Ma. Cristina R. Bargo for the authorization to include this test
// into Scilab
function y = plop(x)
y = 2*x - 1;
endfunction
function y = plip(x)
y = x^2 - 5*x + 2;
endfunction
function y = fcninput(fcnname, x)
// fcnname is the name of the function to be evaluated at x
y = feval(x,fcnname)
endfunction
fcninput(plip,1:10)
ans =
- 2. - 4. - 4. - 2. 2. 8. 16. 26. 38. 52.
fcninput(plop,1:10)
ans =
1. 3. 5. 7. 9. 11. 13. 15. 17. 19.
|
