Đầu tiên xin nói lí do là mình tìm 1 tiếng đồng hồ mà không biết post bài này ở đâu (có thể do mình tìm không kĩ) nên xin mod và các bạn thông cảm, nếu được thì chuyển bài giúp mình.Mình đang học matlad. Thầy cho đề tài viết như sau:
Dùng chương trình matlad tìm nghiệm gần đúng của phương trình (X^3)-X-1=0 bằng phương pháp dây cung (với sai số của x và sai số của f(x) tự nhập vào từ chương tình nhé).
Thuật toán thì mình có nhưng do thầy chẳng chỉ dẫn gì mà bắt dùng matlad làm nên mình chẳng biết phải viết câu lệnh như thế nào cả.
Nếu các bạn không có thời (gian) gian thì có thể tìm giúp cho mình mấy cái ebook hay bất kì chỉ dẫn cơ bản nào về câu lênh cho file.m hay scrip cho matlad để mình có thể tự tìm hiểu vậy. Thank tất cả người rất rất nhiều, mong nhận được sự giúp đỡ sớm nhất từ tất cả người vì 4/12/2011 mình phải nộp bài rồi. Khổ thật, thầy cho vội quá
---------- Bài thêm lúc 13:46 ---------- Bài trước là lúc 13:45 ----------
Đây là thuật toán C++ sách cho tất cả người nhé
int daycung (double (*f) (double), double a, double b,)
{
clrscr();
double &x, double &errx, double &erry, int &buoclap;
if (f(a)==0) {x=a; errx=erry=0; buoclap=0; return true;}
if (f(b)==0) {x=a; errx=erry=0; buoclap=0; return true;}
if (f(a)*f(b)>0)
{ cout<
int k=1; double c,cp,aa,bb;aa=a;bb=b;
kvector xa,xb;
xa[0]=a;xb[0]=b;
c= (a*f(b)-f(a)*b))/(f(b)-f(a));
if (f(a)*f(c)<0) b=c; else a=c;
xa[1]=a;xb[1]=b;
while (true)
{cp=c;
c= (a*f(b)-f(a)*b))/(f(b)-f(a));
if (f(a)*f(c)<0) b=c; else a=c;
k++;
xa[k]=a;xb[k]=b;
if (fabs(c-cp)
if (k>kmax)
{ cout<
delay (1000); return false;
}
}
x=c;
errx=fabs(c-cp);erry=fabs(f(c));buoclap=k;
}
Dùng chương trình matlad tìm nghiệm gần đúng của phương trình (X^3)-X-1=0 bằng phương pháp dây cung (với sai số của x và sai số của f(x) tự nhập vào từ chương tình nhé).
Thuật toán thì mình có nhưng do thầy chẳng chỉ dẫn gì mà bắt dùng matlad làm nên mình chẳng biết phải viết câu lệnh như thế nào cả.
Nếu các bạn không có thời (gian) gian thì có thể tìm giúp cho mình mấy cái ebook hay bất kì chỉ dẫn cơ bản nào về câu lênh cho file.m hay scrip cho matlad để mình có thể tự tìm hiểu vậy. Thank tất cả người rất rất nhiều, mong nhận được sự giúp đỡ sớm nhất từ tất cả người vì 4/12/2011 mình phải nộp bài rồi. Khổ thật, thầy cho vội quá
---------- Bài thêm lúc 13:46 ---------- Bài trước là lúc 13:45 ----------
Đây là thuật toán C++ sách cho tất cả người nhé
int daycung (double (*f) (double), double a, double b,)
{
clrscr();
double &x, double &errx, double &erry, int &buoclap;
if (f(a)==0) {x=a; errx=erry=0; buoclap=0; return true;}
if (f(b)==0) {x=a; errx=erry=0; buoclap=0; return true;}
if (f(a)*f(b)>0)
{ cout<
int k=1; double c,cp,aa,bb;aa=a;bb=b;
kvector xa,xb;
xa[0]=a;xb[0]=b;
c= (a*f(b)-f(a)*b))/(f(b)-f(a));
if (f(a)*f(c)<0) b=c; else a=c;
xa[1]=a;xb[1]=b;
while (true)
{cp=c;
c= (a*f(b)-f(a)*b))/(f(b)-f(a));
if (f(a)*f(c)<0) b=c; else a=c;
k++;
xa[k]=a;xb[k]=b;
if (fabs(c-cp)
if (k>kmax)
{ cout<
delay (1000); return false;
}
}
x=c;
errx=fabs(c-cp);erry=fabs(f(c));buoclap=k;
}