��#�#�################>###��
#################G###########�#######����####O###H#######J###F###���������������
��������������������������������������������������������������������������������
��������������������������������������������������������������������������������
��������������������������������������������������������������������������������
��������������������������������������������������������������������������������
��������������������������������������������������������������������������������
���#7
###�#�###############<�####bjbjU#U################### ###4`##7|##7|
##�_##############################��##########��##########��##################l#
####�#######�###�#######�#######�#######�#######�###############�#######�f######
�f######�f##8###�f##�###�g##,###�#######S�##z###�i##�###�m##"###�m######�m######
�m######�n######�n######�n######?|######A|######A|######A|##7###x|
##�###H ##�####�##$### ##
#
�
###�##�###<�##�###################�#######�n######################�n######�n####
##�n######�n######<�###
�######�n######�n######�n######�n######�#######�m######�#######�m######?|
##############�n######################################################�n######?|
######�n##�###�n######cs##:####t##,###�#######�#################################
##############################ct######�m######�i##
###��`�#��#�###�b##�f######�n######Ot##
###########ct##�####�##0###S�######Yt##
###�######�n######�######ct######�n######�#######�#######�#######�#######�#####
##�#########�###C#H#�#�#N#G# #5#:# #N#�#I# #S#U#Y# #V#�# #X#�#P# #X#�# #H#�#M#
#
#�#1#.# #N#�#I# #S#U#Y# #L#A#G#R#A#N#G#E#
#
#T#r#o#n#g# #t#h#�#c# #t#�# #n#h#i#�#u# #k#h#i# #p#h#�#i# #p#h#�#c#
#h#�#i# #m#�#t# #h#�#m# #y# #=# #f#(#x#)# #t#�#i# #m#�#i# #g#i#�# #t#r#�# #x#
#t#r#o#n#g# #m#�#t# ###o#�#n# #[#a#,# #b#]# #n#�#o# ###�# #m#�# #c#h#�#
#b#i#�#t# #m#�#t# #s#�# #n#h#�#t# ###�#n#h# #c#�#c# #g#i#�# #t#r#�# #c#�#a#
#h#�#m# #t#�#i# #m#�#t# #s#�# ###i#�#m# #c#h#o# #t#r#�#�#c#.# #C#�#c# #g#i#�#
#t#r#�# #n#�#y# ###�#�#c# #c#u#n#g# #c#�#p# #q#u#a# #t#h#�#c# #n#g#h#i#�#m#
#h#a#y# #t#�#n#h# #t#o#�#n#.# #V#�# #v#�#y# #n#�#y# #s#i#n#h# #v#�#n# ###�#
#t#o#�#n# #h#�#c# #l#�# #t#r#�#n# ###o#�#n# #a# #(# #x# #(# #b# #c#h#o# #m#�#t#
#l#o#�#t# #c#�#c# ###i#�#m# #x#i# #(# #i#=# #0#,# #1#,# #2#.#.#.#)# #v#�#
#t#�#i# #c#�#c# ###i#�#m# #x#i# #n#�#y# #g#i#�# #t#r#�# #c#�#a# #h#�#m# #l#�#
#y#i# #=# #f#(#x#i#)# ###�# #b#i#�#t#.# #B#�#y# #g#i#�# #t#a# #c#�#n# #t#�#m#
###a# #t#h#�#c# #:#
#
#P#n#(#x#)# #=# #a#o#x#n# #+# #a#1#x#n#-#1# # #+# #(#+#a#n#-#1#x# # #+#
#a#n# #
#s#a#o# #c#h#o# #P#n#(#x#i#)# #=# #f#(#x#i#)# #=# #y#i#.# ###a# #t#h#�#c#
#P#n#(#x#)# ###�#�#c# #g#�#i# #l#�# ###a# #t#h#�#c# #n#�#i# #s#u#y# #c#�#a#
#h#�#m# #y#=#f#(#x#)#.# #T#a# #c#h#�#n# ###a# #t#h#�#c# ###�# #n#�#i# #s#u#y#
#h#�#m# #y# #=# #f#(#x#)# #v#�# ###a# #t#h#�#c# #l#�# #l#o#�#i# #h#�#m# ###�#n#
#g#i#�#n#,# #l#u#�#n# #c#�# ###�#o# #h#�#m# #v#�# #n#g#u#y#�#n# #h#�#m#.#
#V#i#�#c# #t#�#n#h# #g#i#�# #t#r#�# #c#�#a# #n#�# #t#h#e#o# #t#h#u#�#t#
#t#o#�#n# #H#o#r#n#e#r# #c#i#n#g# ###�#n# #g#i#�#n#.#
#
#B#�#y# #g#i#�# #t#a# #x#�#y# #d#�#n#g# ###a# #t#h#�#c# #n#�#i# #s#u#y#
#k#i#�#u# #L#a#g#r#a#n#g#e#.# #G#�#i# #L#i# #l#�# ###a# #t#h#�#c#:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#2# # #######
#
#R#�# #r#�#n#g# #l#�# #L#i#(#x#)# #l#�# #m#�#t# ###a# #t#h#�#c# #b#�#c#
#n# #v#�# #:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#T#a# #g#�#i# ###a# #t#h#�#c# #n#�#y# #l#�# ###a# #t#h#�#c# #L#a#g#r#a#n#g#e#
#c#�# #b#�#n#.#
#B#�#y# #g#i#�# #t#a# #x#�#t# #b#i#�#u# #t#h#�#c# #:#
# ### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#T#a# #t#h#�#y# #P#n#(#x#)# #l#�# #m#�#t# ###a# #t#h#�#c# #b#�#c# #n# #v#�#
#c#�#c# #L#i#(#x#)# #l#�# #c#�#c# ###a# #t#h#�#c# #b#�#c# #n# #v#�# #t#h#o#�#
#m#�#n# ###i#�#u# #k#i#�#n# #P#n#(#x#i#)# #=# #f#(#x#i#)# #=# #y#i#.# #T#a#
#g#�#i# #n#�# #l#�# ###a# #t#h#�#c# #n#�#i# #s#u#y# #L#a#g#r#a#n#g#e#.#
#V#�#i# #n# #=# #1# #t#a# #c#�# #b#�#n#g# #
#x###x#0###x#1#####y###y#0###y#1#######a# #t#h#�#c# #n#�#i# #s#u#y# #s#�# #l#�#
#:#
#
#P#1#(#x#)# #=# #y#o#L#0#(#x#)# #+# #y#1#L#1#(#x#1#)#
#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # ####### # # # # # # ###
#E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#n#�#n#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#N#h#�# #v#�#y# #P#1#(#x#)# #l#�# #m#�#t# ###a# #t#h#�#c# #b#�#c# #n#h#�#t#
###�#i# #v#�#i# #x#
#V#�#i# #n# #=# #2# #t#a# #c#�# #b#�#n#g# #
#x###x#0###x#1###x#2#####y###y#0###y#1###y#2#######a# #t#h#�#c# #n#�#i# #s#u#y#
#s#�# #l#�# #:#
#
#P#2#(#x#)# #=# #y#o#L#0#(#x#)# #+# #y#1#L#1#(#x#1#)# #+#
#y#2#L#2#(#x#2#)#
#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#
#N#h#�# #v#�#y# #P#1#(#x#)# #l#�# #m#�#t# ###a# #t#h#�#c# #b#�#c# #h#a#i#
###�#i# #v#�#i# #x#
#
#T#r#�#n# #c#�# #s#�# #t#h#u#�#t# #t#o#�#n# #t#r#�#n# #t#a# #c#�#
#c#h#�#�#n#g# #t#r#�#n#h# #t#�#m# ###a# #t#h#�#c# #n#�#i# #s#u#y# #c#�#a#
#m#�#t# #h#�#m# #k#h#i# #c#h#o# #t#r#�#�#c# #c#�#c# ###i#�#m# #v#�# #s#a#u#
###�# #t#�#n#h# #t#r#�# #s#�# #c#�#a# #n#�# #t#�#i# #m#�#t# #g#i#�# #t#r#�#
#n#�#o# ###�# #n#h#�# #s#a#u# #:#
#
#C#h#�#�#n#g# #t#r#�#n#h# #5#-#1#
#
###i#n#c#l#u#d#e# #<#c#o#n#i#o#.#h#>#
###i#n#c#l#u#d#e# #<#s#t#d#i#o#.#h#>#
###i#n#c#l#u#d#e# #<#c#t#y#p#e#.#h#>#
###d#e#f#i#n#e# #m#a#x# # # #2#1#
#
#i#n#t# #m#a#x#k#q#,#n#;#
#f#l#o#a#t#
#x#[#m#a#x#]#,#y#[#m#a#x#]#,#a#[#m#a#x#]#,#x#x#[#m#a#x#]#,#y#y#[#m#a#x#]#;#
#f#l#o#a#t# #x#0#,#p#0#;#
#
#v#o#i#d# #m#a#i#n#(#)#
# # #{#
#
#i#n#t# #i#,#k#;#
#
#c#h#a#r# #o#k# #;#
#
#
#v#o#i#d# #v#a#o#s#o#l#i#e#u#(#v#o#i#d#)#;#
#
#f#l#o#a#t# #l#a#g#r#a#n#g#e#(int,float [],float [],float);
void inkq(void);
}
clrscr();
printf("%24cNOI SUY DA THUC LAGRANGE\n",' ');
vaosolieu();
k=0;
ok='c';
while (ok=='c')
{
printf("Tinh gia tri cua y voi x la x0 = ");
scanf("%f",&x0);
p0=lagrange(n,x,y,x0);
printf("Gia tri cua y = %15.5f\n",p0);
printf("\n");
k=k+1;
maxkq=k;
xx[k]=x0;
yy[k]=p0;
flushall();
printf("Tinh tiep khong(c/k)?");
scanf("%c",&ok);
}
inkq();
void vaosolieu()
{
int i,t;
char ok;
printf("\n");
printf("Ham y = f(x)\n");
printf("So cap (x,y) nhieu nhat la max = 20\n");
printf("So diem da cho truoc n = ");
scanf("%d",&n);
for (i=1;i<=n;i++)
{
printf("x[%d] = ",i);
scanf("%f",&x[i]);
printf("y[%d] = ",i);
scanf("%f",&y[i]);
}
printf("\n");
printf("
SO LIEU BAN VUA NHAP\n");
printf("
x
y\n");
for (i=1;i<=n;i++)
printf("%8.4f
%8.4f\n",x[i],y[i]);
ok=' ';
t=1;
flushall();
while (t)
{
printf("\nCo sua so lieu khong(c/k):?");
scanf("%c",&ok);
if (toupper(ok)=='C')
{
printf("Chi so cua phan tu can sua i = ");
scanf("%d",&i);
printf("Gia tri moi : ");
printf("x[%d] = ",i);
scanf("%f",&x[i]);
printf("y[%d] = ",i);
scanf("%f",&y[i]);
flushall();
}
if (toupper(ok)!='C')
t=0;
}
}
float lagrange(int n,float x[max],float y[max],float x0)
{
int i,k;
float g0;
}
p0=0.0;
for (k=1;k<=n;k++)
{
g0=1.0;
for (i=1;i<=n;i++)
if (i!=k)
g0=g0*(x0-x[i])/(x[k]-x[i]);
p0=p0+y[k]*g0;
}
return(p0);
void inkq()
{
int i,j,k;
printf("\n");
printf("%24cBANG SO LIEU\n",' ');
printf("%18cx %24cy\n",' ',' ');
for (i=1;i<=n;i++)
printf("%20.4f %25.4f\n",x[i],y[i]);
printf("\n");
printf("%24cKET QUA TINH TOAN\n",' ');
printf("%14cx %10cy\n",' ',' ');
for (k=1;k<=maxkq;k++)
printf("%15.5f %15.5f\n",xx[k],yy[k]);
getch();
}
################################################################################
################################################################################
##########G#i#�# #s#�# #t#a# #c#�# #b#�#n#g# #c#�#c# #g#i#�# #t#r#�# #x#,# #y#
#:#
#x###0###3###-#2###2###4#####y###0###-#3#.#7#5###1#0###-#2###4#####v#�#y#
#t#h#e#o# #c#h#�#�#n#g# #t#r#�#n#h# #t#�#i# # #x# #=# #2#.#5# #y# #=# ##3#.#3#5#4#9#.#
#
#�#2#.# #N#�#I# #S#U#Y# #N#E#W#T#O#N#
#
#B#�#y# #g#i#�# #t#a# #x#�#t# #m#�#t# #c#�#c#h# #k#h#�#c# ###�# #x#�#y#
#d#�#n#g# ###a# #t#h#�#c# #n#�#i# #s#u#y# #g#�#i# #l#�# #p#h#�#�#n#g# #p#h#�#p#
#N#e#w#t#o#n#.# #T#r#�#�#c# #h#�#t# #t#a# ###�#a# #v#�#o# #m#�#t# #k#h#�#i#
#n#i#�#m# #m#�#i# #l#�# #t#�# #h#i#�#u# #
#
#G#i#�# #s#�# #h#�#m# #y# #=# #y#(#x#)# #c#�# #g#i#�# #t#r#�# #c#h#o#
#t#r#o#n#g# #b#�#n#g# #s#a#u# #:#
#####x###x#0###x#1###x#2###(###x#n##1###x#n###########y###y#0###y#1###y#2###(###y#n#-#1###y#n#######
#T#�#
#h#i#�#u# #c#�#p# #1# #c#�#a# #y# #t#�#i# #x#i#,# #x#j# #l#�# #:#
#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#
#T#�# #h#i#�#u# #c#�#p# #h#a#i# #c#�#a# #y# #t#�#i# #x#i#,# #x#j#,# #x#k#
#l#�# #:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#v#.#v#.#
#
#
#V#�#i# #y#(#x#)# #=# #P#n#(#x#)# #l#�# #m#�#t# ###a# #t#h#�#c#
#b#�#c# #n# #t#h#�# #t#�# #h#i#�#u# #c#�#p# #1# #t#�#i# #x#,# #x#0# #:#
#
#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#l#�# #m#�#t# ###a# #t#h#�#c# #b#�#c# #(#n# #-# #1#)#.# #T#�# #h#i#�#u# #c#�#p#
#2# #t#�#i# #x#,# #x#0#,# #x#1# #:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#l#�# #m#�#t# ###a# #t#h#�#c# #b#�#c# #(#n#-#2#)# #v#.#v# #v#�# #t#�#i# #t#�#
#h#i#�#u# #c#�#p# #(#n# #+# #1#)# #t#h#�# #:#
#
#
#P#n#[# #x#,# #x#o#,#.#.#,# #x#n#]# #=# # #0#
#T#�# #c#�#c# ###�#n#h# #n#g#h#)#a# #t#�# #h#i#�#u# #t#a# #s#u#y# #r#a# #:#
#
#P#n#(#x#)# #=# #P#n#(#x#0#)# #+# #(# #x#-# #x#0#)#P#n#[#x#,# #x#o#]#
#
#P#n#[#x#,# #x#0#]# #=# #P#n#[#x#0#,# #x#1#]# #+# #(# #x# #-#
#x#1#)#P#n#[#x#,# #x#o#,#x#1#]#
#
#P#n#[#x#,# #x#o#,# #x#1#]# #=# #P#n#[#x#0#,# #x#1#,# #x#2#]# #+# #(# #x#
#-# #x#2#)#P#n#[#x#,# #x#o#,# #x#1#,# #x#2#]#
#
#.#.#.#.#.#.#.#.#.#.#.#.#
#
#P#n#[#x#,# #x#o#,#.#.#,# #x#n#-#1#]# #=# #P#n#[#x#0#,# #x#1#,#.#.#,#
#x#n#]# #+# #(# #x# #-# #x#n#)#P#n#[#x#,# #x#o#,#.#.#,# #x#n#]#
#D#o# #
#P#n#[# #x#,# #x#o#,#.#.#,# #x#n#]# #=# # #0# #n#�#n# #t#�# ###�#
#t#a# #c#�# #:#
#P#n#(#x#)# #=# #P#n#(#x#0#)# #+# #(#x# #-# #x#0#)#P#n#[#x#o#,# #x#1#]# #+#
#(#x# #-# #x#0#)#(#x# #-# #x#1#)#P#n#[#x#0#,# #x#1#,# #x#2#]# #+#(#
#+#(#x# #-# #x#0#)#(#(#x# #-# #x#n#-#1#)#P#n#[#x#0#,#(#,# #x#n#]#
#N#�#u# #P#n#(#x#)# #l#�# ###a# #t#h#�#c# #n#�#i# #s#u#y# #c#�#a# #h#�#m# #y#
#=# #f#(#x#)# #t#h#�#:#
#
#P#n#(#x#i#)# #=# #f#(#x#i#)# #=# #y#i# #v#�#i# #i# #=# #0# #(# #n# #
#
#D#o# ###�# #c#�#c# #t#�# #h#i#�#u# #t#�# #c#�#p# #1# ###�#n# #c#�#p# #n#
#c#�#a# #P#n# #v#�# #c#�#a# #y# #l#�# #t#r#�#n#g# #n#h#a#u# #v#�# #n#h#�#
#v#�#y# #t#a# #c#�# #:#
#
# #P#n#(#x#)# #=# #y#0# #+# #(#x# #-# #x#0#)#y#[#x#0#,# #x#1#]# #+# #(#x#
#-# #x#0#)#(#x# #-# #x#1#)#y#[#x#0#,# #x#1#,# #x#2#]# #+# #.#.# #+#
# # # # # # # # # # # # # # # # # # # # # # # #(#x# #-# #x#0#)#(#x# #-#
#x#1#)#.#.#.#(#x# #-# #x#n#-#1#)#y#[#x#0#,#.#.#,#x#n#]#
###a# #t#h#�#c# #n#�#y# #g#�#i# #l#�# ###a# #t#h#�#c# #n#�#i# #s#u#y#
#N#e#w#t#o#n# #t#i#�#n# #x#u#�#t# #p#h#�#t# #t#�# #n#�#t# #x#0# #c#�#a# #h#�#m#
#y# #=# #f#(#x#)#.# #N#g#o#�#i# ###a# #t#h#�#c# #t#i#�#n# #c#�#n# #c#�# ###a#
#t#h#�#c# #n#�#i# #s#u#y# #N#e#w#t#o#n# #l#�#i# #x#u#�#t# #p#h#�#t# #t#�#
###i#�#m# #x#n# #c#�# #d#�#n#g# #n#h#�# #s#a#u# #:#
#
#P#n#(#x#)# #=# #y#n# #+# #(#x# #-# #x#n#)#y#[#x#n#,# #x#n#-#1#]# #+#
#(#x# #-# #x#n#)#(#x# #-# #x#n#-#1#)#y#[#x#n#,# #x#n#-#1#,#x#n#-#2#]# #+#.#.#+#
# # # # # # # # # # # # # # # # # # # # # #(#x# #-# #x#n#)#(#x# #-# #x#n##1#)#.#.#.#(#x# #-# #x#1#)#y#[#x#n#,#.#.#,# #x#0#]#
#T#r#�#�#n#g# #h#�#p# #c#�#c# #n#�#t# #c#�#c#h# ###�#u# #t#h#�# #x#i# #=# #x#0#
#+# #i#h# #v#�#i# #i# #=# #0#,# #1#,#.#.#,# #n#.# #T#a# #g#�#i# #s#a#i#
#p#h#�#n# #t#i#�#n# #c#�#p# #1# #t#�#i# #i# #l#�# #:#
#
#
#(#y#i# #=# #y#i#+#1# #-# #y#i#
#v#�# #s#a#i# #p#h#�#n# #t#i#�#n# #c#�#p# #h#a#i# #t#�#i# #i#:#
#
#
#(#2#y#i# #=# #(#(#(#y#i#)# #=# #y#i#+#2# #-# #2#y#i#+#1# #+# #y#i#
#
#
#.#.#.#.#.#.#.#.#.#
#v#�# #s#a#i# #p#h#�#n# #t#i#�#n# #c#�#p# #n# #l#�# #:#
#
#
#(#n#y#i# #=# #(#(#(#n#-#1#y#i#)#
#K#h#i# ###�# #t#a# #c#�#:#
#
#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # ####### #
#
#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # ####### #
#
#.#.#.#.#.#.#.#.#.#.#.#
#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#B#�#y# #g#i#�# ###�#t# #x# #=# #x#0# #+# #h#t# # #t#r#o#n#g# ###a# #t#h#�#c#
#N#e#w#t#o#n# #t#i#�#n# #t#a# ###�#�#c#:#
#
#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # ####### #
#t#h#�# #t#a# #n#h#�#n# ###�#�#c# ###a# #t#h#�#c# #N#e#w#t#o#n# #t#i#�#n#
#x#u#�#t# #p#h#�#t# #t#�# #x#0# #t#r#o#n#g# #t#r#�#�#n#g# #h#�#p# #n#�#t#
#c#�#c#h# ###�#u#.# #V#�#i# #n# #=# #1# #t#a# #c#�# #:#
#
#P#1#(#x#0# #+# #h#t#)# #=# #y#0# #+# #(#y#0#
#V#�#i# #n# #=# #2# #t#a# #c#�#:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#M#�#t# #c#�#c#h# #t#�#�#n#g# #t#�# #t#a# #c#�# #k#h#�#i# #n#i#�#m# #c#�#c#
#s#a#i# #p#h#�#n# #l#�#i# #t#�#i# #i#:#
#
#
#(#y#i# #=# #y#i# #-# #y#i#-#1#
#
#
#(#2#y#i# #=# #(#(#(#y#i#)# #=# #y#i# #-# #2#y#i#-#1# #+# #y#i#-#2#
#
#
#.#.#.#.#.#.#.#.#.#
#
#
#(#n#y#i# #=# #(#(#(#n#-#1#y#i#)#
#v#�# ###a# #t#h#�#c# #n#�#i# #s#u#y# #N#e#w#t#o#n# #l#�#i# #k#h#i# #c#�#c#
###i#�#m# #n#�#i# #s#u#y# #c#�#c#h# ###�#u#:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#
#V#�# #d#�#:# #C#h#o# #h#�#m# #n#h#�# #b#�#n#g# #s#a#u#:#
#
#x###0#.#1###0#.#2###0#.#3###0#.#4#####y###0#.#0#9#9#8#3###0#.#1#9#8#6#7###0#.#2
#9#5#5#2###0#.#3#8#9#4#2#####
#T#a# #t#�#n#h# #g#i#�# #t#r#�# #c#�#a# #h#�#m# #t#�#i# #0#.#1#4# #b#�#n#g#
###a# #t#h#�#c# #n#�#i# #s#u#y# #N#e#w#t#o#n# #v#�# #c#�#c# #m#�#c# #c#�#c#h#
###�#u# #h# #=# #0#.#1#.# #T#a# #c#�# #b#�#n#g# #s#a#i# #p#h#�#n# #s#a#u# #:#
#
#i###x###y###(#y###(#2#y###(#3#y#####0###0#.#1###0#.#0#9#9#8#3#################0
#.#0#9#8#8#4#########1###0#.#2###0#.#1#9#8#6#7######0#.#0#0#1#9#9#############0#.#0#9#6#8#5######0#.#0#0#0#9#6#####2###0#.#3###0#.#2#9#5#5#2######0#.#0#0#2#9#5#############0#.#0#9#3#9#0#########3###0#.#4###0#.#3#8#9#4#2######
#####
#T#a# #d#�#n#g# #c#�#n#g# #t#h#�#c# #N#e#w#t#o#n# #t#i#�#n# #v#�#i# ###i#�#m#
#g#�#c# #l#�# #x#0# #=# #0#.#1#;# #h# #=# #0#.#1#.# #V#�#i# #x# #=# #0#.#1#4#
#t#a# #c#�# #0#.#1#4# #=# #0#.#1# #+# #0#.#1#t# #n#�#n# #t# #=# #0#.#4# #v#�#
#k#�#t# #q#u#�# #l#�# #:# ### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#2# #
#######C#h#�#�#n#g# #t#r#�#n#h# #n#�#i# #s#u#y# #N#e#w#t#o#n# #n#h#�# #s#a#u#
#:#
#
#C#h#�#�#n#g# #t#r#�#n#h# #5#-#2#
#
#/#/#N#o#i# #s#u#y# #N#e#w#t#o#n#
###i#n#c#l#u#d#e# #<#c#o#n#i#o#.#h#>#
###i#n#c#l#u#d#e# #<#s#t#d#i#o#.#h#>#
###i#n#c#l#u#d#e# #<#c#t#y#p#e#.#h#>#
###d#e#f#i#n#e# #m#a#x# # # #1#1#
#
#v#o#i#d# #m#a#i#n#(#)#
# # {
int i,j,k,n,t;
float a[max],b[max],x[max],y[max];
char ok;
float x0,p;
clrscr();
printf("So diem da cho n = ");
scanf("%d",&n);
for (i=1;i<=n;i++)
{
printf("x[%d] = ",i);
scanf("%f",&x[i]);
printf("y[%d] = ",i);
scanf("%f",&y[i]);
}
printf("%10cBANG SO LIEU\n",' ');
printf("%8cx%30cy\n",' ',' ');
for (i=1;i<=n;i++)
printf("%4c%8.4f%23c%8.4f\n",' ',x[i],' ',y[i]);
ok=' ';
t=0;
flushall();
while (t)
{
printf("Co sua so lieu khong(c/k): ");
scanf("%c",&ok);
if (toupper(ok)=='C')
{
printf("Chi so cua phan tu can sua i = ");
scanf("%d",&i);
printf("Gia tri moi : ");
printf("x[%d] = ",i);
scanf("%f",&x[i]);
printf("y[%d] = ",i);
scanf("%f",&y[i]);
flushall();
}
if (toupper(ok)!='C')
t=0;
}
a[1]=y[1];
for (j=1;j<=n-1;j++)
{
for (i=1;i<=n-j;i++)
y[i]=(y[i+1]-y[i])/(x[i+j]-x[i]);
a[j+1]=y[1];
}
b[n]=a[n];
for (k=n-1;k>=1;k--)
{
for (j=n-1;j>=1;j--)
b[j]=a[j] ;
for (i=n-1;i>=k;i--)
a[i]=a[i]-b[i+1]*x[k];
}
for (i=n;i>=1;i--)
printf("He so bac %d la :%8.4f\n",i-1,a[i]);
printf("\n");
k=0;
ok='c';
flushall();
while (ok=='c')
{
printf("Tinh gia tri cua y tai x = ");
scanf("%f",&x0);
p=0;
for (k=n;k>=1;k--)
p=p*x0+a[k];
printf("Tri so noi suy tai x0 = %4.2f la : %10.5f\n",x0,p);
getch();
printf("Ban co muon tinh tiep cac diem khac khong(c/k)");
do
scanf("%c",&ok);
while ((ok!='c')&&(ok!='k'));
}
}
D�ng ch#########################################################�#�#n#g#
#t#r#�#n#h# #n#�#y# #n#�#i# #s#u#y# #c#�#c# #g#i#�# #t#r#�# #c#h#o# #t#r#o#n#g#
#b#�#n#g# #s#a#u# #
#
#0###0#.#2###0#.#4###0#.#6###0#.#8###1#.#0#####1###1#.#2#2#1#4#0#2#7#6###1#.#4#9
#1#8#2#4#7###1#.#8#2#2#1#1#8#8###2#.#2#2#5#5#4#0#9#3###2#.#7#1#8#2#8#1#8#3#####t
#a# #c#�# #c#�#c# #h#�# #s#�# #c#�#a# ###a# #t#h#�#c# #n#�#i# #s#u#y#:#
#0#.#0#1#3#9#(#b#�#c# #5#)#,# #0#.#0#3#4#9#(#b#�#c# #4#)#,#
#0#.#1#7#0#4#(#b#�#c#3#)#,# #0#.#4#9#9#1#(#b#�#c# #2#)#,# #1#.#0#0#0#1#(#b#�#c#
#1#)# #v#�# #1#.#0#0#0#0#(#b#�#c# #0#)#.#
#
#�#3#.# #N#�#I# #S#U#Y# #A#I#T#K#E#N#
#M#�#t# #d#�#n#g# #k#h#�#c# #c#�#a# ###a# #t#h#�#c# #n#�#i# #s#u#y# ###�#�#c#
#x#�#c# ###�#n#h# #b#�#n#g# #t#h#u#�#t# #t#o#�#n# #A#i#t#k#e#n#.# #G#i#�# #s#�#
#t#a# #c#�# #n# ###i#�#m# ###�# #c#h#o# #c#�#a# #h#�#m# #f#(#x#)#.# #N#h#�#
#v#�#y# #q#u#a# #h#a#i# ###i#�#m# #x#0# #v#�# #x#1# #t#a# #c#�# ###a# #t#h#�#c#
#n#�#i# #s#u#y# #L#a#g#r#a#n#g#e# #c#�#a# #h#�#m# #f#(#x#)# ###�#�#c# #v#i#�#t#
#d#�#�#i# #d#�#n#g#:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#l#�# #m#�#t# ###a# #t#h#�#c# #b#�#c# #1# #:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#K#h#i# #x# #=# #x#0# #t#h#�#:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#K#h#i# #x# #=# #x#1# #t#h#�#:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
###a# #t#h#�#c# #n#�#i# #s#u#y# #L#a#g#r#a#n#g#e# #c#�#a# #f#(#x#)# #q#u#a# #3#
###i#�#m# #x#0#,# #x#1#,# #x#2# #c#�# #d#�#n#g#:#
#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#v#�# #l#�# #m#�#t# ###a# #t#h#�#c# #b#�#c# #2#:#
# # # # # # #
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#K#h#i# #x# #=# #x#0# #t#h#�#:#
#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#K#h#i# #x# #=# #x#1# #t#h#�#:#
#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#K#h#i# #x# #=# #x#2# #t#h#�#:#
#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#T#�#n#g# #q#u#�#t# ###a# #t#h#�#c# #n#�#i# #s#u#y# #L#a#g#r#a#n#g#e# #q#u#a#
#n# ###i#�#m# #l#�#:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#N#h#�# #v#�#y# #t#a# #c#�# #t#h#�# #d#�#n#g# #p#h#�#p# #l#�#p# ###�# #x#�#c#
###�#n#h# #l#�#n# #l#�#�#t# #c#�#c# ###a# #t#h#�#c# #L#a#g#r#a#n#g#e#.# #S#�#
###�# #t#�#n#h# #t#o#�#n# #n#h#�# #v#�#y# #g#�#i# #l#�# #s#�# ###�#
#N#e#v#i#l#l#e#-#A#i#t#k#e#n#.#
#V#�# #d#�#:# #C#h#o# #c#�#c# #c#�#p# ###i#�#m# #(#0#,# #0#.#4#)#,# #(#1#.#4#,#
#1#.#5#)#,# #(#2#.#6#,# #1#.#8#)#,# #(#3#.#9#,# #2#.#6#)#,#t#�#n#h# #y# #t#�#i#
#x# #=# #2#
#
# ### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # ####### ### #E#M#B#E#D#
#E#q#u#a#t#i#o#n#.#3# # #######
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#C#h#�#�#n#g# #t#r#�#n#h# ###�#�#c# #v#i#�#t# #n#h#�# #s#a#u# #
#
#C#h#�#�#n#g# #t#r#�#n#h# #5#-#3#
#
#/#/#N#o#i# #s#u#y# #A#i#t#k#e#n#
###i#n#c#l#u#d#e# #<#c#o#n#i#o#.#h#>#
###i#n#c#l#u#d#e# #<#s#t#d#i#o#.#h#>#
###i#n#c#l#u#d#e# #<#c#t#y#p#e#.#h#>#
###d#e#f#i#n#e# #m#a#x# # # #1#1#
#
#v#o#i#d# #m#a#i#n#(#)#
# # #{#
#
#f#l#o#a#t# #x#[#m#a#x#]#,#y#[#m#a#x#]#,#y#d#[#m#a#x#]#;#
#
#float x1;
int j,k,n,n1;
#
clrscr();
printf("Cho so diem da co n = ");
scanf("%d",&n1);
n=n1-1 ;
for (k=0;k<=n;k++)
{
printf("x[%d] = ",k+1);
scanf("%f",&x[k]);
printf("y[%d] = ",k+1);
scanf("%f",&y[k]);
}
printf("Cho diem can tinh gia tri cua ham x1 = ");
scanf("%f",&x1);
for (k=0;k<=n-1;k++)
{
yd[k]=(y[k]*(x1-x[k+1])-y[k+1]*(x1-x[k]))/(x[k]-x[k+1]);
if (k!=0)
for (j=k-1;j>=0;j--)
yd[j]=(yd[j]*(x1-x[k+1])-yd[j+1]*(x1-x[j]))/(x[j]-x[k+1]);
}
printf("Gia tri ham tai x = %#6#.#3#f# #l#a# #y# #=# #
%#8#.#4#f#\#n#"#,#x#1#,#y#d#[#0#]#)#;#
#
#g#e#t#c#h#(#)#;#
# # #}#
#
#D#�#n#g# #c#h#�#�#n#g# #t#r#�#n#h# #n#�#y# ###�# #n#�#i# #s#u#y# #c#�#c#
#c#�#p# #s#�# #(#1#,# #3#)#,# #(#2#,# #5#)#,# #(#3#,# #7#)#,# #(#4#,# #9#)#
#v#�# #(#5#,# #1#1#)# #t#�#i# #x# #=# #2#.#5# #t#a# #c#�# #y# #=# #6#.#
#
#�#4#.# #X#�#P# #X#�# #H#�#M# #B#�#N#G# #P#H#�#�#N#G# #P#H#�#P# #B#�#N#H#
#P#H#�#�#N#G# #B#�# #N#H#�#T#
#
#T#r#o#n#g# #c#�#c# #m#�#c# #t#r#�#�#c# #t#a# ###�# #n#�#i# #s#u#y#
#g#i#�# #t#r#�# #c#�#a# #h#�#m#.# #B#�#i# #t#o#�#n# ###�# #l#�# #c#h#o# #m#�#t#
#h#�#m# #d#�#�#i# #d#�#n#g# #b#�#n#g# #s#�# #v#�# #p#h#�#i# #t#�#m# #g#i#�#
#t#r#�# #c#�#a# #h#�#m# #t#�#i# #m#�#t# #g#i#�# #t#r#�# #c#�#a# ###�#i# #s#�#
#k#h#�#n#g# #n#�#m# #t#r#o#n#g# #b#�#n#g#.#
#
#T#r#o#n#g# #t#h#�#c# #t#�#,# #b#�#n# #c#�#n#h# #b#�#i# #t#o#�#n# #n#�#i#
#s#u#y# #t#a# #c#�#n# #g#�#p# #m#�#t# #d#�#n#g# #b#�#i# #t#o#�#n# #k#h#�#c#.#
###�# #l#�# #t#�#m# #c#�#n#g# #t#h#�#c# #t#h#�#c# #n#g#h#i#�#m# #c#�#a# #m#�#t#
#h#�#m#.# #N#�#i# #d#u#n#g# #b#�#i# #t#o#�#n# #l#�# #t#�# #m#�#t# #l#o#�#t#
#c#�#c# ###i#�#m# #c#h#o# #t#r#�#�#c# #(#c#�# #t#h#�# #l#�# #c#�#c# #g#i#�#
#t#r#�# #c#�#a# #m#�#t# #p#h#�#p# ###o# #n#�#o# ###�#)# #t#a# #p#h#�#i# #t#�#m#
#m#�#t# #h#�#m# #x#�#p# #x#�# #c#�#c# #g#i#�# #t#r#�# ###�# #c#h#o#.# #T#a#
#s#�# #d#�#n#g# #p#h#�#�#n#g# #p#h#�#p# #b#�#n#h# #p#h#�#�#n#g# #t#�#i#
#t#h#i#�#u# ###�# #g#i#�#i# #b#�#i# #t#o#�#n#.# #G#i#�# #s#�# #c#�# #m#�#u#
#q#u#a#n# #s#�#t# #(#x#i#,#y#i# #)# #c#�#a# #h#�#m# #y#=#f#(#x#)#.# #T#a#
#c#h#�#n# #h#�#m# #f#(#x#)# #c#�# #d#�#n#g# #:#
#
#f#(#x#)# #=# #a#0#f#0#(#x#)# #+# #a#1#f#1#(#x#)# #+#
#a#2#f#2#(#x#)#.#.#.#
#
#
#
#
#
#(#1#)#
#T#r#o#n#g# ###�# #c#�#c# #h#�#m# #f#0#(#x#)#,# #f#1#(#x#)#,# #f#2#(#x#)#
#v#.#v#.# #l#�# #(#m#+#1#)# #h#�#m# ###�#c# #l#�#p# #t#u#y#�#n# #t#�#n#h# #m#�#
#t#a# #c#�# #t#h#�# #c#h#�#n# #t#u#�# #�# #v#�# #c#�#c# #h#�# #s#�# #a#i# #l#�#
#t#h#a#m# #s#�# #c#h#�#a# #b#i#�#t# #m#�# #t#a# #p#h#�#i# #x#�#c# ###�#n#h#
#d#�#a# #v#�#o# #h#�# #h#�#m# ###�# #c#h#�#n# #v#�# #c#�#c# ###i#�#m# #q#u#a#n#
#s#�#t#.# #S#a#i# #s#�# #g#i#�#a# #t#r#�# ###o# ###�#�#c# #v#�# #t#r#�#
#t#�#n#h# #t#h#e#o# #(#1#)# #l#�# #:#
#
#e#i# #=# #y#i# #-# #f#(#x#i#)#
#
#
#
#
#
#
#
#
#
#(#2#)#
#S#a#i# #s#�# #n#�#y# #c#�# #t#h#�# #�#m# #h#a#y# #d#�#�#n#g# #t#u#�# #t#�#n#g#
#g#i#�# #t#r#�# #c#�#a# #y#i#.# #K#h#i# #d#�#n#g# #p#h#�#�#n#g# #p#h#�#p#
#b#�#n#h# #p#h#�#�#n#g# #b#�# #n#h#�#t# #t#a# #x#�#t# #b#�#n#h# #p#h#�#�#n#g#
#c#�#a# #s#a#i# #s#�# #t#�#i# #m#�#t# ###i#�#m#:#
#
#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # ####### #
#
#
#
#
#
# #
#(#3#)#
#V#�#i# #n# ###i#�#m# #t#�#n#g# #b#�#n#h# #p#h#�#�#n#g# #c#�#a# #s#a#i# #s#�#
#s#�# #l#�# #:#
#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#R#�# #r#�#n#g# #S# #l#�# #h#�#m# #c#�#a# #c#�#c# #g#i#�# #t#r#�# #c#�#n#
#t#�#m# #a#i#.#v#�# #c#h#�#n#g# #t#a# #s#�# #c#h#�#n# #c#�#c# #a#i# #s#a#o#
#c#h#o# #S# ###�#t# #g#i#�# #t#r#�# #m#i#n#,# #n#g#h#)#a# #l#�# #c#�#c# ###�#o#
#h#�#m# ### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######p#h#�#i# #b#�#n#g#
#k#h#�#n#g#.#T#a# #s#�# #x#�#t# #c#�#c# #t#r#�#�#n#g# #h#�#p# #c#�# #t#h#�#.#
#1#.# #H#�#m# #x#�#p# #x#�# #c#�# #d#�#n#g# ###a# #t#h#�#c#:# #T#r#o#n#g#
#t#r#�#�#n#g# #h#�#p# #t#�#n#g# #q#u#�#t# #t#a# #c#h#�#n# #h#�# #h#�#m# #x#�#p#
#x#�# #l#�# #m#�#t# ###a# #t#h#�#c#,# #n#g#h#)#a# #l#�#:#
#
#f#(#x#)# #=# #a#0# #+# #a#1#x# #+# #a#2#x#2# #+#.#.#.#+# #a#m#x#m#
#V#�#y# #h#�#m# #S# #l#�# #:#
#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # ####### #
#T#h#e#o# ###i#�#u# #k#i#�#n# ###�#o# #h#�#m# ### #E#M#B#E#D#
#E#q#u#a#t#i#o#n#.#3# # #######t#a# #n#h#�#n# ###�#�#c# #h#�# #p#h#�#�#n#g#
#t#r#�#n#h#:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#
###�#y# #l#�# #m#�#t# #h#�# #p#h#�#�#n#g# #t#r#�#n#h# #t#u#y#�#n# #t#�#n#h#.#
#G#i#�#i# #n#�# #t#a# #n#h#�#n# ###�#�#c# #c#�#c# #g#�#a# #t#r#�# #a#i#.#
#S#a#u# ###�#y# #l#�# #c#h#�#�#n#g# #t#r#�#n#h# #v#i#�#t# #t#h#e#o# #t#h#u#�#t#
#t#o#�#n# #t#r#�#n#.#
#
#
#C#h#�#�#n#g# #t#r#�#n#h# #5#-#4#
#
#/#/#X#a#p# #x#i# #d#a# #t#h#u#c#
###i#n#c#l#u#d#e# #<#c#o#n#i#o#.#h#>#
###i#n#c#l#u#d#e# #<#s#t#d#i#o#.#h#>#
###include
#define max
11
void main()
{
int i,j,k,m,n,p,kp,t;
float a[max],x[max],y[max],y1[max];
float b[max][max];
char ok;
float s,sx,s1,c,d;
clrscr();
printf("PHUONG PHAP BINH PHUONG TOI THIEU");
printf("\n");
printf("Cho bac cua da thuc xap xi m = ");
scanf("%d",&m);
printf("So diem da cho n = ");
scanf("%d",&n);
for (i=1;i<=n;i++)
{
printf("x[%d] = ",i);
scanf("%f",&x[i]);
printf("y[%d] = ",i);
scanf("%f",&y[i]);
}
x[0]=1;
printf("\n");
printf("%4cBANG SO LIEU\n",' ');
printf("%8cx%30cy\n",' ',' ');
for (i=1;i<=n;i++)
printf("%4c%8.4f%20c%8.4f\n",' ',x[i],' ',y[i]);
ok=' ';
t=1;
flushall();
while (t)
{
printf("Co sua so lieu khong(c/k): ");
scanf("%c",&ok);
if (toupper(ok)=='C')
{
printf("Chi so cua phan tu can sua i = ");
scanf("%d",&i);
printf("Gia tri moi : ");
printf("x[%d] = ",i);
scanf("%f",&x[i]);
printf("y[%d] = ",i);
scanf("%f",&y[i]);
flushall();
}
if (toupper(ok)!='C')
t=0;
}
//for (i=0;i<=n;i++)
//a[i]=0.0;
printf("\n");
printf("CAC GIA TRI DA CHO");
printf("\n");
printf("X = ");
for (i=1;i<=n;i++)
printf("%c%8.3f",' ',x[i]);
printf("\n");
printf("Y = ");
for (i=1;i<=n;i++)
printf("%c%8.3f",' ',y[i]);
printf("\n");
for (p=0;p<=m;p++)
{
y1[p]=0.0;
for (i=1;i<=n;i++)
{
sx=1.0;
for (j=1;j<=p;j++)
sx*=x[i];
y1[p]+=y[i]*sx;
}
}
for (p=0;p<=m;p++)
for (k=0;k<=m;k++)
{
kp=k+p;
b[p][k]=0.0;
for (i=1;i<=n;i++)
{
sx=1.0;
for (j=1;j<=kp;j++)
sx*=x[i];
b[p][k]+=sx;
}
}
}
for (i=0;i<=m-1;i++)
{
c=1.0/b[i][i];
for (k=i+1;k<=m;k++)
{
d=b[i][k];
for (j=i+1;j<=m;j++)
b[k][j]-=b[i][j]*c*d;
y1[k]-=y1[i]*c*d;
b[i][k]*=c;
}
y1[i]*=c;
}
y1[m]/=b[m][m];
for (i=m-1;i>=0;i--)
for (j=i+1;j<=m;j++)
y1[i]-=b[i][j]*y1[j];
printf("\n");
printf("CAC HE SO CUA DA THUC CAN TIM");
printf("\n");
for (i=0;i<=m;i++)
printf("a[%d] = %10.5f\n",i,y1[i]);
getch();
#####################V#�#i# #c#�#c# #g#i#�# #t#r#�# #x#,#y# ###o# ###�#�#c#
#t#h#e#o# #b#�#n#g#
#
#
#x###7###8###9###1#0###1#1###1#2###1#3#####y###7#,#4###8#,#4###9#,#1###9#,#4###9
#,#5###9#,#5###9#,#4#####
#t#a# #c#�# #n# #=# #7# #v#�# #c#h#�#n# #m# #=# #2# #v#�# #t#�#n#h# ###�#�#c#
#t#h#e#o# #c#h#�#�#n#g# #t#r#�#n#h# #c#�#c# #h#�# #s#�# #:# #
#a#0# #=# #-#0#.#1#1#1#9#0#5#;# #a#1# #=# #2#.#5#4#5#2#3#8#;# #a#2# #=# ##4#.#8#5#7#1#4#3#
#v#�# #h#�#m# #x#�#p# #x#�# #s#�# #l#�#:# #f#(#x#)# #=# #-#0#.#1#1#1#9#0#5# #+#
#2#.#5#4#5#2#3#8#x# #-# #4#.#8#5#7#1#4#3#x#2#
#
#2#.#H#�#m# #d#�#n#g# #A#e#c#x#:# #K#h#i# #c#�#c# #s#�# #l#i#�#u# #t#h#�#
#h#i#�#n# #m#�#t# #s#�# #b#i#�#n# ###�#i# ###�#n# ###i#�#u# #t#a# #d#�#n#g#
#h#�#m# #x#�#p# #x#�# #l#�# #y# #=# #A#e#c#x#.# #L#�#y# #l#o#g#a#r#i#t# #h#a#i#
#v#�# #t#a# #c#�# #:#
#
#l#n#y# #=# #l#n#A# #+# #c#x#l#n#e#
#T#h#e#o# ###i#�#u# #k#i#�#n# ###�#o# #h#�#m# ### #E#M#B#E#D#
#E#q#u#a#t#i#o#n#.#3# # #######t#a# #c#�# #h#�# #p#h#�#�#n#g# #t#r#�#n#h# #:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#G#i#�#i# #h#�# #p#h#�#�#n#g# #t#r#�#n#h# #n#�#y# #t#a# #c#�# #c#�#c# #h#�#
#s#�# #A# #v#�# #c# #:#
#
#C#h#�#�#n#g# #t#r#�#n#h# #5#-#5#
#
#/#/#x#a#p#_#x#i#_#e#_#m#u#;#
###i#n#c#l#u#d#e# #<#c#o#n#i#o#.#h#>#
###i#n#c#l#u#d#e# #<#s#t#d#i#o#.#h#>#
###i#n#c#l#u#d#e# #<#c#t#y#p#e#.#h#>#
###i#n#c#l#u#d#e# #<#m#a#t#h#.#h#>#
###d#e#f#i#n#e# #m#a#x# # # #1#1#
#v#o#i#d# #m#a#i#n#(#)#
# # #{#
#
#
#
#
#i#n#t# #i#,#n#,#t#;#
#f#l#o#a#t# #x#[#m#a#x#]#,#y#[#m#a#x#]#;#
#c#h#a#r# #o#k#;#
#f#l#o#a#t# #a#,#b#,#c#,#d#,#e#,#f#,#d#1,d2,d3;
clrscr();
printf("PHUONG PHAP BINH PHUONG TOI THIEU");
printf("\n");
printf("So diem da cho n = ");
scanf("%d",&n);
for (i=1;i<=n;i++)
{
printf("x[%d] = ",i);
scanf("%f",&x[i]);
printf("y[%d] = ",i);
scanf("%f",&y[i]);
}
x[0]=1.0;
printf("%4cBANG SO LIEU\n",' ');
printf("%8cx%30cy\n",' ',' ');
for (i=1;i<=n;i++)
printf("%4c%8.4f%23c%8.4f\n",' ',x[i],' ',y[i]);
ok=' ';
t=1;
while (t)
{
printf("Co sua so lieu khong(c/k): ");
scanf("%c",&ok);
if (toupper(ok)=='C')
{
printf("Chi so cua phan tu can sua i = ");
scanf("%d",&i);
printf("Gia tri moi : ");
printf("x[%d] = ",i);
scanf("%f",&x[i]);
printf("y[%d] = ",i);
scanf("%f",&y[i]);
}
if (toupper(ok)!='C')
t=0;
}
printf("CAC GIA TRI DA CHO");
printf("\n");
printf("X = ");
for (i=1;i<=n;i++)
printf("%c%8.3f",' ',x[i]);
printf("\n");
printf("Y = ");
for (i=1;i<=n;i++)
printf("%c%8.3f",' ',y[i]);
printf("\n");
a=0.0;
for (i=1;i<=n;i++)
a+=x[i];
b=n;
c=0.0;
for (i=1;i<=n;i++)
c+=log(y[i]);
d=0.0;
for (i=1;i<=n;i++)
d+=x[i]*x[i];
e=0.0;
for (i=1;i<=n;i++)
e+=x[i]*log(y[i]);
}
d1=a*a-d*b;
d2=c*a-e*b;
d3=a*e-c*d;
c=d2/d1;
a=d3/d1;
printf("\n");
printf("He so A = %8.4f",exp(a));
printf(" va so mu c = %8.4",c);
printf("\n");
printf("\nBANG CAC GIA TRI TINH TOAN");
printf("\n");
printf("%5cx%28cy\n",' ',' ');
for (i=1;i<=n;i++)
{
printf("%8.4f%21c%8.4f\n",x[i],' ',exp(a)*exp(c*x[i]));
}
getch();
##################################V#�#i# #c#�#c# #g#i#�# #t#r#�# #x#,#y# ###o#
###�#�#c# #t#h#e#o# #b#�#n#g# #
#
#x###0###2###4###6###8###1#0###1#2#####y###1#2#8#0###6#3#5###3#2#4###1#6#2###7#6
###4#3###1#9#####
#t#a# #c#�# #n# #=# #7# #v#�# #t#�#n#h# ###�#�#c# #t#h#e#o# #c#h#�#�#n#g#
#t#r#�#n#h# #c#�#c# #h#�# #s#�#:# #A# #=# #1#2#8#5#.#4#4# #v#a# #c# #=# ##0#.#3#4#7#6# #v#�# #h#�#m# #x#�#p# #x#�# #s#�# #l#�#:# #f#(#x#)# #=#
#1#2#8#5#.#4#4#
#
#3#.# #H#�#m# #d#�#n#g# #A#x#q#:# #K#h#i# #c#�#c# #s#�# #l#i#�#u# #t#h#�#
#h#i#�#n# #m#�#t# #s#�# #b#i#�#n# ###�#i# ###�#n# ###i#�#u# #t#a# # #c#i#n#g#
#c#�# #t#h#�# #d#�#n#g# #h#�#m# #x#�#p# #x#�# #l#�# #y# #=# #A#x#q#.# #L#�#y#
#l#o#g#a#r#i#t# #h#a#i# #v#�# #t#a# #c#�#:#
#
#l#n#y# #=# #l#n#A# #+# #q#l#n#x#
#T#h#e#o# ###i#�#u# #k#i#�#n# ###�#o# #h#�#m# #t#r#i#�#t# #t#i#�#u# #t#a# #c#�#
#h#�# #p#h#�#�#n#g# #t#r#�#n#h# #:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#G#i#�#i# #h#�# #p#h#�#�#n#g# #t#r#�#n#h# #n#�#y# #t#a# #c#�# #c#�#c# #h#�#
#s#�# #A# #v#�# #q#:#
#
#C#h#�#�#n#g# #t#r#�#n#h# #5#-#6#
#
#//xap_xi_x_mu;
#include
#include
#include
#include
#define max
11
void main()
{
int i,n,t;
float x[max],y[max];
char ok;
float a,b,c,d,e,f,d1,d2,d3;
clrscr();
printf("PHUONG PHAP BINH PHUONG TOI THIEU");
printf("\n");
printf("So diem da cho n = ");
scanf("%d",&n);
for (i=1;i<=n;i++)
{
printf("x[%d] = ",i);
scanf("%f",&x[i]);
printf("y[%d] = ",i);
scanf("%f",&y[i]);
}
x[0]=1.0;
printf("%4cBANG SO LIEU\n",' ');
printf("%8cx%30cy\n",' ',' ');
for (i=1;i<=n;i++)
printf("%4c%8.4f%23c%8.4f\n",' ',x[i],' ',y[i]);
ok=' ';
flushall();
t=1;
while (t)
{
printf("Co sua so lieu khong(c/k): ");
scanf("%c",&ok);
if (toupper(ok)=='C')
{
printf("Chi so cua phan tu can sua i = ");
scanf("%d",&i);
printf("Gia tri moi : ");
printf("x[",i,"] = ");
scanf("%f",&x[i]);
printf("y[%d] = ",i);
scanf("%f",&y[i]);
}
if (toupper(ok)!='C')
t=0;
}
printf("\n");
printf("\nCAC GIA TRI DA CHO");
printf("\n");
printf("X = ");
for (i=1;i<=n;i++)
printf("%c%8.3f",' ',x[i]);
printf("\n");
printf("Y = ");
for (i=1;i<=n;i++)
printf("%c%8.3f",' ',y[i]);
printf("\n");
a=0.0;
for (i=1;i<=n;i++)
a+=log(x[i]);
b=n;
c=0.0;
for (i=1;i<=n;i++)
c+=log(y[i]);
d=0.0;
for (i=1;i<=n;i++)
d+=log(x[i])*log(x[i]);
e=0.;
for (i=1;i<=n;i++)
e+=log(x[i])*log(y[i]);
d1=a*a-d*b;
d2=c*a-e*b;
d3=a*e-c*d;
c=d2/d1;
a=d3/d1;
printf("\n");
printf("He so A = %8.4f",exp(a));
printf(" va so mu q = %8.4f\n",c);
print#f#(#"#\#n#"#)#;#
#
#p#r#i#n#t#f#(#"#\#n#B#A#N#G# #C#A#C# #G#I#A# #T#R#I# #T#I#N#H#
#T#O#A#N#\#n#"#)#;#
#
#p#r#i#n#t#f#(#"#%#5#c#x#%#2#7#c#y#\#n#"#,#'# #'#,#'# #'#)#;#
#
#f#o#r# #(#i#=#1#;#i#<#=#n#;#i#+#+#)#
#
# # #{#
#
#
#p#r#i#n#t#f#(#"#%#8#.#4#f#%#2#0#c#%#8#.#4#f#\#n#"#,#x#[#i#]#,#'#
#'#,#e#x#p#(#a#)#*#e#x#p#(#c#*#l#o#g#(#x#[#i#]#)#)#)#;#
#
# # #}#
#
#g#e#t#c#h#(#)#;#
# # #}#
#
#V#�#i# #c#�#c# #g#i#�# #t#r#�# #x#,#y# ###o# ###�#�#c# #t#h#e#o# #b#�#n#g# #
#
#x###1###2###4###5###6#####y###7#.#1###2#7#.#8###6#2#.#1###1#1#0###1#6#1#####
#t#a# #c#�# #n# #=# #5# #v#�# #t#�#n#h# ###�#�#c# #t#h#e#o# #c#h#�#�#n#g#
#t#r#�#n#h# #c#�#c# #h#�# #s#�#:# #A# #=# #7#.#1#6#4#1# #v#�# #q# #=#
#1#.#9#5#3#1# #v#�# #h#�#m# #x#�#p# #x#�# #s#�# #l#�#:# #f#(#x#)# #=#
#1#2#8#5#.#4#4#x#1#.#9#5#3#1#
#
#4#.# #H#�#m# #l#�#�#n#g# #g#i#�#c#:# #K#h#i# #q#u#a#n# #h#�# #y# #=# #f#(#x#)#
#c#�# #d#�#n#g# #t#u#�#n# #h#o#�#n# #t#a# #d#�#n#g# #h#�#m# #x#�#p# #x#�# #l#�#
#t#�# #h#�#p# #t#u#y#�#n# #t#�#n#h# #c#�#a# #c#�#c# #h#�#m# #s#i#n# #v#�#
#c#o#s#i#n# #d#�#n#g#:#
#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
###�# ###�#n# #g#i#�#n# #t#r#�#�#c# #h#�#t# #t#a# #x#�#t# #h#�#m# #c#h#�# #c#�#
#m#�#t# #s#�# #h#�#n#g# #s#i#n#-#c#o#s#,#n#g#h#)#a# #l#�# #:#
#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#H#�#m# #S# #s#�# #c#�# #d#�#n#g# #:#
#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#T#h#e#o# ###i#�#u# #k#i#�#n# ###�#o# #h#�#m# #t#r#i#�#t# #t#i#�#u# #t#a# #c#�#
#h#�# #p#h#�#�#n#g# #t#r#�#n#h# ###�#i# #v#�#i# #c#�#c# #h#�# #s#�# #d#�#n#g#:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#D#o# #:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#n#�#n# #h#�# #p#h#�#�#n#g# #t#r#�#n#h# #c#�# #d#�#n#g# ###�#n# #g#i#�#n# #:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#G#i#�#i# #h#�# #t#a# #c#�# #:#
#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#T#r#o#n#g# #t#r#�#�#n#g# #h#�#p# #t#�#n#g# #q#u#�#t#,#m#�#t# #c#�#c#h#
#t#�#�#n#g# #t#�# #t#a# #c#�#:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#C#h#�#�#n#g# #t#r#�#n#h# #t#�#m# #c#�#c# #h#�# #s#�# #a#i# #v#�# #b#i#
###�#�#c# #t#h#�# #h#i#�#n# #n#h#�# #s#a#u#:#
#
#C#h#�#�#n#g# #t#r#�#n#h# #5#-#7#
#
#/#/#x#a#p#_#x#i#_#s#i#n#_#c#o#s#;#
###i#n#c#l#u#d#e# #<#c#o#n#i#o#.#h#>#
###i#n#c#l#u#d#e# #<#s#t#d#i#o#.#h#>#
###i#n#c#l#u#d#e# #<#c#t#y#p#e#.#h#>#
###i#n#c#l#u#d#e# #<#m#a#t#h#.#h#>#
###d#e#f#i#n#e# #m#a#x# # # #1#1#
###d#e#f#i#n#e# #p#i# # #3#.#1#5#1#5#9#
#
#v#o#i#d# #m#a#i#n#(#)#
# # #{#
#
#i#n#t# #i#,#j#,#m#,#n#,#t#;#
#
#f#l#o#a#t# #x#[#m#a#x#]#,#y#[#m#a#x#]#,#a#[#m#a#x#]#,#b#[#m#a#x#]#;#
#
#char ok;
float omg,t1;
clrscr();
printf("PHUONG PHAP BINH PHUONG TOI THIEU");
printf("\n");
printf("Cho so so hang sin-cos m = ");
scanf("%d",&m);
printf("Cho chu ki T = ");
scanf("%f",&t1);
printf("So diem da cho n = ");
scanf("%d",&n);
for (i=1;i<=n;i++)
{
printf("x[%d] = ",i);
scanf("%f",&x[i]);
printf("y[%d] = ",i);
scanf("%f",&y[i]);
}
x[0]=1.0;
printf("%4cBANG SO LIEU\n",' ');
printf("%8cx%30cy\n",' ',' ');
for (i=1;i<=n;i++)
printf("%4c%8.4f%23c%8.4f\n",' ',x[i],' ',y[i]);
ok=' ';
t=1;
flushall();
while (t)
{
printf("Co sua so lieu khong(c/k): ");
scanf("%c",&ok);
if (toupper(ok)=='C')
{
printf("Chi so cua phan tu can sua i = ");
scanf("%d",&i);
printf("Gia tri moi : ");
printf("x[%d] = ",i);
scanf("%f",&x[i]);
printf("y[%d] = ",i);
scanf("%f",&y[i]);
flushall();
}
if (toupper(ok)!='C')
t=0;
}
printf("\nCAC GIA TRI DA CHO\n");
printf("\n");
printf("
X
Y\n");
for (i=1;i<=n;i++)
printf("%c%8.3f%c%8.3f\n",' ',x[i],' ',y[i]);
printf("\n");
a[0]=0.0;
omg=2*pi/t1;
for (i=1;i<=n;i++)
a[0]+=y[i];
a[0]/=n;
for (j=1;j<=m;j++)
{
a[j]=0.0;
for (i=1;i<=n;i++)
a[j]+=y[i]*cos(j*omg*x[i]);
a[j]=2*a[j]/n;
}
for (j=1;j<=m;j++)
{
b[j]=0.0;
for (i=1;i<=n;i++)
b[j]+=y[i]*sin(j*omg*x[i]);
b[j]=2*b[j]/n;
}
printf("\n");
printf("TAN SO GOC OMEGA = %10.5f\n",omg);
printf("HE SO HANG\n");
printf("a[0] = %8.4f\n",a[0]);
printf("CAC HE SO BAC CAO\n");
printf("%5ccos%25csin\n",'# #'#,#'# #'#)#;#
#
#f#o#r# #(#i#=#1#;#i#<#=#m#;#i#+#+#)#
#
# # #p#r#i#n#t#f#(#"#%#8#.#4#f#%#2#1#c#%#6#.#4#f#\#n#"#,#a#[#i#]#,#'#
#'#,#b#[#i#]#)#;#
#
#g#e#t#c#h#(#)#;#
# # #}#
#
#V#�#i# #h#�#m# #c#h#o# #b#�#n#g# #b#�#n#g# #s#�# #:#
#
#x###0###0#.#1#5###0#.#3###0#.#4#5###0#.#6###0#.#7#5###0#.#9###1#.#0#5###1#.#2##
#1#.#3#####y###2#.#2###1#.#5#9#5###1#.#0#3#1###0#.#7#2#2###0#.#7#8#6###1#.#2###1
#.#8#1###2#.#3#6#9###2#.#6#7#8###2#.#6#1#4#####C#h#�#n# #s#�# #h#�# #s#�#
#s#i#n#-#c#o#s# #m# #=# #1#,# #s#�# ###i#�#m# #c#h#o# #t#r#�#�#c# #n# #=#
#1#0#,# #c#h#u# #k#�# # #T# #=# #1#5# #t#a# #n#h#�#n# ###�#�#c# #k#�#t# #q#u#�#
#t#�#n#h# #a#0# #=# #1#.#7#;# #a#1# #=# #0#.#5#;# #b#1# #=# #-#0#.#8#6#6#1#
#v#�# #(# #=# #4#.#1#8#8#7#9#.# #N#h#�# #v#�#y# #h#�#m# #x#�#p# #x#�# #c#�#
#d#�#n#g#:#
#
#f#(#x#)# #=# #1#.#7# #+# #0#.#5#c#o#s#(#4#.#1#8#8#7#9#x#)# #-#
#0#.#8#6#6#1#s#i#n#(#4#.#1#8#8#7#9#x#)#
#
#5# #.#H#�#m# #h#�#u# #t#�#:# #K#h#i# #q#u#a#n# #h#�# #y# #=# #f#(#x#)# #c#�#
#d#�#n#g# ###�#�#n#g# #c#o#n#g# #b#�#o# #h#o#�# #h#a#y# #d#�#n#g#
#a#r#c#t#a#n#,# #t#a#n# #v#.#v# #t#a# #d#�#n#g# #h#�#m# #x#�#p# #x#�# #l#�#
#h#�#m# #h#�#u# #t#�# #d#�#n#g# ###�#n# #g#i#�#n#:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#L#�#y# #n#g#h#�#c#h# ###�#o# #c#�#a# #n#�# #t#a# #c#�# #:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
###�#t# #1#/#y# #=# #Y#,# #1#/#x# #=# #X#,# #b#/#a# #=# #B# #v#�# #1#/#a# #=#
#A# #p#h#�#�#n#g# #t#r#�#n#h# #t#r#�#n# #s#�# #c#�# #d#�#n#g#:#
#
#
#Y# #=# #A# #+# #B#X#
#v#�# #l#�# #m#�#t# ###a# #t#h#�#c# #b#�#c# #m#�#t#.#D#o# #v#�#y# #t#a# #c#�#
#h#�# #p#h#�#�#n#g# #t#r#�#n#h# ###�#i# #v#�#i# #c#�#c# #h#�# #s#�# #A# #v#�#
#B# #l#�#:#
### #E#M#B#E#D# #E#q#u#a#t#i#o#n#.#3# # #######
#v#�# #t#�# ###�# #t#�#n#h# ###�#�#c# #a# #v#�# #b#.# #C#h#�#�#n#g# #t#r#�#n#h#
#s#a#u# #m#�# #t#�# #t#h#u#�#t# #t#o#�#n# #t#r#�#n#
#
#C#h#�#�#n#g# #t#r#�#n#h# #5#-#8#
#
#/#/#x#a#p# #x#i# #h#u#u#_#t#y#;#
###i#n#c#l#u#d#e# #<#c#o#n#i#o#.#h#>#
###i#n#c#l#u#d#e# #<#s#t#d#i#o#.#h#>#
###i#n#c#l#u#d#e# #<#c#t#y#p#e#.#h#>#
###i#n#c#l#u#d#e# #<#m#ath.h>
#define k 11
void main()
{
float x[k],y[k];
float a,b,a1,b1,c,d,e;
int i,n,t;
char ok;
clrscr();
printf("PHUONG PHAP BINH PHUONG TOI THIEU");
printf("\n");
printf("So diem da cho n = ");
scanf("%d",&n);
for (i=1;i<=n;i++)
{
printf("x[%d] = ",i);
scanf("%f",&x[i]);
printf("y[%d] = ",i);
scanf("%f",&y[i]);
}
}
printf("%4cBANG SO LIEU\n",' ');
printf("%8cx%30cy\n",' ',' ');
for (i=1;i<=n;i++)
printf("%4c%8.4f%23c%8.4f\n",' ',x[i],' ',y[i]);
ok=' ';
t=1;
flushall();
while (t)
{
printf("Co sua so lieu khong(c/k): ");
scanf("%c",&ok);
if (toupper(ok)=='C')
{
printf("Chi so cua phan tu can sua i = ");
scanf("%d",&i);
printf("Gia tri moi : ");
printf("x[%d] = ",i);
scanf("%f",&x[i]);
printf("y[%d] = ",i);
scanf("%f",&y[i]);
flushall();
}
if (toupper(ok)!='C')
t=0;
}
printf("CAC GIA TRI DA CHO\n");
printf("\n");
printf("X = ");
for (i=1;i<=n;i++)
printf("%c%8.3f",' ',x[i]);
printf("\n");
printf("Y = ");
for (i=1;i<=n;i++)
printf("%c%8.3f",' ' ,y[i]);
printf("\n");
a=n;
b=0.0;
c=0.0;
d=0.0;
e=0.0;
for (i=1;i<=n;i++)
{
b+=1/x[i];
c+=1/y[i];
d+=1/(x[i]*x[i]);
e+=1/(x[i]*y[i]);
}
a1=(c*d-b*e)/(a*d-b*b);
b1=(a*e-b*c)/(a*d-b*b);
a=1/a1;
b=b1*a;
printf("\n");
printf("Cac he so cua ham huu ty\n");
printf("a = %10.5f b = %10.5f",a,b);
getch();
################################################################################
#######################################################################V#�#i#
#d#�#y# #s#�# #l#i#�#u# ###�# #c#h#o# #:#
#
#x###1###2###3###4###5#####y###0#.#3#3#3#3#3#3#3###0#.#5###0#.#6###0#.#6#6#6#6#6
###0#.#7#1#4#2#8#5#7#####
#t#a# #n#h#�#n# ###�#�#c# #t#�# #c#h#�#�#n#g# #t#r#�#n#h# #t#r#�# #s#�# #a# #=#
#1# #v#�# #b# #=# #2#
###P#A#G#E# # ###
#
#
###P#A#G#E# # ###1#3#0###
#
#
#
#
################################################################################
################################################################################
################################################################################
################################################################################
################################################################################
################################################################################
################################################################################
################################################################################
################################################################################
################################################################################
################################################################################
################################################################################
###############################################################
##�
##�
##(###X###�###�####
- Xem thêm -
.DOC 
452 
65 
140 
