<Text-field style="Heading 1" size="14" layout="Heading 1"><Font size="14">pag. 46, #17.</Font><Font size="18"> </Font><Font size="14">solving the initial value problem for a 1st order separable equation</Font></Text-field> you cannot use dsolve, you should integrate. 1-y; LCYiIiJGI0kieUc2IiEiIg== int(1/%,y); LCQtSSNsbkc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYjLCYiIiJGK0kieUdGKCEiIkYt As y(0)=3, 1-y is negative: we have to change the sign under the argument of the logaritm. -ln(y-1); LCQtSSNsbkc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYjLCZJInlHRigiIiIhIiJGLEYt %=int(x^3,x)+C; LywkLUkjbG5HNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHNiI2IywmSSJ5R0YpIiIiISIiRi1GLiwmKiRJInhHRikiIiUjRi1GMkkiQ0dGKUYt sol:=y(x)=solve(%,y); Ly1JInlHNiI2I0kieEdGJSwmLUkkZXhwRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliR0YlNiMsJiokRiciIiUjISIiRjFJIkNHRiVGMyIiIkY1RjU= 3=subs(x=0,rhs(sol)); LyIiJCwmLUkkZXhwRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiNiMsJEkiQ0dGKiEiIiIiIkYvRi8= C=solve(%,C); L0kiQ0c2IiwkLUkjbG5HNiQlKnByb3RlY3RlZEdJKF9zeXNsaWJHRiQ2IyIiIyEiIg== sol:=subs(%,sol); Ly1JInlHNiI2I0kieEdGJSwmLUkkZXhwRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliR0YlNiMsJiokRiciIiUjISIiRjEtSSNsbkdGKzYjIiIjIiIiRjhGOEY4 simplify(%); Ly1JInlHNiI2I0kieEdGJSwmLUkkZXhwRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliR0YlNiMsJCokRiciIiUjISIiRjEiIiMiIiJGNQ==
<Text-field style="Heading 1" size="14" layout="Heading 1"><Font size="14">pag. 54, #14.</Font><Font size="18"> </Font><Font size="14">obtaining the general solution of a 1st order linear equation</Font></Text-field> restart: The equation is not in standard form. It is imperative that the linear equation be in standard form before we attempt to determine the coefficient p of y. If the coefficient of y' is not one, then p and the resulting integrating factor will be incorrect; and the "solution" will not check. It is also necessary that the term involving y be on the same side as y' and that all x's be on the right hand side. Here is the original differential equation. odiffeq:=x*diff(y(x),x)+3*y(x)+2*x^2=x^3+4*x; LywoKiZJInhHNiIiIiItSSVkaWZmRyUqcHJvdGVjdGVkRzYkLUkieUdGJjYjRiVGJUYnRidGLCIiJCokRiUiIiNGMSwmKiRGJUYvRidGJSIiJQ== To rewrite it in standard form, we subtract 2*x^2 from both sides, then we divide by x. (In this simple case, it may be easier to retype the equation by hand.) ndiffeq:=diff(y(x),x)+(3*y(x))/x=x^2+4-2*x; LywmLUklZGlmZkclKnByb3RlY3RlZEc2JC1JInlHNiI2I0kieEdGKkYsIiIiKiZGKEYtRiwhIiIiIiQsKCokRiwiIiNGLSIiJUYtRiwhIiM= The expression p is the coefficient of y, including sign, when the equation is in standard form. p:=3/x; LCQqJEkieEc2IiEiIiIiJA== int(p,x); LCQtSSNsbkc2JCUqcHJvdGVjdGVkR0koX3N5c2xpYkc2IjYjSSJ4R0YoIiIk The integrating factor is the exponential of the antiderivative of p. mu:=exp(%); KiRJInhHNiIiIiQ= After the standard form equation is multiplied by the integrating factor, the left hand side is a perfect derivative. Hence, when we integrate, the left hand side will always be mu*y(x). To find y(x), we need only divide by mu. Note that Maple does not automatically insert a constant of integration, so we must do it manually. mu*ndiffeq; LyomSSJ4RzYiIiIkLCYtSSVkaWZmRyUqcHJvdGVjdGVkRzYkLUkieUdGJTYjRiRGJCIiIiomRixGL0YkISIiRiZGLyomRiRGJiwoKiRGJCIiI0YvIiIlRi9GJCEiI0Yv expand(%); LywmKiZJInhHNiIiIiQtSSVkaWZmRyUqcHJvdGVjdGVkRzYkLUkieUdGJjYjRiVGJSIiIkYvKiZGJSIiI0YsRi9GJywoKiRGJSIiJkYvKiRGJUYnIiIlKiRGJUY2ISIj mu*y(x)=int(rhs(%),x)+c; LyomSSJ4RzYiIiIkLUkieUdGJTYjRiQiIiIsKiokRiQiIicjRipGLSokRiQiIiVGKiokRiQiIiYjISIjRjJJImNHRiVGKg== solve(%,y(x)): sol:=y(x)=expand(%); Notice that x is different from 0 as we divided from it in one of the previous steps. Ly1JInlHNiI2I0kieEdGJSwqKiRGJyIiJCMiIiIiIidGJ0YsKiRGJyIiIyMhIiMiIiYqJkYnISIkSSJjR0YlRixGLA== Checking (in the original differential equation): subs(sol,odiffeq): simplify(%); LyomSSJ4RzYiIiIiLCYqJEYkIiIjRiYiIiVGJkYmLCYqJEYkIiIkRiZGJEYq
<Text-field style="Heading 1" size="14" layout="Heading 1"><Font size="14">pag. 54, #31. discontinuous coefficient</Font></Text-field> There are two approaches to solve this kind of situations. The 1st one is to consider two different differential equations. We use the 2nd one as dsolve command accepts the piecewise entry, producing a piecewise answer. restart: P:=piecewise(0<=x and x<=2,1,x>2,3); LUkqcGllY2V3aXNlRyUqcHJvdGVjdGVkRzYmMzEiIiFJInhHNiIxRikiIiMiIiIyRixGKSIiJA== diffeq:=diff(y(x),x)+P*y(x)=x; LywmLUklZGlmZkclKnByb3RlY3RlZEc2JC1JInlHNiI2I0kieEdGKkYsIiIiKiYtSSpwaWVjZXdpc2VHRiY2JjMxIiIhRiwxRiwiIiNGLTJGNkYsIiIkRi1GKEYtRi1GLA== inits:=y(0)=1; Ly1JInlHNiI2IyIiISIiIg== sol:=dsolve({diffeq,inits},y(x)); Notice that the program is computing the backward solution as well, just by putting P(x)=0 for x<0. The backward solution is meaningless fopr us, as our differential problem lives in the halfline [0,+infty). Ly1JInlHNiI2I0kieEdGJS1JKnBpZWNld2lzZUclKnByb3RlY3RlZEc2KDJGJyIiISwmIiIiRi8qJEYnIiIjI0YvRjEyRidGMSwoLUkkZXhwRzYkRipJKF9zeXNsaWJHRiU2IywkRichIiJGMSooRidGLy1GNkYmRi9GNUYvRi8qJkY9Ri9GNUYvRjsxRjFGJywqLUY2NiMsJkYnISIkIiIlRi9GMSomLUY2NiNGMUYvRkFGLyNGRSIiKiomLUY2NiMsJkYnIiIkISIlRi9GL0ZBRi8jRjtGSiooRidGL0ZMRi9GQUYvI0YvRk8= plot(rhs(sol),x=0..5); NigtJSdDVVJWRVNHNiQ3Xm83JCQiIiEhIiIkIiM1ISIiNyQkIjJrbW07YSlHXGEhIz0kIjEkKXA1d1BCJVsqISM7NyQkIjJMTEwkM3gmKSozIiEjPCQiMW5ca2E3c0MhKiEjOzckJCIuRDFSKCpSYyIhIzgkIjA0KlJjRE1vJykhIzo3JCQiMm5tO0gyUCJRPyEjPCQiMVImZk4yQi9OKSEjOzckJCIyTExMZVJ3WDUkISM8JCIxaUopKm9SIW93KCEjOzckJCIxTUwkM3glM3lUISM7JCIxWk4jKjQ5JnlNKCEjOzckJCIyMSt2JGZ5RzdaISM8JCIxdE8/eUYpcD4oISM7NyQkIjFubSJ6JTRcWV8hIzskIjEpPTk0Y108MyghIzs3JCQiLXYkZmw8dSYhIzckIjBIZnQnUTgwcSEjOjckJCIxTUxlUi0vUGkhIzskIjEmKSoqZkI6OmNwISM7NyQkIjFubSJ6V28pXG4hIzskIjFYWDdBMzhMcCEjOzckJCItRGNtcGlzISM3JCIxVjZbTXMqbyRwISM7NyQkIjFubSJIS2tJeighIzskIjFtWXdtYWFucCEjOzckJCIxTUxlKik+VkIkKSEjOyQiMS1nSE1tK0NxISM7NyQkIjFubVRnKCk0XykpISM7JCIxTmxWSitsL3IhIzs3JCQiLURKYnchUSohIzckIjFzI0h0dGokM3MhIzs3JCQiMm9tbVRJT28vIiEjOyQiMSMzKyE9XUgqWyghIzs3JCQiMk5MJDNfPmpVNiEjOyQiMWBOTERyKWUheSEjOzckJCIyLCsrRDt2L0QiISM7JCIxOSVmQVtEQEIpISM7NyQkIix2PWgoZTghIzUkIjEkPU1oVyE+RigpISM7NyQkIix2JFs2ajkhIzUkIjF0I0dKIilHOUUqISM7NyQkIjJNTGUqW3ooeWIiISM7JCIxJCkqNEJ4Ry96KiEjOzckJCIyb21tVFhnMG4iISM7JCIyJillZW9qVm8vIiEjOzckJCIyb21tbUo8Z3ciISM7JCIxbThTa0YvMzYhIzo3JCQiLkQxTWNxKD0hIzckIjJNKDRDQWk4JD0iISM7NyQkIjJMJGU5OzA/RT4hIzskIjInSEtoKCk+Zzw3ISM7NyQkIjJtbW07cFdgKD4hIzskIjIibzchSG5yRkQiISM7NyQkIjJPJGU5VE9FSD8hIzskIjJvQl9eUUAuQSIhIzs3JCQiMi4rRDFmIz0kMyMhIzskIjJvYyRlTHJZUzYhIzs3JCQiLUQiPUVYOCMhIzYkIjIydjFEQFQheTUhIzs3JCQiLnY9eHBlPSMhIzckIjIqelJfUCsocC0iISM7NyQkIjJQJGVSQTlXUkEhIzskIjFiJ0hQIypwLiUpKiEjOzckJCIxbm0iSDI4SUgjISM6JCIxWSdRdSg9SCwmKiEjOzckJCIyb20ienBTUyJSIyEjOyQiMVckZiQ+clFxISohIzs3JCQiL3YkNDFvV1cjISM4JCIxUShcTiJcK0EqKSEjOzckJCIyTEwkM18/YChcIyEjOyQiMWlUUCpvIVxAKSkhIzs3JCQiMSQzLSlRODRERCEjOiQiMFlRKD1seiZ5KSEjOjckJCIxTDNfRDFsX0QhIzokIjFfJUgjM1dBZygpISM7NyQkIjJNZVJBIio0LWUjISM7JCIwMGcpWzYoUnUpISM6NyQkIjJPTGUqKT5weGcjISM7JCIxXEI5b3NITygpISM7NyQkIjFubTt6K3ZiRSEjOiQiMTk+dThAVFQoKSEjOzckJCIudiRmNHQuRiEjNyQiMWQwOz5GR24oKSEjOzckJCIybjt6V2leYnYjISM7JCIxS1kxKlJQYCIpKSEjOzckJCIyUEwkZSpHc3QhRyEjOyQiMSV6TyV6T00iKSkpISM7NyQkIitEUlc5SCEiKiQiMXdcIzRLPVIxKiEjOzckJCItREpFPj5JISM2JCIxOT0uI1J5KilHKiEjOzckJCIyLitEMVJVMDckISM7JCIxWycqUl8rcFEmKiEjOzckJCItdj1TMkxLISM2JCIxNDkqKkc6dVUpKiEjOzckJCIyam1tO3ApPU1MISM7JCIyTDZGNE1cTCwiISM7NyQkIjItKyt2PV1AVyQhIzskIjI7QSNleD9zWDUhIzs3JCQiMUwkZSpbJHoqUk4hIzokIjI6Oio0QyRHZjIiISM7NyQkIjIxKytEWUtwayQhIzskIjIudl9jS1gnNDYhIzs3JCQiMmxtIkgycWNaUCEjOyQiMiY+S3BneiY9OSIhIzs3JCQiMi8rREo1ZkYmUSEjOyQiMnAkeV4nKSopKmU8IiEjOzckJCIybG1tVGcuYyZSISM7JCIyPDpINGJbJTQ3ISM7NyQkIi5EY0VzSzElISM3JCIyQm8nKXk1elpDIiEjOzckJCIxTUxMJCkqcHA7JSEjOiQiMk8hcFJQSSYqeTchIzs3JCQiMU1MM3hlLHRVISM6JCIyI3BuWmsqM1NKIiEjOzckJCIxbjtIZE89eVYhIzokIjJ0JXl4Wk0mKVs4ISM7NyQkIjIwKytdIz4jW1olISM7JCIyJkd3aVtGIzRRIiEjOzckJCIxbm1UJkchZSZlJSEjOiQiMVdZdTY7czw5ISM6NyQkIjFNTEwkKVFrJW8lISM6JCIyNyhcOlxUbV05ISM7NyQkIi5EMU1tLXolISM3JCIyUChwbzEqNGVbIiEjOzckJCIyMSt2JDQwTyIqWyEjOyQiMlA1MCFbWVk+OiEjOzckJCIjXSEiIiQiMV8qZXghUWtiOiEjOi0lJkNPTE9SRzYmJSRSR0JHJCIjNSEiIiQiIiEhIiIkIiIhISIiLSUlVklFV0c2JDskIiIhISIiJCIjXSEiIiUoREVGQVVMVEctJStBWEVTTEFCRUxTRzYnSSJ4RzYiUSE2Ii0lJUZPTlRHNiUlIUclIUciIzUlK0hPUklaT05UQUxHJStIT1JJWk9OVEFMRy0lKkFYRVNTVFlMRUc2IyUnTk9STUFMRy0lKFNDQUxJTkdHNiMlLlVOQ09OU1RSQUlORURHLSUlUk9PVEc2Jy0lKUJPVU5EU19YRzYjJCIkNSQhIiItJSlCT1VORFNfWUc2IyQiIyEqISIiLSUtQk9VTkRTX1dJRFRIRzYjJCIlNU8hIiItJS5CT1VORFNfSEVJR0hURzYjJCIlXU0hIiItJSlDSElMRFJFTkc2Ig==