לדעתי בשורה האחרונה
בשורה האחרונה הפחתת המסה צריכה להיות:
zz -a*t+b*t^2/2 zz ולא כמו שרשמת (כי זו תוצאה של אינטגרציה על a-bt), כלומר חסרה חלוקה ב - 2 של b.
מעבר לזה, הדרך שלך נראית לי נכונה אם כי לא מתקבלת התשובה שאני קיבלתי (לעומת זה הפתרון שלך יותר דומה למשוואת הרקטה המקורית).
אני עשיתי על שני צידי המשוואה אינטגרציה לפי זמן, והצבתי את הנתון של dm/dt (ה - a-bt), עקרונית זה צריך להיות נכון אבל לא ברור לי למה התשובות שונות.