My OperaKurva
Saturday, November 15, 2008 3:32:26 AM
> with (student):
> p1:=[a-h,sin(a-h)];
> p2:=[a,sin(a)];
> p3:=[a+h,sin(a+h)];
> m1:=slope(p1,p2);
> m2:=slope(p2,p3);
> midp1:=midpoint(p1,p2);
> midp2:=midpoint(p2,p3);
> line1:=y=-1/m1*(x-midp1[1])+midp1[2];
> line2:=y=-1/m2*(x-midp2[1])+midp2[2];
> c:=solve({line1,line2},{x,y});
> assign(c);
x,y;
Error, (in assign) invalid arguments
> xCent:=limit(x,h=0);
> yCent:=limit(y,h=0);
> rad:=distance(p2,[xCent,yCent]);
> radi:=simplify(rad,trig);
> radiu:=factor(radi);
> radius:=unapply(radiu,a);
> xCentA:=simplify(subs(a=3*Pi/4,xCent));
> yCentA:=simplify(subs(a=3*Pi/4,yCent));
> oscCircle:=[xCentA+radius(3*Pi/4)*cos(t),yCentA+radius(3*Pi/4)*sin(t),t=0..2*Pi];
> plt1:=plot(sin,-Pi..2*Pi):
> plt2:=plot({oscCircle,[[xCentA,yCentA],[3*Pi/4,sin(3*Pi/4)]]}):
> plots[display]({plt1,plt2},scaling=constrained);
> curvature:=1/radius:
> plot({sin,curvature},0..2*Pi);
>
> p1:=[a-h,sin(a-h)];
> p2:=[a,sin(a)];
> p3:=[a+h,sin(a+h)];
> m1:=slope(p1,p2);
> m2:=slope(p2,p3);
> midp1:=midpoint(p1,p2);
> midp2:=midpoint(p2,p3);
> line1:=y=-1/m1*(x-midp1[1])+midp1[2];
> line2:=y=-1/m2*(x-midp2[1])+midp2[2];
> c:=solve({line1,line2},{x,y});
> assign(c);
x,y;
Error, (in assign) invalid arguments
> xCent:=limit(x,h=0);
> yCent:=limit(y,h=0);
> rad:=distance(p2,[xCent,yCent]);
> radi:=simplify(rad,trig);
> radiu:=factor(radi);
> radius:=unapply(radiu,a);
> xCentA:=simplify(subs(a=3*Pi/4,xCent));
> yCentA:=simplify(subs(a=3*Pi/4,yCent));
> oscCircle:=[xCentA+radius(3*Pi/4)*cos(t),yCentA+radius(3*Pi/4)*sin(t),t=0..2*Pi];
> plt1:=plot(sin,-Pi..2*Pi):
> plt2:=plot({oscCircle,[[xCentA,yCentA],[3*Pi/4,sin(3*Pi/4)]]}):
> plots[display]({plt1,plt2},scaling=constrained);
> curvature:=1/radius:
> plot({sin,curvature},0..2*Pi);
>
