C

>> C=[0 1;5 5;10 1]
C=
0 1
5 5
10 1
>> [A, b]=maprox(C)
A=
1 0 0
1 5 25
1 10 100
b=
>> [ a ] = met_gass(A, b)
1 0 0
0 5 25
0 0 50
1
4
-8
1.0000
1.6000
-0.1600
a=
1.0000
1.6000
-0.1600
f=
@(x)a(1,1)+a(2,1)*x+a(3,1)*x.^2
>> C=[0 -1;5 4;10 -1]
C=
0 -1
5 4
10 -1
>> [A, b]=maprox(C)
A=
1 0 0
1 5 25
1 10 100
b=
-1
-1
1 0 0
0 5 25
0 0 50
-1
-10
-1.0000
2.0000
-0.2000
a=
-1.0000
2.0000
-0.2000
>> g=@(x) a(1,1) + a(2,1)*x + a(3,1)*x.^2
g=
@(x)a(1,1)+a(2,1)*x+a(3,1)*x.^2
>> pc(10,1,20,f,g)
>> C=[0 -1;5 4;10 -1]
C=
0 -1
5 4
10 -1
>> [A, b]=maprox(C)
A=
1 0 0
1 5 25
1 10 100
b=
-1
4
-1
1 0 0
0 5 25
0 0 50
-1
-10
-1.0000
2.0000
-0.2000
a=
-1.0000
2.0000
-0.2000
>> j=@(x) a(1,1) + a(2,1)*x + a(3,1)*x.^2
j=
@(x)a(1,1)+a(2,1)*x+a(3,1)*x.^2
>> [cc,cn]=pl(10,10,j)
>> [cc, cn]=puente1(10,.5,30)

10
-2
0 2 4 6 8 10 12 14 16 18
C=[0 -2;5 5;10 -2]
C=
0 -2
5 5
10 -2
>> [A, b]=maprox(C)
A=
1 0 0
1 5 25
1 10 100
b=
-2
-2
1 0 0
0 5 25
0 0 50
-2
-14
-2.0000
2.8000
-0.2800
a=
-2.0000
2.8000
-0.2800
>> f=@(x) a(1,1) + a(2,1)*x + a(3,1)*x.^2
f=
@(x)a(1,1)+a(2,1)*x+a(3,1)*x.^2
>> C=[1 -2;5 4;9 -2]

C=
1 -2
5 4
9 -2
>> [A, b]=maprox(C)
A=
1 1 1
1 5 25
1 9 81
b=
-2
-2
1 1 1
0 4 24
0 0 32
-2
-12
-5.3750
3.7500
-0.3750
a=
-5.3750
3.7500
-0.3750
>> g=@(x) a(1,1) + a(2,1)*x + a(3,1)*x.^2
g=
@(x)a(1,1)+a(2,1)*x+a(3,1)*x.^2
>> [C,cn]=pc(10,1,40,f,g)
10
-2
0 2 4 6 8 10 12 14 16 18
-2
-4
-6
-4 -2 0 2 4 6 8 10 12 14
grid on
hold on
x1=[-5 31];
y1=[0 0];
plot(x1,y1)
x=[0 4 8 .5 4 7.5 1 4 7 1.5 4 6.5 2 4 6 2.5 4 5.5 3 4 5];
y=[0 6 0 0 5.5 0 0 5 0 0 4.5 0 0 4 0 0 3.5 0 0 3 0];
plot(x,y)
a=[9 13 17 9.5 13 16.5 10 13 16 10.5 13 15.5 11 13 15 11.5 13 14.5 12 13

14];
b=[0 6 0 0 5.5 0 0 5 0 0 4.5 0 0 4 0 0 3.5 0 0 3 0];
plot(a,b)
c=[ 18 22 26 18.5 22 25.5 19 22 25 19.5 22 24.5 20 22 24 20.5 22 23.5 21

22 23];
d=[0 6 0 0 5.5 0 0 5 0 0 4.5 0 0 4 0 0 3.5 0 0 3 0];
plot(c,d)

C

Uploaded by

Copyright:

Available Formats

C

Uploaded by

Document Information

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

C

Uploaded by

Copyright:

Available Formats

>> C=[0 1;5 5;10 1]

>> [A, b]=maprox(C)

>> C=[0 -1;5 4;10 -1]

>> [A, b]=maprox(C)

>> C=[0 -1;5 4;10 -1]

>> [A, b]=maprox(C)

>> j=@(x) a(1,1) + a(2,1)x + a(3,1)x.^2

>> [cc, cn]=puente1(10,.5,30)

C=[0 -2;5 5;10 -2]

>> [A, b]=maprox(C)

>> f=@(x) a(1,1) + a(2,1)x + a(3,1)x.^2

>> C=[1 -2;5 4;9 -2]

>> [A, b]=maprox(C)

>> g=@(x) a(1,1) + a(2,1)x + a(3,1)x.^2

a=[9 13 17 9.5 13 16.5 10 13 16 10.5 13 15.5 11 13 15 11.5 13 14.5 12 13

c=[ 18 22 26 18.5 22 25.5 19 22 25 19.5 22 24.5 20 22 24 20.5 22 23.5 21

You might also like