ট্র্যাপিজয়েডাল নিয়ম ক্যালকুলেটর

ট্র্যাপিজয়েডাল নিয়ম ক্যালকুলেটর সংখ্যাগতভাবে চাপের দৈর্ঘ্য অনুমান করে। এটি আপনার দেওয়া ফাংশনের চাপ-দৈর্ঘ্য ইন্টিগ্রালে কম্পোজিট ট্র্যাপিজয়েডাল নিয়ম প্রয়োগ করে। এক মিনিটের মধ্যে প্রথম অনুমান পেতে এই ধাপগুলো অনুসরণ করুন।

1. ফাংশনটি লিখুন। ফাংশন ঘরে f(x) টাইপ করুন, উদাহরণস্বরূপ x^2, sin(x), বা ln(x+1)। ঘাতের জন্য ^ এবং sin, cos, exp, ln এর মতো প্রমাণ ফাংশন নাম ব্যবহার করুন।
2. শুরুর মান লিখুন। নিম্ন সমাকলন সীমা a নির্ধারণ করুন। x = 0 থেকে x = 1 পর্যন্ত বক্ররেখা খণ্ডের জন্য 0 লিখুন।
3. শেষের মান লিখুন। ঊর্ধ্ব সমাকলন সীমা b নির্ধারণ করুন। ক্যালকুলেটরের জন্য b > a হওয়া আবশ্যক।
4. উপবিভাগ বেছে নিন। সমান উপ-ব্যবধানের সংখ্যা n লিখুন। মসৃণ বক্ররেখার জন্য n = 10 দিয়ে প্রথম রান ভালো কাজ করে।
5. ট্র্যাপিজয়েডাল পদ্ধতি নির্বাচন করুন। সংখ্যাসূচক সমাকলন ইন্টারফেসে ট্র্যাপিজয়েডাল নিয়ম বেছে নিন, যাতে নিচের প্যানেল সূত্রটিই প্রয়োগ হয়।
6. Calculate ক্লিক করুন। টুলটি n + 1টি নমুনা বিন্দুতে চাপ-দৈর্ঘ্য ইন্টিগ্র্যান্ড মূল্যায়ন করে, প্রান্তবিন্দুর ওজন প্রয়োগ করে এবং অনুমানটি দেখায়।
7. ফলাফল ব্যাখ্যা করুন। প্রদর্শিত সংখ্যাটি x = a থেকে x = b পর্যন্ত বক্ররেখার দৈর্ঘ্যের একটি অনুমান, আপনার x মানের একই দৈর্ঘ্য এককে।
8. নির্ভুলতা বাড়ান। n দ্বিগুণ করে পুনরায় গণনা করুন। পরপর দুটি অনুমান যখন প্রায় না বদলায়, তখন ব্যবহারিক উদ্দেশ্যে ফলাফল অভিসারী হয়েছে।

প্রতিটি ইনপুট কী নিয়ন্ত্রণ করে তা এখানে দেওয়া হলো:

f(x) = x², a = 0, b = 1 এবং n = 10 দিলে 1.479-এর কাছাকাছি চাপ-দৈর্ঘ্য অনুমান পাওয়া যায়, যা সঠিক মান 1.478943-এর কাছাকাছি। এই গাইডের বাকি অংশে ব্যাখ্যা করা হয়েছে ক্যালকুলেটর ঠিক কী গণনা করে এবং এর উত্তরের মান কীভাবে বিচার করবেন।

ট্র্যাপিজয়েডাল নিয়ম একটি সংখ্যাসূচক সমাকলন পদ্ধতি। এটি বক্ররেখার নিচের ক্ষেত্রফলকে ইন্টিগ্র্যান্ডের নিচে আঁকা ট্র্যাপিজয়েডের ক্ষেত্রফল দিয়ে প্রতিস্থাপন করে একটি নির্দিষ্ট ইন্টিগ্রাল অনুমান করে। প্রতিটি ট্র্যাপিজয়েড দুটি পাশাপাশি নমুনা বিন্দু থেকে তৈরি হয়, এবং সব ট্র্যাপিজয়েডের ক্ষেত্রফলের যোগফল ইন্টিগ্রালটির অনুমান দেয়।

নামটি এসেছে জ্যামিতি থেকে। দুটি নমুনা বিন্দু x_i এবং x_i+1-এর মধ্যে ইন্টিগ্র্যান্ড g(x)-কে একটি সরলরেখা দিয়ে প্রতিস্থাপন করা হয়। সেই সরলরেখার নিচের অঞ্চলটি একটি ট্র্যাপিজয়েড: g(x_i) ও g(x_i+1) উচ্চতার দুটি সমান্তরাল উল্লম্ব বাহু এবং h = x_i+1 - x_i অনুভূমিক প্রস্থ।

একটি ট্র্যাপিজয়েডের ক্ষেত্রফল h/2 · [g(x_i) + g(x_i+1)]। সব nটি ট্র্যাপিজয়েডের ক্ষেত্রফল যোগ করলে ইন্টিগ্রালের কম্পোজিট ট্র্যাপিজয়েডাল অনুমান পাওয়া যায়। "কম্পোজিট" শব্দের অর্থ হলো ব্যবধানটিকে একটির বদলে কয়েকটি প্যানেলে ভাগ করা হয়েছে।

এই নিয়মটি সংখ্যাসূচক কোয়াড্রেচার পদ্ধতির পরিবারের অন্তর্গত। এটি স্বচ্ছতার জন্য মূল্যবান: প্রতিটি ধাপের একটি স্পষ্ট জ্যামিতিক অর্থ আছে, এবং একমাত্র প্যারামিটার হলো উপবিভাগের সংখ্যা n। মূল ধারণাটি সহজে বলা যায়: সরলতর, সংকীর্ণতর প্যানেল ইন্টিগ্র্যান্ডকে আরও ঘনিষ্ঠভাবে অনুসরণ করে, তাই n বাড়লে অনুমান উন্নত হয়।

চাপের দৈর্ঘ্যের ক্ষেত্রে ইন্টিগ্র্যান্ড মূল ফাংশন f(x) নয়। এটি চাপ-দৈর্ঘ্য ইন্টিগ্র্যান্ড g(x) = √(1 + (f'(x))²), যা বক্ররেখা বরাবর দৈর্ঘ্য কত দ্রুত জমা হয় তা পরিমাপ করে। এই পার্থক্য স্পষ্ট রাখলে পদ্ধতিটির সবচেয়ে সাধারণ ভুল বোঝাবুঝি এড়ানো যায়।

ট্র্যাপিজয়েডাল নিয়মের সূত্রের দুটি রূপ আছে: একটি একক-প্যানেল সংস্করণ এবং বাস্তবে ব্যবহৃত কম্পোজিট সংস্করণ। দুটোই একই ট্র্যাপিজয়েড জ্যামিতি থেকে আসে।

একক-প্যানেল সূত্র। পুরো ব্যবধান জুড়ে একটি ট্র্যাপিজয়েড দিয়ে:

প্রস্থ (b - a)-কে দুটি প্রান্তবিন্দুর উচ্চতার গড় দিয়ে গুণ করা হয়। এটি ঠিক একটি ট্র্যাপিজয়েডের ক্ষেত্রফল।

কম্পোজিট সূত্র। [a, b]-কে h = (b - a)/n প্রস্থের nটি সমান উপ-ব্যবধানে ভাগ করুন, যেখানে নমুনা বিন্দু x₀ = a, x₁ = a + h, ..., x_n = b:

প্রতিটি প্রতীকের একটি নির্দিষ্ট অর্থ আছে:

• T_n হলো nটি প্যানেলসহ সংখ্যাসূচক অনুমান।
• h হলো উপ-ব্যবধানের প্রস্থ, (b - a)/n।
• a ও b হলো নিম্ন ও ঊর্ধ্ব সমাকলন সীমা।
• n হলো উপ-ব্যবধানের সংখ্যা, একটি ধনাত্মক পূর্ণসংখ্যা।
• x_i হলো নমুনা বিন্দুগুলো, h দূরত্বে ব্যবধানযুক্ত।
• g(x_i) হলো সেই বিন্দুগুলোতে ইন্টিগ্র্যান্ডের মান।

প্রান্তবিন্দুগুলো অর্ধেক ওজন পায় কারণ প্রতিটি ভেতরের নমুনা বিন্দু দুটি পাশাপাশি ট্র্যাপিজয়েডের মধ্যে ভাগ হয়, অথচ প্রতিটি প্রান্তবিন্দু কেবল একটির অন্তর্গত। বীজগাণিতিকভাবে, পাশাপাশি h/2 · [g(x_i) + g(x_i+1)] পদগুলো এমনভাবে মিলে যায় যে ভেতরের মানগুলো পূর্ণ গণনা পায় এবং দুটি বাইরের মান অর্ধেক গণনা পায়।

একটি সংক্ষিপ্ত হাতে-কলমে গণনা। n = 2 দিয়ে ∫₀¹ x² dx অনুমান করুন। তাহলে h = 0.5, x₀ = 0, x₁ = 0.5, x₂ = 1, এবং T₂ = 0.5 · [0/2 + 0.25 + 1/2] = 0.5 · 0.75 = 0.375। সঠিক ইন্টিগ্রাল হলো 1/3 ≈ 0.333, তাই মোটা দুই-প্যানেল অনুমানটি সঠিক এলাকায় আছে, পরিশোধনের জায়গা রেখে।

চাপের দৈর্ঘ্য একটি ইন্টিগ্রাল, তাই ট্র্যাপিজয়েডাল নিয়ম সরাসরি এতে প্রযোজ্য। x = a থেকে x = b পর্যন্ত y = f(x) বক্ররেখার চাপের দৈর্ঘ্য হলো:

ক্যালকুলেটরটি চাপ-দৈর্ঘ্য ইন্টিগ্র্যান্ড g(x) = √(1 + (f'(x))²) গঠন করে এবং কম্পোজিট ট্র্যাপিজয়েডাল সূত্রে ঢুকিয়ে দেয়। সম্পূর্ণ প্রক্রিয়াটি হলো:

1. y = f(x) দিয়ে শুরু করুন। নিশ্চিত করুন ফাংশনটি [a, b]-তে অন্তরকলনযোগ্য।
2. অন্তরজ f'(x) নির্ণয় করুন। টুলটি ভেতরে ভেতরে দেওয়া রাশিটির অন্তরকলন করে।
3. অন্তরজের বর্গ করুন। (f'(x))² গঠন করুন। বর্গ করলে চিহ্নের তথ্য মুছে যায় কারণ ঢালের দিক যাই হোক দৈর্ঘ্য জমা হয়।
4. 1 যোগ করুন। ইন্টিগ্র্যান্ড হলো 1 + (f'(x))², যা কখনো 1-এর চেয়ে ছোট নয়।
5. বর্গমূল নিন। এতে g(x) পাওয়া যায়, x-এর প্রতি এককে দৈর্ঘ্য জমার স্থানীয় হার।
6. ব্যবধান ভাগ করুন। n বেছে নিন এবং h = (b - a)/n গণনা করুন।
7. ইন্টিগ্র্যান্ড মূল্যায়ন করুন। নমুনা বিন্দুগুলোতে g(x₀), ..., g(x_n) গণনা করুন।
8. কম্পোজিট নিয়ম প্রয়োগ করুন। প্রান্তবিন্দুতে অর্ধেক ওজন দিয়ে মানগুলো একত্র করুন।
9. ফলাফল ব্যাখ্যা করুন। যোগফলটি বক্ররেখা খণ্ডের দৈর্ঘ্যের আনুমানিক মান।

অন্তরজের বর্গ করে তারপর বর্গমূলের নিচে বসানোর একটি জ্যামিতিক কারণ আছে: dx প্রস্থ ও dy উচ্চতার একটি ক্ষুদ্র বক্ররেখা খণ্ডের দৈর্ঘ্য পিথাগোরাসের উপপাদ্য অনুযায়ী √(dx² + dy²) = √(1 + (dy/dx)²) dx। [a, b]-তে সেই স্থানীয় দৈর্ঘ্য সমাকলন করলে মোট দৈর্ঘ্য পাওয়া যায়, এবং ট্র্যাপিজয়েডাল নিয়ম ইন্টিগ্রালটি অনুমান করে।

ইন্টিগ্র্যান্ড সবসময় অঋণাত্মক কারণ 1 + (f'(x))² কমপক্ষে 1। এর সর্বনিম্ন মান 1, যেখানে বক্ররেখা সমতল সেখানে পৌঁছায়। খাড়া ঢাল g(x)-কে বড় করে, আর ঠিক এই কারণেই একই x-ব্যবধানে খাড়া বক্ররেখা অগভীর বক্ররেখার চেয়ে লম্বা হয়।

এই সম্পূর্ণ হাতে-কলমে উদাহরণে ট্র্যাপিজয়েডাল নিয়মের প্রতিটি সংখ্যা দেখানো হয়েছে। [0, 1]-তে f(x) = x² ব্যবহার করুন, n = 4 সহ।

1. অন্তরজ। f'(x) = 2x।

2. ইন্টিগ্র্যান্ড। g(x) = √(1 + (2x)²) = √(1 + 4x²)।

3. ধাপের আকার। h = (1 - 0)/4 = 0.25।

4. নমুনা বিন্দু। 0, 0.25, 0.5, 0.75, 1।

5. ইন্টিগ্র্যান্ড মান ও ওজন।

6. ওজনযুক্ত যোগফল। 1.000000/2 + 1.118034 + 1.414214 + 1.802776 + 2.236068/2 = 5.953057।

7. h দিয়ে গুণ করুন।

8. সঠিক মানের সাথে তুলনা করুন। সঠিক চাপের দৈর্ঘ্য হলো L = √5/2 + asinh(2)/4 ≈ 1.478943। n = 4 অনুমান 1.488264 এটির চেয়ে 0.009321 বেশি, প্রায় 0.63% ত্রুটি। n = 8 দিয়ে পুনরাবৃত্তি করলে 1.481272 পাওয়া যায়, যা ইতিমধ্যে সঠিক মান থেকে 0.002330-এর মধ্যে।

লক্ষ্য করুন অনুমানটি সঠিক মানের চেয়ে বড়। এই উত্তল ইন্টিগ্র্যান্ডের ক্ষেত্রে সরলরেখা জ্যাগুলো বক্ররেখার ওপরে বসে, তাই ট্র্যাপিজয়েডাল নিয়ম অতিরিক্ত অনুমান করে। এই নিয়মতান্ত্রিক আচরণ স্বাভাবিক, এবং n বাড়লে অতিরিক্ত অনুমান কমে যায়।

উপবিভাগের সংখ্যা n এবং ধাপের আকার h, h = (b - a)/n সম্পর্ক দিয়ে যুক্ত। দুটো মিলে নিয়ন্ত্রণ করে ইন্টিগ্র্যান্ড কত সূক্ষ্মভাবে নমুনা নেওয়া হয় এবং অনুমান কতটা নির্ভুল।

n কী বোঝায়। এটি ট্র্যাপিজয়েড প্যানেলের সংখ্যা, একটি ধনাত্মক পূর্ণসংখ্যা। নমুনা বিন্দুর সংখ্যা সবসময় n + 1 কারণ পরপর প্যানেল তাদের সীমানা বিন্দু ভাগ করে।

n = 1 কেন একক-প্যানেল নিয়ম। n = 1 হলে h = b - a এবং T₁ = h/2 · [g(a) + g(b)]। এটি পদ্ধতির সবচেয়ে মোটা রূপ এবং খুব কমই যথেষ্ট নির্ভুল।

[0, 1]-তে প্যানেল প্রস্থ। নিচের সারণিতে দেখানো হয়েছে n বাড়লে h কীভাবে কমে:

ছোট প্যানেল ইন্টিগ্র্যান্ডকে আরও ঘনিষ্ঠভাবে অনুসরণ করে, তাই প্রতিটি প্যানেলের ভেতরে সরলরেখা অনুমান প্রকৃত বক্ররেখা থেকে কম বিচ্যুত হয়। এই কারণে n বাড়ালে সাধারণত নির্ভুলতা বাড়ে। বিনিময় হলো গণনার খরচ: প্রতিটি নতুন প্যানেলের জন্য আরও একটি ফাংশন মূল্যায়ন দরকার।

বেশি উপবিভাগ প্রতিটি ইন্টিগ্র্যান্ডের জন্য কঠোরভাবে ভালো ফলাফলের নিশ্চয়তা দেয় না। অমসৃণ বা দোদুল্যমান ইন্টিগ্র্যান্ড অনিয়মিত ত্রুটি আচরণ দেখাতে পারে, এবং খুব বড় n-এ ভাসমান-বিন্দু গোলকরণ শেষ পর্যন্ত লাভ সীমিত করে। তবে এই গাইডের মতো মসৃণ বক্ররেখার ক্ষেত্রে উন্নতি স্থিতিশীল ও অনুমানযোগ্য।

একটি ব্যবহারিক কৌশল: n = 10 বা n = 20 দিয়ে শুরু করুন, তারপর আপনার প্রয়োজনীয় নির্ভুলতায় অনুমান বদলানো বন্ধ না হওয়া পর্যন্ত n দ্বিগুণ করুন। n-এর যেকোনো একক মানের ওপর আস্থা রাখার চেয়ে এই অভিসারণ পরীক্ষা বেশি নির্ভরযোগ্য।

ছয়টি উদাহরণ সঠিক ক্ষেত্র থেকে এমন ফাংশন পর্যন্ত বিস্তৃত যাদের চাপের দৈর্ঘ্যের কোনো প্রাথমিক বদ্ধ রূপ নেই। প্রতিটিতে অন্তরজ, ইন্টিগ্র্যান্ড, প্রতিস্থাপন এবং সংখ্যাসূচক ফলাফল দেখানো হয়েছে।

উদাহরণ 1: সরলরেখা, [0, 3]-তে f(x) = 2x। অন্তরজ হলো f'(x) = 2, একটি ধ্রুবক, তাই ইন্টিগ্র্যান্ড হলো g(x) = √(1 + 4) = √5, তাও ধ্রুবক। একটি প্যানেলই সঠিক কারণ ইন্টিগ্র্যান্ড কখনো বাঁকে না: T₁ = 3 · √5 ≈ 6.708204। সঠিক দৈর্ঘ্য 3√5 মিলে যায়। এতে নিয়মটি নিশ্চিত হয়: ইন্টিগ্র্যান্ড রৈখিক বা ধ্রুবক হলে ট্র্যাপিজয়েডাল নিয়ম সঠিক।

উদাহরণ 2: দ্বিঘাত, [0, 1]-তে f(x) = x², n = 4। f'(x) = 2x এবং g(x) = √(1 + 4x²) সহ, আগের অংশের সম্পূর্ণ গণনায় L ≈ 1.488264 পাওয়া যায়, সঠিক 1.478943-এর তুলনায়। ছোট ব্যবধানটি আসে উত্তল ইন্টিগ্র্যান্ড থেকে; বক্ররেখার ওপরে টানা জ্যাগুলো অতিরিক্ত ক্ষেত্রফল যোগ করে।

উদাহরণ 3: দীর্ঘতর দ্বিঘাত, [0, 2]-তে f(x) = x², n = 8। এখানে h = 0.25 আবার, কিন্তু ব্যবধান দ্বিগুণ লম্বা এবং ইন্টিগ্র্যান্ড খাড়া। অনুমান হলো T₈ = 4.656890। উচ্চ-রেজোলিউশন গণনায় 4.646784 পাওয়া যায়, তাই ত্রুটি প্রায় 0.010106। ব্যবধান বাড়ালে উত্তর শুধু স্কেল হয় না: ইন্টিগ্র্যান্ড x-এর সাথে বাড়ে, তাই চাপের দৈর্ঘ্য ব্যবধানের দৈর্ঘ্যের চেয়ে দ্রুত বাড়ে।

উদাহরণ 4: সাইন বক্ররেখা, [0, π]-তে f(x) = sin x। অন্তরজ হলো f'(x) = cos x, ফলে g(x) = √(1 + cos²x)। n = 8-এ অনুমান 3.820198; n = 16-এও 3.820198, ছয় দশমিক স্থানে অপরিবর্তিত। প্রতিসম ইন্টিগ্র্যান্ড দ্রুত অভিসারী হয় কারণ এর বক্রতা মৃদু ও সমানভাবে বণ্টিত। (কোণ রেডিয়ানে লিখুন; ডিগ্রি দিলে ভিন্ন বক্ররেখা হবে।)

উদাহরণ 5: লগারিদম, [0, 2]-তে f(x) = ln(x + 1), n = 10। অন্তরজ হলো f'(x) = 1/(x + 1), তাই g(x) = √(1 + 1/(x + 1)²)। ডোমেনের জন্য x > -1 দরকার, যা ব্যবধানটি মেনে চলে। অনুমান হলো T₁₀ = 2.304215, অথচ পরিশোধিত গণনায় 2.301988, প্রায় 0.002227 ত্রুটি। x = 0-এর কাছে ঢাল সবচেয়ে খাড়া, তাই ব্যবধানের শুরুর দিকে বেশিরভাগ দৈর্ঘ্য জমা হয়।

উদাহরণ 6: ঘন, [0, 1]-তে f(x) = x³, n = 10। এখানে f'(x) = 3x² এবং g(x) = √(1 + 9x⁴)। এই চাপের দৈর্ঘ্যের কোনো প্রাথমিক প্রতিঅন্তরক নেই, আর ঠিক তখনই সংখ্যাসূচক সমাকলন তার জায়গা করে নেয়। অনুমান হলো T₁₀ = 1.552609, পরিশোধিত মান 1.547866-এর বিপরীতে, প্রায় 0.004743 ত্রুটি। সংখ্যাসূচক পদ্ধতি ছাড়া এই দৈর্ঘ্য পাওয়া খুব কঠিন হতো।

উপরের প্রতিটি ফলাফল কম্পোজিট সূত্র থেকে সরাসরি গণনা করা হয়েছে, ছয় দশমিক স্থানে গোল করে। ক্যালকুলেটরে একই n দিয়ে এগুলোর যেকোনোটি পুনরাবৃত্তি করলে একই অঙ্ক পাওয়া যায়।

অভিসারণ মানে n বাড়লে অনুমানগুলো প্রকৃত মানের দিকে স্থির হয়। নিচের সারণিতে [0, 1]-তে f(x) = x² অনুসরণ করা হয়েছে, যার সঠিক দৈর্ঘ্য 1.478943:

n দ্বিগুণ হলে প্রতিবার ত্রুটি প্রায় চার গুণ কমে: 0.037181 থেকে 0.009321 থেকে 0.002330 থেকে 0.000582 থেকে 0.000146 থেকে 0.000036। এটি দ্বিতীয়-ক্রমের অভিসারণ, মসৃণ ইন্টিগ্র্যান্ডে কম্পোজিট ট্র্যাপিজয়েডাল নিয়মের স্বাক্ষর। h অর্ধেক করলে প্রধান ত্রুটি পদ প্রায় এক-চতুর্থাংশে নেমে আসে।

"আগেরটি থেকে পরিবর্তন" কলামটি হলো ব্যবহারিক অভিসারণ পরীক্ষা যা ক্যালকুলেটর কর্মপ্রবাহ সুপারিশ করে। n দ্বিগুণ করলে উত্তর আপনার প্রয়োজনীয় নির্ভুলতার চেয়ে কম বদলালে অনুমান স্থিতিশীল। এই বক্ররেখার জন্য n = 16 ইতিমধ্যে দৈর্ঘ্যকে প্রকৃত মান থেকে 0.000582-এর মধ্যে আটকে দেয়।

কম্পোজিট ট্র্যাপিজয়েডাল নিয়মের ত্রুটির একটি সুনির্দিষ্ট গাণিতিক রূপ আছে। [a, b]-তে পর্যাপ্ত মসৃণ ইন্টিগ্র্যান্ড g-এর জন্য:

[a, b]-তে কোনো ξ-এর জন্য। চিহ্নটি গুরুত্বপূর্ণ: g'' ধনাত্মক হলে (উত্তল ইন্টিগ্র্যান্ড) ত্রুটি ঋণাত্মক হয়, অর্থাৎ T_n ইন্টিগ্রালকে অতিরিক্ত অনুমান করে। এটি x² উদাহরণের সাথে মেলে, যেখানে প্রতিটি অনুমান সঠিক 1.478943-কে ছাড়িয়ে গিয়েছিল।

g''(ξ)-কে [a, b]-তে এর সর্বোচ্চ মান M = max|g''(x)| দিয়ে প্রতিস্থাপন করলে একটি নিশ্চিত সীমা পাওয়া যায়:

মূল পদ। পরম ত্রুটি হলো |I - T_n|, প্রকৃত মান ও অনুমানের মধ্যে দূরত্ব। আপেক্ষিক ত্রুটি এটিকে |I| দিয়ে ভাগ করে। ছাঁটাই ত্রুটি হলো ইন্টিগ্রালকে সসীম যোগফল দিয়ে অনুমান করার ত্রুটি, সসীম-নির্ভুলতা পাটিগণিতের গোলকরণ ত্রুটির বিপরীতে। দ্বিতীয় অন্তরজ g'' ইন্টিগ্র্যান্ডের বক্রতা পরিমাপ করে, যা সরলরেখা প্যানেল ধরতে ব্যর্থ হয়।

x² উদাহরণের জন্য সুনির্দিষ্ট সীমা। [0, 1]-তে g''(x) = 4/(1 + 4x²)^3/2, x = 0-তে M = 4 সহ সবচেয়ে বড়। সীমা দেয় |I - T₄| ≤ 4/(12 · 16) = 0.020833। প্রকৃত ত্রুটি 0.009321 আরামে এর ভেতরে। সীমা নকশা অনুযায়ী রক্ষণশীল; এটি একটি ছাদ প্রতিশ্রুতি দেয়, সঠিক ত্রুটি নয়।

দুটি সতর্কতা। প্রথমত, g হলো চাপ-দৈর্ঘ্য ইন্টিগ্র্যান্ড, মূল f নয়। বর্গমূলের মধ্য দিয়ে দুবার অন্তরকলন করলে যা পাওয়া যায়, M সেটাই বর্ণনা করে। দ্বিতীয়ত, জটিল ফাংশনের জন্য ত্রুটি সীমা গণনা কষ্টকর হতে পারে কারণ M-এর নিজেরই সর্বোচ্চকরণ দরকার হতে পারে। বাস্তবে, n দ্বিগুণ করার অভিসারণ পরীক্ষাই বেশি উপযোগী টুল, এবং ক্যালকুলেটর স্বয়ংক্রিয় ত্রুটি সীমা দেখায় না।

সিম্পসনের নিয়ম ট্র্যাপিজয়েডাল নিয়মের নিকটতম আত্মীয়। দুটোই [a, b]-কে সমান উপ-ব্যবধানে ভাগ করে, কিন্তু ট্র্যাপিজয়েডাল নিয়ম দুটি বিন্দুর মধ্য দিয়ে একটি সরল জ্যা টানে যেখানে সিম্পসন পরপর তিনটি বিন্দুর মধ্য দিয়ে একটি পরাবৃত্ত বসায়।

মসৃণ ইন্টিগ্র্যান্ডে সিম্পসনের চতুর্থ-ক্রমের অভিসারণ সাধারণত জেতে: n দ্বিগুণ করলে এর ত্রুটি প্রায় ষোল গুণ কমে, ট্র্যাপিজয়েডাল নিয়মের চার গুণের বিপরীতে। কিন্তু সিম্পসন সার্বজনীনভাবে ভালো নয়। কর্কশ বা বিচ্ছিন্ন ইন্টিগ্র্যান্ডের ক্ষেত্রে দুটি পদ্ধতিই খারাপ হয়, এবং অতিরিক্ত নির্ভুলতা ক্রম বাস্তবায়িত নাও হতে পারে। ট্র্যাপিজয়েডাল নিয়ম বিজোড় n-ও মেনে নেয়, যা দ্রুত পরীক্ষার জন্য সুবিধাজনক।

একটি ভালো কর্মপ্রবাহ দুটোই ব্যবহার করে: স্বচ্ছ ভিত্তিরেখা হিসেবে ট্র্যাপিজয়েডাল অনুমান, এবং বেশি নির্ভুলতার জন্য সিম্পসনের নিয়ম। দুটো আপনার প্রয়োজনীয় অঙ্কে মিলে গেলে আস্থা বেশি। না মিললে কোনোটিকে বিশ্বাস করার আগে n বাড়ান।

দুটি ভিন্ন "ট্র্যাপিজয়েড-সদৃশ" ধারণা সহজে গুলিয়ে যায়। একই বক্ররেখার জন্য এগুলো ভিন্ন সংখ্যা দেয়, তাই পার্থক্যটি গুরুত্বপূর্ণ।

পদ্ধতি 1: চাপ-দৈর্ঘ্য ইন্টিগ্র্যান্ডের ট্র্যাপিজয়েডাল সমাকলন। ক্যালকুলেটর এটাই করে। এটি ইন্টিগ্র্যান্ড বক্ররেখার নিচে ট্র্যাপিজয়েড প্যানেল দিয়ে g(x) = √(1 + (f'(x))²) সমাকলন করে:

পদ্ধতি 2: মূল বক্ররেখার বহুভুজীয় (জ্যা-যোগফল) অনুমান। এটি বক্ররেখার ওপর পাশাপাশি বিন্দুগুলোকে সরল খণ্ড দিয়ে যুক্ত করে এবং তাদের দৈর্ঘ্য যোগ করে:

জ্যা-যোগফল মূল বক্ররেখার ওপর (x_i, f(x_i)) বিন্দুগুলো যুক্ত করে। ট্র্যাপিজয়েডাল প্যানেল বরং ইন্টিগ্র্যান্ড g(x)-এর নিচের ক্ষেত্রফল অনুমান করে। দুটি ছবিতেই সরলরেখা থাকলেও প্যানেলগুলো মূল বক্ররেখার সরল খণ্ড উপস্থাপন করে না।

n বাড়লে দুটি পদ্ধতিই প্রকৃত দৈর্ঘ্যের দিকে অভিসারী হয়, এবং দুটোই বৈধ সংখ্যাসূচক টুল। এই পৃষ্ঠার ক্যালকুলেটর পদ্ধতি 1 বাস্তবায়ন করে। স্থানাঙ্ক তথ্য থেকে জ্যা-যোগফল দরকার হলে বিন্দু থেকে চাপের দৈর্ঘ্য টুলটি বরং পদ্ধতি 2 অনুসরণ করে।

নির্ভুলতার শেষ অঙ্কগুলো চেপে বের করার চেয়ে সরলতা ও স্বচ্ছতা বেশি গুরুত্বপূর্ণ হলে ট্র্যাপিজয়েডাল নিয়মই সঠিক পছন্দ।

• সংখ্যাসূচক সমাকলন শেখা। প্রতিটি ধাপের একটি দৃশ্যমান জ্যামিতিক অর্থ আছে, যা এটিকে ক্যালকুলাস কোর্সে প্রমাণ প্রথম পদ্ধতি করে।
• বক্ররেখার দৈর্ঘ্য দ্রুত অনুমান করা। মুষ্টিমেয় প্যানেল মসৃণ বক্ররেখার জন্য প্রায়শই দুই বা তিনটি সঠিক অঙ্ক দেয়।
• বিশ্লেষণমূলক সমাধান পরীক্ষা করা। হাতে বের করা চাপের দৈর্ঘ্যকে একটি স্বাধীন সংখ্যাসূচক অনুমানের বিপরীতে তুলনা করুন।
• প্রকৌশল ও পদার্থবিজ্ঞান মডেল। তার, খিলান ও গতিপথের দৈর্ঘ্য অনুমান যেখানে আনুমানিক মান বড় গণনায় ঢোকে।
• সংখ্যাসূচক পদ্ধতির তুলনা। এটি সেই ভিত্তিরেখা যার বিপরীতে সিম্পসনের নিয়ম ও উচ্চতর-ক্রমের পদ্ধতি বিচার করা হয়।
• মসৃণ নমুনাযুক্ত ফাংশন। f মসৃণ হলে এবং n মাঝারি রাখা গেলে পদ্ধতিটি নির্ভরযোগ্য।
• শিক্ষামূলক যাচাই। এর অনুমানযোগ্য O(h²) আচরণ অভিসারণ পরীক্ষা ব্যাখ্যা সহজ করে।

ইন্টিগ্র্যান্ড দ্রুত বদলালে, ভাঁজ বা বিচ্ছিন্নতা থাকলে, বা কম খরচে উচ্চ নির্ভুলতা দরকার হলে ভিন্ন পদ্ধতি নিন। মসৃণ বক্ররেখার জন্য সিম্পসনের নিয়ম সাধারণত প্রতি প্যানেলে বেশি অঙ্ক দেয়। প্যারামেট্রিক, পোলার বা ত্রিমাত্রিক বক্ররেখার জন্য নিবেদিত প্যারামেট্রিক, পোলার বা 3D চাপের দৈর্ঘ্য ক্যালকুলেটর ব্যবহার করুন, যা একই সংখ্যাসূচক ধারণা সঠিক স্থানাঙ্ক পদ্ধতিতে প্রয়োগ করে।

ছাত্রদের কাজে ও সহায়তা প্রশ্নে এগুলোই সবচেয়ে বেশি দেখা যায়। একবার জানা থাকলে প্রতিটি এড়ানো সহজ।

1. ভুল সমাকলন সীমা লেখা। [0, 2] খণ্ড [2, 0] খণ্ডের মতো একই বক্ররেখা অংশ নয়। সবসময় a < b নিশ্চিত করুন এবং সীমা সমস্যার সাথে মিলিয়ে নিন।
2. ভুল ফাংশন ব্যবহার করা। x^2-এর বদলে x^3-এর মতো টাইপো অন্তরজ এবং পরের প্রতিটি মান বদলে দেয়।
3. চাপ-দৈর্ঘ্য সূত্রে অন্তরজ ভুলে যাওয়া। চাপের দৈর্ঘ্য √(1 + (f')²) সমাকলন করে, কখনো f নিজেকে নয়।
4. চাপের দৈর্ঘ্যের জন্য নিয়মটি সরাসরি f(x)-তে প্রয়োগ করা। প্যানেলগুলোকে g(x)-এর নিচে বসতে হবে। f(x) সমাকলন করলে চিহ্নযুক্ত ক্ষেত্রফল গণনা হয়, দৈর্ঘ্য নয়।
5. খুব কম উপবিভাগ বেছে নেওয়া। বাঁকা ইন্টিগ্র্যান্ডে n = 2 কয়েক শতাংশ ভুল হতে পারে; x² উদাহরণে n = 2-এ 2.5% ত্রুটি ছিল।
6. শূন্য উপবিভাগ ব্যবহার করা। n = 0 মানে h = (b - a)/0, যা অসংজ্ঞায়িত। n-কে ধনাত্মক পূর্ণসংখ্যা হতেই হবে।
7. অপূর্ণসংখ্যক উপবিভাগ সংখ্যা লেখা। n = 2.5-এর প্যানেল গণনায় কোনো অর্থ নেই; পূর্ণ সংখ্যায় গোল করুন।
8. n-কে নমুনা বিন্দুর সংখ্যার সাথে গুলিয়ে ফেলা। nটি প্যানেলের জন্য n + 1টি নমুনা বিন্দু দরকার। n = 4 হলে বিন্দু 5টি, 4টি নয়।
9. প্রান্তবিন্দুর অর্ধেক ওজন ভুলে যাওয়া। g(x₀) ও g(x_n)-কে পূর্ণ ওজন দিলে অনুমান h/2 · (g(x₀) + g(x_n)) পরিমাণ ফুলে যায়।
10. ভুল ধাপের আকার ব্যবহার করা। h = (b - a)/n, b - a নয় এবং 1/n নয়। [0, 2]-তে n = 8 হলে h = 0.25।
11. ডিগ্রি ও রেডিয়ান মেশানো। ডিগ্রিতে x সহ sin(x) রেডিয়ানে x সহ sin(x)-এর চেয়ে ভিন্ন বক্ররেখা; চাপের দৈর্ঘ্যও তদনুযায়ী ভিন্ন হয়।
12. ফাংশনের ডোমেন সীমাবদ্ধতা উপেক্ষা করা। ln(x + 1)-এর জন্য x > -1 দরকার। ডোমেনের বাইরে মূল্যায়ন করলে অর্থহীন মান আসে।
13. অনুমানকে সঠিক সমাকলনের সাথে গুলিয়ে ফেলা। T_n L-এর দিকে এগোয়; এটি L-এর সমান হয় কেবল ধ্রুবকের মতো বিশেষ ইন্টিগ্র্যান্ডে।
14. ধরে নেওয়া বেশি প্যানেল সবসময় ভালো নির্ভুলতার নিশ্চয়তা দেয়। মসৃণ ইন্টিগ্র্যান্ডের উপসর্গীয় অঞ্চলে সত্য, কিন্তু চরম পর্যায়ে গোলকরণ ত্রুটি ও অমসৃণতা ধারা ভেঙে দিতে পারে।
15. অভিসারণ পরীক্ষা উপেক্ষা করা। একটি অনুমান হলো আন্দাজ; n ও 2n-এ মিলে যাওয়া দুটি অনুমান হলো প্রমাণ।
16. সিম্পসনের নিয়মকে ট্র্যাপিজয়েডাল নিয়মের সাথে গুলিয়ে ফেলা। সিম্পসন 1-4-2-4-...-1 ওজন ব্যবহার করে এবং জোড় n দরকার; ট্র্যাপিজয়েডাল নিয়ম 1/2-1-...-1-1/2 ব্যবহার করে এবং যেকোনো ধনাত্মক n চলে।
17. ট্র্যাপিজয়েডাল সমাকলনকে জ্যা-যোগফলের সাথে গুলিয়ে ফেলা। ইন্টিগ্র্যান্ডের নিচের প্যানেল বক্ররেখার ওপর জ্যা নয়; দুটি পদ্ধতি ভিন্ন সংখ্যা দেয়।
18. অতিরিক্ত দশমিক নির্ভুলতা জানানো। ত্রুটি 0.002-এর কাছাকাছি হলে ছয় দশমিক স্থান লেখা জানা তথ্যকে বাড়িয়ে বলে। যাচাইকৃত নির্ভুলতার সাথে মিলিয়ে গোল করুন।

চাপের দৈর্ঘ্যের পেছনের একই প্যানেল পাটিগণিত সেখানেই দেখা যায় যেখানে একটি মসৃণ রাশিকে সংখ্যাগতভাবে জমা করতে হয়।

• প্রকৌশল বক্ররেখা-দৈর্ঘ্য অনুমান। তৈরির আগে নকশা ফাংশন থেকে তার, ঝুলন্ত বক্ররেখা ও পাইপ বাঁকের দৈর্ঘ্য গণনা করা হয়।
• যান্ত্রিক নকশা। ক্যাম প্রোফাইল ও গিয়ার দাঁতের বক্ররেখার উপাদান পরিকল্পনার জন্য নির্ভুল পরিধি অনুমান দরকার।
• সড়ক ও পথের জ্যামিতি। সড়ক ও রেলপথের সারিবদ্ধ বক্ররেখা তাদের কেন্দ্ররেখা বরাবর মাপা হয়, যা খুব কমই প্রাথমিক ফাংশন।
• কম্পিউটার-সহায়ত নকশা। রেন্ডারিং, টুলপাথ তৈরি ও মাত্রা নির্ধারণের জন্য CAD কার্নেল বক্ররেখার দৈর্ঘ্য অনুমান করে।
• পদার্থবিজ্ঞান সিমুলেশন। সংখ্যাসূচক মডেল থেকে কণার গতিপথ ও ক্ষেত্ররেখার দৈর্ঘ্য জমা করা হয়।
• সংখ্যাসূচক মডেলিং। সম্ভাবনা ঘনত্ব থেকে সংকেত শক্তি পর্যন্ত বদ্ধ রূপহীন যেকোনো ইন্টিগ্রাল ট্র্যাপিজয়েডাল প্যানেল দিয়ে শুরু হতে পারে।
• গাণিতিক শিক্ষা। নিয়মটি সংখ্যাসূচক বিশ্লেষণের প্রমাণ প্রবেশদ্বার কারণ এর ত্রুটি আচরণ প্রদর্শন সহজ।
• বক্ররেখা পরিমাপ ও যাচাই। অনিয়মিত বক্ররেখার ভৌত পরিমাপ সংখ্যাসূচক পূর্বাভাসের বিপরীতে ক্রস-চেক করা হয়।

বিশেষ জ্যামিতির জন্য পদ্ধতিটি খাপ খায়: একই কম্পোজিট নিয়ম প্যারামেট্রিক গতি √((x')² + (y')²) ও পোলার ইন্টিগ্র্যান্ড সমাকলন করে। এখানকার ক্যালকুলেটর কার্তেসীয় ক্ষেত্র সামলায়; সম্পর্কিত টুলগুলো অন্য স্থানাঙ্ক পদ্ধতি কভার করে। এটি CAD ফাইল বা জরিপ ডেটাসেট সরাসরি প্রক্রিয়া করে না; সেগুলোকে আগে একটি ফাংশন ও সীমায় নামিয়ে আনতে হয়।

এগুলো হাতে করুন, তারপর ক্যালকুলেটরের সাথে মিলিয়ে নিন। প্রতিটি সমাধানে প্রতিস্থাপন দেখানো হয়েছে যাতে নিজের পাটিগণিতের ভুল খুঁজে পান।

সমস্যা 1. n = 1 সহ ট্র্যাপিজয়েডাল নিয়ম ব্যবহার করে ∫₀² 3 dx অনুমান করুন, এবং ফলাফল সঠিক কিনা বলুন।
সমাধান। h = 2, T₁ = 2/2 · [3 + 3] = 6। ইন্টিগ্র্যান্ড ধ্রুবক, তাই একক ট্র্যাপিজয়েড আয়তক্ষেত্রের সাথে ঠিক মিলে যায়। সঠিক ইন্টিগ্রালও 6, তাই ফলাফল সঠিক।

সমস্যা 2. n = 1 সহ ট্র্যাপিজয়েডাল নিয়ম ব্যবহার করে [0, 1]-তে f(x) = x-এর চাপের দৈর্ঘ্য নির্ণয় করুন।
সমাধান। f'(x) = 1, তাই g(x) = √2, একটি ধ্রুবক। T₁ = 1 · √2 ≈ 1.414214, যা সঠিক কারণ ইন্টিগ্র্যান্ড ধ্রুবক। জ্যামিতিকভাবে এটি একটি একক বর্গের কর্ণ।

সমস্যা 3. n = 2 সহ [0, 1]-তে f(x) = x²-এর চাপের দৈর্ঘ্য অনুমান করুন।
সমাধান। h = 0.5; g(0) = 1, g(0.5) = 1.414214, g(1) = 2.236068। T₂ = 0.5 · [1/2 + 1.414214 + 2.236068/2] = 0.5 · 3.032249 = 1.516124। সঠিক মান 1.478943, তাই ত্রুটি 0.037181।

সমস্যা 4. n = 4 সহ [0, π]-তে f(x) = sin x-এর চাপের দৈর্ঘ্য অনুমান করুন।
সমাধান। g(x) = √(1 + cos²x), h = π/4। 0, π/4, π/2, 3π/4, π-তে মূল্যায়ন করে ওজন প্রয়োগ করলে T₄ ≈ 3.819944 পাওয়া যায়। n = 8-এ মান 3.820198, তাই n = 4 ফলাফল ইতিমধ্যে অভিসারী মান থেকে 0.000254-এর মধ্যে।

সমস্যা 5. a = 2, b = 5, n = 6-এর জন্য ধাপের আকার গণনা করুন।
সমাধান। h = (5 - 2)/6 = 3/6 = 0.5। নমুনা বিন্দুগুলো হলো 2, 2.5, 3, 3.5, 4, 4.5, 5: ছয়টি প্যানেলের জন্য সাতটি বিন্দু।

সমস্যা 6. n = 10-এর জন্য কয়টি নমুনা বিন্দু আছে, এবং তাদের ওজন কী?
সমাধান। n + 1 = 11টি নমুনা বিন্দু আছে। প্রান্তবিন্দু x₀ ও x₁₀ 1/2 ওজন বহন করে; নয়টি ভেতরের বিন্দু x₁ থেকে x₉ 1 ওজন বহন করে।

সমস্যা 7. [0, 1]-তে f(x) = x²-এর জন্য T₂ ও T₄ তুলনা করুন এবং পরিবর্তন বলুন।
সমাধান। T₂ = 1.516124 এবং T₄ = 1.488264, 0.027860 পরিবর্তন। সঠিক মান 1.478943, তাই n দ্বিগুণ করলে ত্রুটি 0.037181 থেকে 0.009321-এ নেমেছে, প্রায় চার গুণ, যেমন দ্বিতীয়-ক্রমের অভিসারণ ভবিষ্যদ্বাণী করে।

সমস্যা 8. T₈ = 1.481272, যেখানে [0, 1]-তে f(x) = x²। সঠিক 1.478943-এর বিপরীতে পরম ও আপেক্ষিক ত্রুটি জানান।
সমাধান। পরম ত্রুটি = |1.481272 - 1.478943| = 0.002330। আপেক্ষিক ত্রুটি = 0.002330 / 1.478943 ≈ 0.001576, প্রায় 0.16%। অনুমান ও ত্রুটি দুটোই এখানে তিন দশমিক স্থানে উল্লেখ করা উচিত; 1.481272-এর অতিরিক্ত অঙ্কগুলো এখনো বিশ্বাসযোগ্য নয়।