Spline SNOOPY - Rolando Quispe Mamani

CODIGO MATLAB 2007
x0=1:0.1:2
x1=2:0.1:5
x2=5:0.1:6
x3=6:0.1:7
x4=7:0.1:8
x5=8:0.1:10
x6=10:0.1:13
x7=13:0.1:17
s0=(3)+((0.7857)*(x0-1))+((0)*(x0-1).^2)+((-0.0857)*(x0-1).^3)
s1=(3.7)+((0.5285667)*(x1-2))+((-0.2571)*(x1-2).^2)+((0.034378)*(x1-2).^3)
s2=(3.9)+((-0.08583)*(x2-5))+((0.0523)*(x2-5).^2)+((0.3335333)*(x2-5).^3)
s3=(4.2)+((1.0194)*(x3-6))+((1.0529)*(x3-6).^2)+((-0.5723)*(x3-6).^3)
s4=(5.7)+((1.4083)*(x4-7))+((-0.664)*(x4-7).^2)+((0.1557)*(x4-7).^3)
s5=(6.6)+((0.547333)*(x5-8))+((-0.1969)*(x5-8).^2)+((0.0241167)*(x5-8).^3)
s6=(7.1)+((0.049166)*(x6-10))+((-0.0522)*(x6-10).^2)+((-0.0028778)*(x6-10).^3)
s7=(6.7)+((-0.34173)*(x7-13))+((-0.0781)*(x7-13).^2)+((0.00650833)*(x7-13).^3)

x8=17:0.1:20
x9=20:0.1:23
x10=23:0.1:24
x11=24:0.1:25
x12=25:0.1:27
x13=27:0.1:27.7
s8=(4.5)+((1.105733)*(x8-17))+((0)*(x8-17).^2)+((-0.0302667)*(x8-17).^3)
s9=(7)+((0.2885)*(x9-20))+((-0.2724)*(x9-20).^2)+((0.0254111111111111)*(x9-20).^3)
s10=(6.1)+((-0.6598667)*(x10-23))+((-0.0437)*(x10-23).^2)+((0.203567)*(x10-23).^3)
s11=(5.6)+((-0.136533)*(x11-24))+((0.567)*(x11-24).^2)+((-0.230467)*(x11-24).^3)
s12=(5.8)+((0.306)*(x12-25))+((-0.1244)*(x12-25).^2)+((-0.0893)*(x12-25).^3)
s13=(5.2)+((-1.26333523809524)*(x13-27))+((-0.6602)*(x13-27).^2)+((0.314380952380953)*(x13-27).^3)

x14=27.7:0.1:28
x15=28:0.1:29
x16=29:0.1:30
s14=(4.1)+((0.748576)*(x14-27.7))+((0)*(x14-27.7).^2)+((-0.910111)*(x14-27.7).^3)
s15=(4.3)+((0.5028)*(x15-28))+((-0.8191)*(x15-28).^2)+((0.1163)*(x15-28).^3)
s16=(4.1)+((-0.78653)*(x16-29))+((-0.4702)*(x16-29).^2)+((0.15673)*(x16-29).^3)

x17=8:0.1:9
x18=9:0.1:11
x19=11:0.1:19
x20=19:0.1:26
x21=26:0.1:27.2
x22=27.2:0.1:27.7
s17=(1.8)+((-0.383433)*(x17-8))+((0)*(x17-8).^2)+((0.083433)*(x17-8).^3)
s18=(1.5)+((-0.1332)*(x18-9))+((0.2503)*(x18-9).^2)+((-0.05435)*(x18-9).^3)
s19=(1.8)+((0.216133)*(x19-11))+((-0.0758)*(x19-11).^2)+((0.00531667)*(x19-11).^3)
s20=(1.4)+((0.0243952380952381)*(x20-19))+((0.0518)*(x20-19).^2)+((-0.00848095238095238)*(x20-19).^3)
s21=(1.2)+((-0.497173)*(x21-26))+((-0.1263)*(x21-26).^2)+((0.9134722)*(x21-26).^3)
s22=(2)+((3.145933)*(x22-27.2))+((3.1622)*(x22-27.2).^2)+((-2.10813)*(x22-27.2).^3)

x23=1:0.1:5
x24=5:0.1:7
x25=7:0.1:8.5
s23=(3)+((0.0020666)*(x23-1))+((0)*(x23-1).^2)+((-0.00169167)*(x23-1).^3)
s24=(2.9)+((-0.0792)*(x24-5))+((-0.0203)*(x24-5).^2)+((0.01745)*(x24-5).^3)
s25=(2.8)+((0.048933)*(x25-7))+((0.0844)*(x25-7).^2)+((-0.0187556)*(x25-7).^3)

x26=6.8:0.1:8
x27=8:0.1:9
x28=9:0.1:11
x29=11:0.1:12
x30=12:0.1:12.5
x31=12.5:0.1:12.7
s26=(4.5)+((0.907093)*(x26-6.8))+((0)*(x26-6.8).^2)+((-0.05122)*(x26-6.8).^3)
s27=(5.5)+((0.6858)*(x27-8))+((-0.1844)*(x27-8).^2)+((-0.001399)*(x27-8).^3)
s28=(6)+((0.31273)*(x28-9))+((-0.1886)*(x28-9).^2)+((0.01611667)*(x28-9).^3)
s29=(6)+((-0.24833)*(x29-11))+((-0.0919)*(x29-11).^2)+((0.140233)*(x29-11).^3)
s30=(5.8)+((-0.01144)*(x30-12))+((0.3288)*(x30-12).^2)+((-3.0118)*(x30-12).^3)
s31=(5.5)+((-1.94148)*(x31-12.5))+((-4.1889)*(x31-12.5).^2)+((6.981500)*(x31-12.5).^3)

x32=11:0.1:12
x33=12:0.1:12.5
x34=12.5:0.1:12.7
s32=(3)+((0.973433)*(x32-11))+((0)*(x32-11).^2)+((-0.273433)*(x32-11).^3)
s33=(3.7)+((0.1531833)*(x33-12))+((-0.8203)*(x33-12).^2)+((5.0278667)*(x33-12).^3)
s34=(4.2)+((3.103800)*(x34-12.5))+((6.7215)*(x34-12.5).^2)+((-11.2025)*(x34-12.5).^3)

x35=11:0.1:11.5
x36=11.5:0.1:12
x37=12:0.1:13
x38=13:0.1:13.8
x39=13.8:0.1:14.2
x40=14.2:0.1:14.5
s35=(3)+((-0.5580833)*(x35-11))+((0)*(x35-11).^2)+((0.63233)*(x35-11).^3)
s36=(2.8)+((-0.0838166)*(x36-11.5))+((0.9485)*(x36-11.5).^2)+((-0.761733)*(x36-11.5).^3)
s37=(2.9)+((0.2933667)*(x37-12))+((-0.1941)*(x37-12).^2)+((0.0007333)*(x37-12).^3)
s38=(3)+((-0.0926533)*(x38-13))+((-0.1919)*(x38-13).^2)+((0.18933)*(x38-13).^3)
s39=(2.9)+((-0.0361866)*(x39-13.8))+((0.2625)*(x39-13.8).^2)+((-3.555083335)*(x39-13.8).^3)
s40=(2.7)+((-1.5326133)*(x40-14.2))+((-4.0036)*(x40-14.2).^2)+((4.44844444)*(x40-14.2).^3)

x41=15.5:0.1:16
x42=16:0.1:17
x43=17:0.1:19
x44=19:0.1:22
x45=22:0.1:24.5
x46=24.5:0.1:25
s41=(1.5)+((2.0978167)*(x41-15.5))+((0)*(x41-15.5).^2)+((-1.191267)*(x41-15.5).^3)
s42=(2.4)+((1.20433)*(x42-16))+((-1.7869)*(x42-16).^2)+((0.682567)*(x42-16).^3)
s43=(2.5)+((-0.321867)*(x43-17))+((0.2608)*(x43-17).^2)+((-0.062433)*(x43-17).^3)
s44=(2.4)+((-0.0279666666666666)*(x44-19))+((-0.1138)*(x44-19).^2)+((0.03363)*(x44-19).^3)
s45=(2.2)+((0.1974167)*(x45-22))+((0.1889)*(x45-22).^2)+((-0.0943467)*(x45-22).^3)
s46=(2.4)+((-0.6271)*(x46-24.5))+((-0.5187)*(x46-24.5).^2)+((0.3458)*(x46-24.5).^3)


plot(x0,s0,'k','LineWidth',2);grid;title('Grafico');hold on;
plot(x0,s0,'k','LineWidth',2);
plot(x1,s1,'k','LineWidth',2);
plot(x2,s2,'k','LineWidth',2);
plot(x3,s3,'k','LineWidth',2);
plot(x4,s4,'k','LineWidth',2);
plot(x5,s5,'k','LineWidth',2);
plot(x6,s6,'k','LineWidth',2);
plot(x7,s7,'k','LineWidth',2);

plot(x8,s8,'r','LineWidth',2);
plot(x9,s9,'r','LineWidth',2);
plot(x10,s10,'r','LineWidth',2);
plot(x11,s11,'r','LineWidth',2);
plot(x12,s12,'r','LineWidth',2);
plot(x13,s13,'r','LineWidth',2);

plot(x14,s14,'g','LineWidth',2);
plot(x15,s15,'g','LineWidth',2);
plot(x16,s16,'g','LineWidth',2);

plot(x17,s17,'k','LineWidth',2);
plot(x18,s18,'k','LineWidth',2);
plot(x19,s19,'k','LineWidth',2);
plot(x20,s20,'k','LineWidth',2);
plot(x21,s21,'k','LineWidth',2);
plot(x22,s22,'k','LineWidth',2);

plot(x23,s23,'b','LineWidth',2);
plot(x24,s24,'b','LineWidth',2);
plot(x25,s25,'b','LineWidth',2);

plot(x26,s26,'k','LineWidth',2);
plot(x27,s27,'k','LineWidth',2);
plot(x28,s28,'k','LineWidth',2);
plot(x29,s29,'k','LineWidth',2);
plot(x30,s30,'k','LineWidth',2);
plot(x31,s31,'k','LineWidth',2);

plot(x32,s32,'r','LineWidth',2);
plot(x33,s33,'r','LineWidth',2);
plot(x34,s34,'r','LineWidth',2);

plot(x35,s35,'k','LineWidth',2);
plot(x36,s36,'k','LineWidth',2);
plot(x37,s37,'k','LineWidth',2);
plot(x38,s38,'k','LineWidth',2);
plot(x39,s39,'k','LineWidth',2);
plot(x40,s40,'k','LineWidth',2);

plot(x41,s41,'b','LineWidth',2);
plot(x42,s42,'b','LineWidth',2);
plot(x43,s43,'b','LineWidth',2);
plot(x44,s44,'b','LineWidth',2);
plot(x45,s45,'b','LineWidth',2);
plot(x46,s46,'b','LineWidth',2);

hold off;