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

[ravegoat]

این کار را کردم اما مقدار 0 را می دهد !!

آیا در برنامه ای نوشتم مشکلی هست؟

متاسفانه بنده میپل نصب ندارم که سورس شما رو بررسی کنم ولی در یه نگاه کلی که بهش انداختم مشکلی نداشت.

شما اول بررسی کنید که به ازای dd(0) در خروجی آیا مشتق g در x=0 ظاهر میشه.

در صورتی که این ترم وجود داشت سورس زیر رو امتحان کنید؛ شاید جواب داد:

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 = operator)
dg := eval(D(g)(x), dd);
dg(0)

در اولین فرصت هم اگر میپل پیدا کردم سورس ها رو مجددا بررسی می کنم.

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

[aseman1369_s]

متاسفانه بنده میپل نصب ندارم که سورس شما رو بررسی کنم ولی در یه نگاه کلی که بهش انداختم مشکلی نداشت.

شما اول بررسی کنید که به ازای dd(0) در خروجی آیا مشتق g در x=0 ظاهر میشه.

در صورتی که این ترم وجود داشت سورس زیر رو امتحان کنید؛ شاید جواب داد:

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 = operator)
dg := eval(D(g)(x), dd);
dg(0)

در اولین فرصت هم اگر میپل پیدا کردم سورس ها رو مجددا بررسی می کنم.

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

بالاخره تونستم خروجی بگیرم
اینم دستورش

dg := eval(diff(g(x), x), dd)
print(dd, eval(dg(x), x = 0))

با تشکر از پیگیری شما