মূলত আমি যা করতে চাইছি তা হল:

১. নিচের সমীকরণ থেকে L21ও LA1 এর সমাধান বের করতে হবে। আসল উদ্দেশ্য LA1এর মান বের করা।
২. তারপর সমাধানকে s এর সাপেক্ষে ডিফারনশিয়েট করে সেই ফলাফলে s=0,lambda=1,mu=2,a=2,b=3 বসিয়ে মান বের করা।
খুবই সহজ সমস্যা। কিন্তু আমি যে ফল পাচ্ছি (-.০০০০সামথিং) তা আশানুরুপ নয়।

তাই ভাবছি কেউ যদি (বিষেষত: স্বপ্নচারী ভাই) একটু চেক করে দেখতেন আমি কোন ভুল করছি কী-না। নিচের কোডগুলো ম্যাপলের জন্য। আশা করি বুঝতে অসুবিধা হবেনা।

ধন্যবাদ।

ebmus := exp(-b*(mu+s));
eamus := exp(-a*(mu+s));
ebs := exp(-b*s);
eas := exp(-a*s);
eamu := exp(-a *mu);
ebmu := exp(-b*mu);
esb := exp(-s*b);
lmus := lambda + mu +s;
ls := lambda+s;

e21 := L21 =    eamus * lambda * (L21/lmus) +
                               (eamus * mu * lambda * LA1)/(lmus * ls ) + eas * (1-eamu) * lambda* LA1/ls ;

ea1 := LA1 =    ebmus * (lambda /lmus) +
                               (ebmus * mu * lambda * L21/(lmus * lmus)) +
                               (ebmus * mu * mu * lambda * LA1) /(lmus* lmus* ls) +
                               (ebs * mu*b * ebmu * lambda * L21)/lmus +
                               (ebs * mu*b * esb * mu*lambda * LA1)/(lmus * ls) +
                               ebs * (1-ebmu - (mu*b * ebmu))*lambda * LA1 /ls;

ans := solve({e21, ea1},{LA1, L21});

df := diff(ans, s);

df0 :=subs(s=0,lambda=1,mu=2,a=2,b=3, df);

#simplify(df0);

evalf(df0);