المعطيات
المعادلة الأولى:dN1\dt=-c1N1
المعادلة الثانيه:dN2\dt=c1N1-c2N2
المعادلة الثالثه:dN3\dt=c2N2-c3N3
المعادلة الرابعه:dN4\dt=c3N3-c4N4
وباقي المعادلات نفس الطريقة حيث c1,c2.c3,c4 كلها ثوابت
الى أن نصل الى المعادلة العاشره:dN10\dt=c9N9-c10N10
والمعادلة الأخيره تكون على الشكل:dN11\dt=c10N10
ويعتبر ثابت التكامل الناتج عن عملية التكامل مساويا للصفر ماعدا N0 تكون له قيمه ولتكن 1000
المطلوب:
عمل برنامج الغرض منه حساب قيم N1,N2,N3,N4,N5,…و(في كل مره يتم ادخال قيم مختلفه لــN0 يعني تترك للمستخدم حرية ادخال الرقم وكذلك قيمة tولكن تكون على شكل مدى مثل [0 100] ). و في كل مره يتم رسم منحنى بياني يضم كل ثلاث معادلات مع بعض.
انا استخدمت طريقه ثانيه باستخدام m-file
function dN=new(t,N
a,b,c لها قيم ثابته
[(dN=[-aN(1);aN(1)-bN(2);bN(2)-cN(3);cN(3
ثم استدعائه في نافذة الأوامر من الأمر التالي:
([قيمة t],[قيمة الثابت ], t,N]=ode45(@new ]
وكتبتها بالنسبه لجميع المعادلات
لكن انا ابغى طريقه مايكون فيها أمر استدعاء
كل عام وأنتم بألف خير
اتمنى أن تردون علي بأسرع وقت ممكن وجزاكم الله خيرا
البرمجه التي عملها المهندس سعود مشكورا صحيحه ولكن ناقصه N11 والرسم كيف أقدر أكملها
أطلب اساسيات معادلات التفاضلية من الرتبة الاولى-معادلات تفاضلية العاديةفى عامل التكامل-"" “” ““الخطية العادية من الرتبة الثانية (المحددة)
”” " المتجانسة ذات المعاملات
"" “” غير متجانسة ذات المعاملات الثابتة -"" ""غير محددة