سورس ها و آموزش هاي Maple

[ravegoat]

سلام

من می خوام 2 معادله را باهم حل کنم (البته با شرایط مرزی)
برنامه زیر را نوشته ام ;

eq1 := diff(f(x), x, x, x)+(1/2)*f(x)*(diff(f(x), x, x)) = 0;
eq2 := diff(g(x), x, x)+.4*f(x)*(diff(g(x), x)) = 0;
eq3 := f(0) = 0, (D(f))(0) = 0, (D(f))(10) = 1;
;eq4 := g(0) = 1, g(10) = 0
dd := dsolve({eq1, eq2, eq3, eq4}, {f(x), g(x)}, numeric, output = listprocedure)
g := eval(g(x), dd);
f := eval(f(x), dd)
for k from 5.3 by 0.5e-1 to 5.4 do
print(dd, eval(g(x), x = k))
end do

حالا می خواهم مشتق تابع g را در x=0 بدست بیاورم اما نمی دونم چجوری

لطفا راهنمایی کنید

با تشکر

با سلام!
دوست گرامی عضویت شما در شهر سخت افزار تبریک می گم.برای مشاهده این لینک/عکس می بایست عضو شوید ! برای عضویت اینجا کلیک کنید

اگر شما ضابطه ی تابع g رو داشته باشید، با دستور diff می تونید مشتق اون رو به دست بیارید و با دستور eval مقدار مشتق رو در نقطه ی دلخواه محاسبه کنید:

dg := diff(g, x);
eval(dg, x = 0)

همچنین می تونید با داشتن مقدار g در همسایگی نقطه ی صفر مشتق رو با روش Five-point stencil تخمین بزنید.

موفق باشید
آرمین