Advanced Search

Ronan

1506 Reputation

16 Badges

14 years, 224 days
East Grinstead, United Kingdom
Contact Ronan

MaplePrimes Activity


These are answers submitted by Ronan

How did you install...

Ronan 1506
March 16 2022
2 2

 If you are using windows

Always with this programs like Maple I install them using the hidden administrator account in windows. i learn't this the hard way using the Solidworks CAD systtem. 

When I first used Maple in 2012 I forgot to do the above and had all sorts of problems until I reinstalled as above. 

Might help.

an equation or implicitplot or intersect...

Ronan 1506
February 22 2022
4 0 
f-0.001-g #gives the intersection curve.

there ae=re two basic ways to show it implicitplot or in 3D intersectplot.


 

NULL

> 

with(plots)

[animate, animate3d, animatecurve, arrow, changecoords, complexplot, complexplot3d, conformal, conformal3d, contourplot, contourplot3d, coordplot, coordplot3d, densityplot, display, dualaxisplot, fieldplot, fieldplot3d, gradplot, gradplot3d, implicitplot, implicitplot3d, inequal, interactive, interactiveparams, intersectplot, listcontplot, listcontplot3d, listdensityplot, listplot, listplot3d, loglogplot, logplot, matrixplot, multiple, odeplot, pareto, plotcompare, pointplot, pointplot3d, polarplot, polygonplot, polygonplot3d, polyhedra_supported, polyhedraplot, rootlocus, semilogplot, setcolors, setoptions, setoptions3d, shadebetween, spacecurve, sparsematrixplot, surfdata, textplot, textplot3d, tubeplot]

 (1)
> 

restart
> 

NULL
> 

f := -.497068470024476*x^3*y^4-6.70447066551140*x^5*y^3-4.75254543704658*x^5*y^2+29.2165392673953*x^5*y^4+19.7458656447309*x^5*y^6-34.1911496490261*x^5*y^5+1.01958235276613*x^5*y^8-6.20249607007329*x^5*y^7+.122625097102966*x^3*y^5-.557878506354632*x^4*y+1.10377462803280*10^(-9)*y^9-3.01639944258456*10^(-8)*y^8+2.70433066580774*10^(-7)*y^7-58.9547371078715*x^6*y^4+70.6881059337016*x^6*y^5-41.4723627878841*x^6*y^6+13.2184205014647*x^6*y^7+1.88972684858187*x^5*y-2.20596265158541*x^6*y^8+0.415520691947419e-1*x*y-0.264590704119876e-1*x^9*y^9-0.118773599704818e-5*y^6+0.290568461642441e-5*y^5-0.409894118008625e-5*y^4+0.324884379509192e-5*y^3+11.9305450198352*x^6*y^3+10.8462050916126*x^6*y^2-4.10654876148122*x^6*y+2.54724280276503*x^7*y^8-14.9818956625447*x^7*y^7+46.0951349392065*x^7*y^6-0.109255815129257e-2*x^9+0.788567804811102e-2*x^3*y^8-0.131109942751000e-5*y^2+0.249200320311903e-4*x^2+0.514872391346379e-2*x^8-0.100662629994939e-1*x^7+0.105239075112828e-1*x^6-0.630639949679202e-2*x^5-0.375713420739403e-3*x^3+0.214317754102793e-2*x^4+.876088098407044*x^9*y-2.70457990868458*x^9*y^2+.273231436604042*x^9*y^3+6.75578062548269*x^9*y^4-9.57907483795682*x^9*y^5+6.09218732418035*x^9*y^6+.722761182584203*x^3*y^3-.466894994882040*x^3*y^2+.101139100881401*x^3*y-.223464159932998*x^4*y^8+1.35561916171540*x^4*y^7-4.28474908826935*x^4*y^6+7.29078271720222*x^4*y^5-5.90416839043862*x^4*y^4+.791654444114778*x^4*y^3+1.52828070029702*x^4*y^2+.364411288643973*x^9*y^8-76.7100686748060*x^7*y^5-3.35171664262170*x^8*y+9.91733485134206*x^8*y^2-0.148435770331285e-1*x^2*y-.276512584364289*x^2*y^3+.104481013510196*x^2*y^2-.279034248189274*x^2*y^5+.370033041604670*x^2*y^4+0.399016415926658e-2*x^2*y^8-0.307519757322687e-1*x^2*y^7+.122644475629609*x^2*y^6-0.261755569520174e-2*x*y^3-0.956445131268424e-2*x*y^2-0.525240272678308e-2*x*y^6+0.294149187341931e-1*x*y^5-0.364419367372761e-1*x*y^4-0.690607476536226e-1*x^5*y^9+0.151457171844930e-1*x^4*y^9+0.671188193729842e-3*x*y^8-0.156037383968479e-2*x*y^7+.108086167003366*x^8*y^9-.178621376287247*x^7*y^9-0.705741073722833e-3*x^3*y^9+.151841215581478*x^6*y^9+2.06864668194457*10^(-7)*x+2.06863293440197*10^(-7)*y-0.199666299235104e-3*x^2*y^9-0.673706850700663e-4*x*y^9+2.00727312663099*x^8*y^3-32.1473826670758*x^8*y^4+42.6266489705831*x^8*y^5-26.3247561285797*x^8*y^6+8.73691621490782*x^8*y^7-2.06514882458964*x^9*y^7+0.363726776841009e-1*x^3*y^6-0.317849887375473e-1*x^3*y^7-1.51380636269074*x^8*y^8+5.14796137189361*x^7*y-14.4626792261666*x^7*y^2-8.76022418809059*x^7*y^3+61.1998773291695*x^7*y^4-1.79614673481709*10^(-9)-0.1e-2

0.2068646682e-6*x+0.2068632934e-6*y-0.1000001796e-2-.557878506354632*x^4*y+1.52828070029702*x^4*y^2+.791654444114778*x^4*y^3-5.90416839043862*x^4*y^4+7.29078271720222*x^4*y^5-4.28474908826935*x^4*y^6+1.35561916171540*x^4*y^7-.223464159932998*x^4*y^8+.101139100881401*x^3*y-.466894994882040*x^3*y^2+.722761182584203*x^3*y^3-.497068470024476*x^3*y^4+.122625097102966*x^3*y^5+0.363726776841009e-1*x^3*y^6-0.317849887375473e-1*x^3*y^7+0.788567804811102e-2*x^3*y^8+.876088098407044*x^9*y-2.70457990868458*x^9*y^2+.273231436604042*x^9*y^3+6.75578062548269*x^9*y^4-9.57907483795682*x^9*y^5+6.09218732418035*x^9*y^6-2.06514882458964*x^9*y^7+.364411288643973*x^9*y^8-3.35171664262170*x^8*y+9.91733485134206*x^8*y^2+2.00727312663099*x^8*y^3-32.1473826670758*x^8*y^4+42.6266489705831*x^8*y^5-26.3247561285797*x^8*y^6+8.73691621490782*x^8*y^7-1.51380636269074*x^8*y^8+5.14796137189361*x^7*y-14.4626792261666*x^7*y^2-8.76022418809059*x^7*y^3+61.1998773291695*x^7*y^4-76.7100686748060*x^7*y^5+46.0951349392065*x^7*y^6-14.9818956625447*x^7*y^7+2.54724280276503*x^7*y^8-4.10654876148122*x^6*y+11.9305450198352*x^6*y^3+10.8462050916126*x^6*y^2-58.9547371078715*x^6*y^4+70.6881059337016*x^6*y^5-41.4723627878841*x^6*y^6+13.2184205014647*x^6*y^7-2.20596265158541*x^6*y^8+0.415520691947419e-1*x*y+1.88972684858187*x^5*y-0.264590704119876e-1*x^9*y^9-4.75254543704658*x^5*y^2+29.2165392673953*x^5*y^4-6.70447066551140*x^5*y^3-34.1911496490261*x^5*y^5-6.20249607007329*x^5*y^7+19.7458656447309*x^5*y^6+1.01958235276613*x^5*y^8+.104481013510196*x^2*y^2-0.148435770331285e-1*x^2*y+.370033041604670*x^2*y^4-.276512584364289*x^2*y^3+.122644475629609*x^2*y^6-.279034248189274*x^2*y^5+0.399016415926658e-2*x^2*y^8-0.307519757322687e-1*x^2*y^7-0.261755569520174e-2*x*y^3-0.956445131268424e-2*x*y^2-0.525240272678308e-2*x*y^6+0.294149187341931e-1*x*y^5-0.364419367372761e-1*x*y^4+0.151457171844930e-1*x^4*y^9+0.671188193729842e-3*x*y^8-0.156037383968479e-2*x*y^7+.151841215581478*x^6*y^9-0.690607476536226e-1*x^5*y^9-0.199666299235104e-3*x^2*y^9-0.673706850700663e-4*x*y^9+.108086167003366*x^8*y^9-.178621376287247*x^7*y^9-0.705741073722833e-3*x^3*y^9-0.109255815129257e-2*x^9-0.409894118008625e-5*y^4+0.324884379509192e-5*y^3+0.290568461642441e-5*y^5+0.2704330666e-6*y^7-0.118773599704818e-5*y^6-0.3016399443e-7*y^8-0.630639949679202e-2*x^5+0.214317754102793e-2*x^4+0.105239075112828e-1*x^6-0.131109942751000e-5*y^2-0.100662629994939e-1*x^7-0.375713420739403e-3*x^3+0.514872391346379e-2*x^8+0.249200320311903e-4*x^2+0.1103774628e-8*y^9

 (2)
> 

g := 0.486130607155690e-1*x^8*y^9-0.793915710864734e-1*x^7*y^9+.166316351180282*x^9*y^2+0.100163106986489e-1*x^9*y-0.725705040036393e-1*x^5*y+9.99632884281232*10^(-10)-14.7626936165794*x^6*y^6-0.257557973826478e-3*x^3*y^9-5.72759908720050*10^(-8)*y^2+0.662263463130918e-1*x^6*y^9-0.292171279505531e-1*x^5*y^9+0.608295954591597e-2*x^4*y^9-0.128182657077761e-3*x*y^8+0.942633285021959e-3*x*y^7-0.350081319685295e-2*x*y^6+0.688776508883605e-2*x*y^5-0.759510525544657e-2*x*y^4+0.779592933988801e-2*x*y^3-0.157071356900974e-2*x*y^2+0.120360948207261e-2*x^2*y^8-0.763275503500010e-2*x^2*y^7+0.261578297444302e-1*x^2*y^6-0.525360371999116e-1*x^2*y^5+0.625949969717321e-1*x^2*y^4-0.425238644768537e-1*x^2*y^3-3.86106645787476*10^(-10)*x-3.85902759075275*10^(-10)*y-0.790340436745485e-4*x^2*y^9+0.686404280878835e-5*x*y^9+1.61849630014731*10^(-9)*y^9-2.09189804289573*10^(-8)*y^8+1.12047595681372*10^(-7)*y^7+0.206051274192785e-2*x^6-0.653096669907143e-4*x^3-14.1011477616247*x^8*y^4+16.1206664094461*x^8*y^5+24.3848818418075*x^6*y^5-0.119825252833206e-1*x^9*y^9-0.194739348540728e-1*x*y-22.5326440659846*x^6*y^4-10.1171356900718*x^8*y^6+3.54607571091852*x^8*y^7+24.8991084141403*x^7*y^4+0.108822319123503e-2*x^8-0.204991873524550e-2*x^7-0.130932246100715e-1*x^3*y^2+0.366339579352754e-2*x^3*y^8-0.205724334437174e-1*x^3*y^7+1.18133952115280*x^5*y^2-5.37560745909542*x^5*y^3-11.4126688531743*x^5*y^5+10.8065136845075*x^5*y^4+17.0670410168283*x^7*y^6+6.78593938914942*x^5*y^6-5.90525014814763*x^7*y^7+1.07208303147736*x^7*y^8+.109044949802749*x^6*y-2.20370941827329*x^6*y^2+0.111947605945080e-1*x^3*y+10.7889789185019*x^6*y^3+.159056987714934*x^9*y^8+5.03727993859085*x^6*y^7-2.40269859450631*x^4*y^4+1.20635516320513*x^4*y^3-.265458418460343*x^4*y^2+0.150663088709287e-1*x^4*y-1.31283787245192*x^9*y^3-.903861947690635*x^6*y^8-0.239122454967206e-3*x^9-0.165983942264791e-2*x^3*y^3-3.21279186038477*10^(-7)*y^6+0.540484515816569e-1*x^3*y^4-0.842888013414850e-1*x^3*y^5+0.579687010994481e-1*x^3*y^6-27.6818094812177*x^7*y^5-2.28138893951844*x^5*y^7-0.192976732163787e-2*x^2*y+.403933247361425*x^5*y^8+0.147245143634235e-1*x^2*y^2+2.59098767602453*10^(-7)*y^3-.858862881780697*x^9*y^7+2.42075981205875*x^9*y^6-0.308703623001616e-2*x^8*y-.961875822502716*x^8*y^2+2.50971144938558*x^4*y^5-1.47470431176274*x^4*y^6+.489488747286588*x^4*y^7-0.854522394871459e-1*x^4*y^8+0.386516741572877e-3*x^4+0.431899528468855e-5*x^2+5.31031847990512*10^(-7)*y^5-0.118522142744809e-2*x^5-3.79079342069720*x^9*y^5+3.22173828058960*x^9*y^4-5.05713571671948*10^(-7)*y^4+6.10514710675470*x^8*y^3+2.08327870192822*x^7*y^2-0.610537804766400e-1*x^7*y-.650514845456873*x^8*y^8-11.3732791811746*x^7*y^3

-0.3861066458e-9*x-0.3859027591e-9*y+0.150663088709287e-1*x^4*y-.265458418460343*x^4*y^2+1.20635516320513*x^4*y^3-2.40269859450631*x^4*y^4+2.50971144938558*x^4*y^5-1.47470431176274*x^4*y^6+.489488747286588*x^4*y^7-0.854522394871459e-1*x^4*y^8+0.111947605945080e-1*x^3*y-0.130932246100715e-1*x^3*y^2-0.165983942264791e-2*x^3*y^3+0.540484515816569e-1*x^3*y^4-0.842888013414850e-1*x^3*y^5+0.579687010994481e-1*x^3*y^6-0.205724334437174e-1*x^3*y^7+0.366339579352754e-2*x^3*y^8+0.100163106986489e-1*x^9*y+.166316351180282*x^9*y^2-1.31283787245192*x^9*y^3+3.22173828058960*x^9*y^4-3.79079342069720*x^9*y^5+2.42075981205875*x^9*y^6-.858862881780697*x^9*y^7+.159056987714934*x^9*y^8-0.308703623001616e-2*x^8*y-.961875822502716*x^8*y^2+6.10514710675470*x^8*y^3-14.1011477616247*x^8*y^4+16.1206664094461*x^8*y^5-10.1171356900718*x^8*y^6+3.54607571091852*x^8*y^7-.650514845456873*x^8*y^8-0.610537804766400e-1*x^7*y+2.08327870192822*x^7*y^2-11.3732791811746*x^7*y^3+24.8991084141403*x^7*y^4-27.6818094812177*x^7*y^5+17.0670410168283*x^7*y^6-5.90525014814763*x^7*y^7+1.07208303147736*x^7*y^8+.109044949802749*x^6*y+10.7889789185019*x^6*y^3-2.20370941827329*x^6*y^2-22.5326440659846*x^6*y^4+24.3848818418075*x^6*y^5-14.7626936165794*x^6*y^6+5.03727993859085*x^6*y^7-.903861947690635*x^6*y^8-0.194739348540728e-1*x*y-0.725705040036393e-1*x^5*y-0.119825252833206e-1*x^9*y^9+1.18133952115280*x^5*y^2+10.8065136845075*x^5*y^4-5.37560745909542*x^5*y^3-11.4126688531743*x^5*y^5-2.28138893951844*x^5*y^7+6.78593938914942*x^5*y^6+.403933247361425*x^5*y^8+0.147245143634235e-1*x^2*y^2-0.192976732163787e-2*x^2*y+0.625949969717321e-1*x^2*y^4-0.425238644768537e-1*x^2*y^3+0.261578297444302e-1*x^2*y^6-0.525360371999116e-1*x^2*y^5+0.120360948207261e-2*x^2*y^8-0.763275503500010e-2*x^2*y^7+0.779592933988801e-2*x*y^3-0.157071356900974e-2*x*y^2-0.350081319685295e-2*x*y^6+0.688776508883605e-2*x*y^5-0.759510525544657e-2*x*y^4+0.608295954591597e-2*x^4*y^9-0.128182657077761e-3*x*y^8+0.942633285021959e-3*x*y^7+0.662263463130918e-1*x^6*y^9-0.292171279505531e-1*x^5*y^9-0.790340436745485e-4*x^2*y^9+0.686404280878835e-5*x*y^9+0.486130607155690e-1*x^8*y^9-0.793915710864734e-1*x^7*y^9-0.257557973826478e-3*x^3*y^9-0.239122454967206e-3*x^9-0.5057135717e-6*y^4+0.2590987676e-6*y^3+0.5310318480e-6*y^5+0.1120475957e-6*y^7-0.3212791860e-6*y^6-0.2091898043e-7*y^8-0.118522142744809e-2*x^5+0.386516741572877e-3*x^4+0.206051274192785e-2*x^6-0.5727599087e-7*y^2-0.204991873524550e-2*x^7-0.653096669907143e-4*x^3+0.108822319123503e-2*x^8+0.431899528468855e-5*x^2+0.1618496300e-8*y^9+0.9996328843e-9

 (3)
> 

P1 := plots:-display(plot3d(f-0.1e-2, x = 0 .. 1, y = 0 .. 2, color = red, style = surface), plot3d(g, x = 0 .. 1, y = 0 .. 2, color = green, style = wireframe))

 
> 

fg := simplify(f-0.1e-2-g)

(HFloat(0.8660717876999807)*y-HFloat(2.870896259999945)*y^2+HFloat(1.5860693089999662)*y^3+HFloat(3.5340423439999995)*y^4-HFloat(5.788281417)*y^5+HFloat(3.671427512)*y^6-HFloat(1.2062859430000001)*y^7+HFloat(0.20535430090000004)*y^8-0.1447654513e-1*y^9-HFloat(8.534356959999808e-4))*x^9+(-HFloat(3.3486296069999995)*y+HFloat(10.879210669999999)*y^2-HFloat(4.097873979999998)*y^3-HFloat(18.04623491)*y^4+HFloat(26.505982559999996)*y^5-HFloat(16.20762044)*y^6+HFloat(5.190840504)*y^7-HFloat(0.8632915174999999)*y^8+0.5947310628e-1*y^9+HFloat(0.004060500722000001))*x^8+(HFloat(5.209015151999995)*y-HFloat(16.54595792999999)*y^2+HFloat(2.613054992000008)*y^3+HFloat(36.30076891999998)*y^4-HFloat(49.02825919)*y^5+HFloat(29.028093920000003)*y^6-HFloat(9.076645512)*y^7+HFloat(1.475159772)*y^8-0.9922980521e-1*y^9-HFloat(0.008016344264999992))*x^7+(-HFloat(4.215593710999997)*y+HFloat(1.1415660999999866)*y^3+HFloat(13.049914509999999)*y^2-HFloat(36.42209304000001)*y^4+HFloat(46.303224090000015)*y^5-HFloat(26.70966917)*y^6+HFloat(8.181140561000001)*y^7-HFloat(1.3021007040000003)*y^8+0.8561486929e-1*y^9+HFloat(0.008463394767999994))*x^6+(HFloat(1.9622973529999983)*y-HFloat(5.9338849579999975)*y^2+HFloat(18.410025589999996)*y^4-HFloat(1.3288632069999955)*y^3-HFloat(22.778480799999997)*y^5-HFloat(3.92110713)*y^7+HFloat(12.95992625)*y^6+HFloat(0.6156491056)*y^8-0.3984361970e-1*y^9-HFloat(0.005121178069999996))*x^5+(-HFloat(0.5729448153000001)*y+HFloat(1.7937391179999997)*y^2-HFloat(0.4147007189000007)*y^3-HFloat(3.5014697949999984)*y^4+HFloat(4.781071268)*y^5-HFloat(2.8100447759999994)*y^6+HFloat(0.8661304146999997)*y^7-HFloat(0.1380119204)*y^8+0.9062757634e-2*y^9+HFloat(0.0017566607989999996))*x^4+(HFloat(0.08994434031000005)*y-HFloat(0.4538017703000001)*y^2+HFloat(0.7244210220000001)*y^3-HFloat(0.5511169216)*y^4+HFloat(0.2069138984)*y^5-HFloat(0.021596023419999985)*y^6-HFloat(0.011212555300000001)*y^7+HFloat(0.004222282254000001)*y^8-0.4481830999e-3*y^9-HFloat(3.1040375370000004e-4))*x^3+(HFloat(0.08975649913999999)*y^2-HFloat(0.01291380970999996)*y+HFloat(0.3074380445999999)*y^4-HFloat(0.23398871989999998)*y^3+HFloat(0.09648664585999998)*y^6-HFloat(0.22649821099999995)*y^5+HFloat(0.0027865546769999996)*y^8-HFloat(0.023119220699999995)*y^7-0.1206322555e-3*y^9+HFloat(2.060103673999994e-5))*x^2+(HFloat(0.06102600403999991)*y-HFloat(0.010413485039999991)*y^3-HFloat(0.007993737744000015)*y^2-HFloat(0.001751589530000003)*y^6+HFloat(0.02252715363999998)*y^5-HFloat(0.028846831479999953)*y^4+HFloat(7.993708507999997e-4)*y^8-HFloat(0.002503007124999997)*y^7-0.7423472788e-4*y^9+HFloat(2.072507747999997e-7))*x-HFloat(0.002000002796195685)-HFloat(3.5932276091581405e-6)*y^4+HFloat(2.989745027101643e-6)*y^3+HFloat(2.3746527682378423e-6)*y^5+HFloat(1.5838547089901542e-7)*y^7-HFloat(8.664568110100262e-7)*y^6-HFloat(9.245014000011583e-9)*y^8-HFloat(1.2538234244762572e-6)*y^2-0.514721672e-9*y^9+HFloat(2.0724919623361815e-7)*y

 (4)
> 

P2 := plots:-implicitplot(fg, x = 0 .. 1, y = 0 .. 2)

 

``

P3 := plots:-intersectplot(f-0.1e-2, g, x = 0 .. 1, y = 0 .. 2, colour = black, thickness = 5)

  

plots:-display(P1, P3)

  

NULL


 

Download 22-02-2022_contact.mw

Here is an way...

Ronan 1506
February 14 2022
1 0

This is a something I did a few years ago. I have crudely edited it to suit. 

restartNULL

with(plots)

with(ColorTools)

Clrs := GetPalette("Patchwork")

Color(Clrs[4])

NULL

``

"f(x):=x^(3)-1 ; frts:=[fsolve(f(x),x,complex)]; R1:=plot([Re,Im]~(frts), style=point, color=black, symbol=solidcircle, symbolsize=20, scaling=constrained) ; "

NULL

g := unapply(diff(f(x), x), x)

rts := frts

nops(rts)

Listx := Vector(); Listy := Vector(); C := Array(datatype = float[8])
 

 
xa := 2*min(`~`[Re](rts)); xb := 1.251*max(`~`[Re](rts)); ya := 1.25*min(`~`[Im](rts)); yb := 1.25*max(`~`[Im](rts)); incriment := 401

``

Newton := proc (XA := -2, XB := 2, YA := -2, YB := 2, INC := 400) local x1, xn, `&epsilon;`, j, x, y, itter, err; global rts; j := 1; for x from XA by (XB-XA)/INC to XB do for y from YA by (YB-YA)/INC to YB do Listx(j) := x; Listy(j) := y; `&epsilon;` := 1; x1 := x+0.1e-6+I*y; itter := 0; while 0.5e-2 < `&epsilon;` and itter < 25 do xn := evalhf(x1-f(x1)/g(x1)); err := xn-x1; `&epsilon;` := abs(err); x1 := xn; itter := itter+1 end do; C(j) := min(min[index](`~`[abs]([seq(rts[j]-x1, j = 1 .. nops(rts))]))/nops(rts), 1.0); j := j+1 end do end do end proc

````

``

````

NULLNewton(xa, xb, ya, yb, incriment)

p2 := plots:-pointplot(Listx, Listy, symbol = point, symbolsize = 1, color = COLOR(HUE, C), size = [incriment, incriment])

display(p2, R1)``

Download 14-2-22_Newton_complex.mw

Using algebra and equating coefficients...

Ronan 1506
February 13 2022
0 0

Long handed way of solving the problem
 

restart

See https://en.wikipedia.org/wiki/Parabola

 

with(plots)

NULL

``

NULLNULLParabola := proc (x, y) options operator, arrow; a*x^2+2*b*x*y+c*y^2+d*x+e*y+g end proc

proc (x, y) options operator, arrow; a*x^2+2*b*x*y+c*y^2+d*x+e*y+g end proc

 (1)

cond := a*c-b^2 = 0

a*c-b^2 = 0

 (2)

c := solve(cond, c)

b^2/a

 (3)

Parabola(x, y)

a*x^2+2*b*x*y+b^2*y^2/a+d*x+e*y+g

 (4)

``

``

A parabola can be written in terms of the distance from directrix=distance from focus

``

Par1 := sort(numer(expand((A*x+B*y+C)^2/(A^2+B^2)-(x-f[1])^2-(y-f[2])^2)), [x, y])

-B^2*x^2+2*A*B*x*y-A^2*y^2+2*A^2*f[1]*x+2*B^2*f[1]*x+2*A*C*x+2*A^2*f[2]*y+2*B^2*f[2]*y+2*B*C*y-A^2*f[1]^2-A^2*f[2]^2-B^2*f[1]^2-B^2*f[2]^2+C^2

 (5)

Equate the coefficients of Parabola and Par1

eq1 := coeff(Par1, x, 2) = coeff(Parabola(x, y), x, 2)

-B^2 = a

 (6)

eq2 := coeff(coeff(Par1, x, 1), y, 1) = coeff(coeff(Parabola(x, y), x, 1), y, 1)

2*A*B = 2*b

 (7)

NULL

eq3 := coeff(Par1, x, 1) = coeff(Parabola(x, y), x, 1)

2*A^2*f[1]+2*A*B*y+2*B^2*f[1]+2*A*C = 2*b*y+d

 (8)

eq4 := coeff(Par1, y, 1) = coeff(Parabola(x, y), y, 1)

2*A^2*f[2]+2*A*B*x+2*B^2*f[2]+2*B*C = 2*b*x+e

 (9)

eq5 := coeff(coeff(Par1, x, 0), y, 0) = coeff(coeff(Parabola(x, y), x, 0), y, 0)

-A^2*f[1]^2-A^2*f[2]^2-B^2*f[1]^2-B^2*f[2]^2+C^2 = g

 (10)

sol := solve({eq1, eq2, eq3, eq4, eq5}, [A, B, C, f[1], f[2]], explicit)

[[A = -b*(-a)^(1/2)/a, B = (-a)^(1/2), C = (1/4)*(4*a^2*g-a*d^2-a*e^2+4*b^2*g)*(-a)^(1/2)/(a*(a*e-b*d)), f[1] = (1/4)*(4*a^2*b*g-2*a^2*d*e+a*b*d^2-a*b*e^2+4*b^3*g)/((a^2+b^2)*(a*e-b*d)), f[2] = -(1/4)*a*(4*a^2*g-a*d^2+a*e^2+4*b^2*g-2*b*d*e)/(a^3*e-a^2*b*d+a*b^2*e-b^3*d)], [A = b*(-a)^(1/2)/a, B = -(-a)^(1/2), C = -(1/4)*(4*a^2*g-a*d^2-a*e^2+4*b^2*g)*(-a)^(1/2)/(a*(a*e-b*d)), f[1] = (1/4)*(4*a^2*b*g-2*a^2*d*e+a*b*d^2-a*b*e^2+4*b^3*g)/((a^2+b^2)*(a*e-b*d)), f[2] = -(1/4)*a*(4*a^2*g-a*d^2+a*e^2+4*b^2*g-2*b*d*e)/(a^3*e-a^2*b*d+a*b^2*e-b^3*d)]]

 (11)

assign(sol[1])

directrix := A*x+B*y+C

-b*(-a)^(1/2)*x/a+(-a)^(1/2)*y+(1/4)*(4*a^2*g-a*d^2-a*e^2+4*b^2*g)*(-a)^(1/2)/(a*(a*e-b*d))

 (12)

[f[1], f[2]]

[(1/4)*(4*a^2*b*g-2*a^2*d*e+a*b*d^2-a*b*e^2+4*b^3*g)/((a^2+b^2)*(a*e-b*d)), -(1/4)*a*(4*a^2*g-a*d^2+a*e^2+4*b^2*g-2*b*d*e)/(a^3*e-a^2*b*d+a*b^2*e-b^3*d)]

 (13)

axis := simplify(B*(x-f[1])-A*(y-f[2]))

-(2*a^3*x+2*a^2*b*y+2*a*b^2*x+2*b^3*y+a^2*d+a*b*e)/((-a)^(1/2)*(2*a^2+2*b^2))

 (14)

parameters := [a = 1, b = 2, d = 3, e = -1, g = 4]

[a = 1, b = 2, d = 3, e = -1, g = 4]

 (15)

display(implicitplot([eval(Parabola(x, y), parameters), eval(directrix, parameters), eval(axis, parameters)], x = -7 .. 1, y = -1 .. 5, colour = [red, green, blue]), pointplot([eval(f[1], parameters), eval(f[2], parameters)]))

  

NULL


 

Download 13-2-22_Parabola_axis_focus.mw

My variant on a solution...

Ronan 1506
February 11 2022
1 0

My variant on solving it. 

@vv  why did you define two parabolas? I see it gives more different solutions.

Edit:- corrected cond to b^2-a c=0


restart

with(plots)

CircleParm := proc (t) options operator, arrow; [(-t^2+1)/(t^2+1), 2*t/(t^2+1)] end proc

proc (t) options operator, arrow; [(-t^2+1)/(t^2+1), 2*t/(t^2+1)] end proc

 (1)

P1 := CircleParm(1/4)

[15/17, 8/17]

 (2)

P2 := CircleParm(5)

[-12/13, 5/13]

 (3)

P3 := CircleParm(-1/10)

[99/101, -20/101]

 (4)

P4 := CircleParm(-3/2)

[-5/13, -12/13]

 (5)

Pts := [P1, P2, P3, P4]

[[15/17, 8/17], [-12/13, 5/13], [99/101, -20/101], [-5/13, -12/13]]

 (6)

para := proc (x, y) options operator, arrow; a*x^2+c*y^2+2*b*x*y+e*x+f*y+1 end proc

proc (x, y) options operator, arrow; a*x^2+c*y^2+2*b*x*y+e*x+f*y+1 end proc

 (7)

cond := -a*c+b^2 = 0

-a*c+b^2 = 0

 (8)

``

NULL

NULL

sys1 := [seq(para(op(Pts[i])), i = 1 .. 4), cond]

[(225/289)*a+(64/289)*c+(240/289)*b+(15/17)*e+(8/17)*f+1, (144/169)*a+(25/169)*c-(120/169)*b-(12/13)*e+(5/13)*f+1, (9801/10201)*a+(400/10201)*c-(3960/10201)*b+(99/101)*e-(20/101)*f+1, (25/169)*a+(144/169)*c+(120/169)*b-(5/13)*e-(12/13)*f+1, -a*c+b^2 = 0]

 (9)

sol := solve(sys1, explicit)

{a = -400221/504721-(4940/504721)*3434^(1/2), b = 106975/504721-(5057/504721)*3434^(1/2), c = -760471/504721+(12090/504721)*3434^(1/2), e = -71500/504721+(3380/504721)*3434^(1/2), f = -107250/504721+(5070/504721)*3434^(1/2)}, {a = -400221/504721+(4940/504721)*3434^(1/2), b = 106975/504721+(5057/504721)*3434^(1/2), c = -760471/504721-(12090/504721)*3434^(1/2), e = -71500/504721-(3380/504721)*3434^(1/2), f = -107250/504721-(5070/504721)*3434^(1/2)}

 (10)

Parab1 := eval(para(x, y), sol[1])

(-400221/504721-(4940/504721)*3434^(1/2))*x^2+2*(106975/504721-(5057/504721)*3434^(1/2))*x*y+(-760471/504721+(12090/504721)*3434^(1/2))*y^2+(-71500/504721+(3380/504721)*3434^(1/2))*x+(-107250/504721+(5070/504721)*3434^(1/2))*y+1

 (11)

NULL

Parab2 := eval(para(x, y), sol[2])

(-400221/504721+(4940/504721)*3434^(1/2))*x^2+2*(106975/504721+(5057/504721)*3434^(1/2))*x*y+(-760471/504721-(12090/504721)*3434^(1/2))*y^2+(-71500/504721-(3380/504721)*3434^(1/2))*x+(-107250/504721-(5070/504721)*3434^(1/2))*y+1

 (12)

display(implicitplot([Parab1, Parab2, x^2+y^2-1], x = -2 .. 2, y = -4 .. 2, colour = [red, green, blue], scaling = constrained), pointplot(Pts, symbolsize = 16))

  

NULL

Download Parabolas_CoCyclic_Points.mw

Another approach. Shift and Truncate...

Ronan 1506
February 09 2022
1 0

Just in the spirit of another way to find a tangent without differentiation,

restart

with(plots)

conic := proc (x, y) options operator, arrow; 7*x^2-7*y^2-12*x+9*y+2.25-2*x*y end proc

proc (x, y) options operator, arrow; 7*x^2-7*y^2-12*x+9*y+2.25-2*y*x end proc

 (1)

x0, y0 := .504244923, .3781836925

.504244923, .3781836925

 (2)

display(implicitplot(conic(x, y), x = -1 .. 2, y = -2 .. 2), pointplot([x0, y0]))

  

x0, y0 := .504244923, .3781836925

.504244923, .3781836925

 (3)

conic(x0, y0)

0.22e-8

 (4)

Translate the conic at point [x0,y0] to the origin

 

conic1 := expand(conic(x+x0, y+y0))

7*x^2-5.696938463*x+0.22e-8-7*y^2+2.696938459*y-2*y*x

 (5)

NULL````````

implicitplot(conic1, x = -1 .. 1, y = -2 .. 2)

  

Truncate the conic to 1st order. this gives the tangent line to this conic at the origin

tanline1 := mtaylor(conic1, [x, y], 2)

0.22e-8-5.696938463*x+2.696938459*y

 (6)

implicitplot([conic1, tanline1], x = -1 .. 1, y = -1 .. 1, colour = [red, green])

  

Translate the line to [x0,y0]. This gives the tangent line to the original conic

 

tanline := eval(tanline1, [x = x-x0, y = y-y0])

1.852714154-5.696938463*x+2.696938459*y

 (7)

eval(tanline, [x = x0, y = y0])

0.2e-8

 (8)

display(implicitplot([conic(x, y), tanline], x = 0 .. 2, y = -1 .. 1, colour = [red, blue]), pointplot([x0, y0]))

  

NULL

Download Tangent_to_Conic.mw

Is this what you want?I could not get th...

Ronan 1506
December 27 2021
0 4

complex_log_plot.mw

Is this what you want?

I could not get the document to upload.

NULL;

z := x + y*I;

Plts := seq(plot3d([Re(k*log(z)), Im(k*log(z))], x = -3 .. 3, y = -3 .. 3, labels = ["Re(z)", "Im(z)", " ln(z)"], size = [800, 800]), k = [1, 2, 4, 6]);

NULL;
Plts[1];
Plts[2];
Plts[3];
Plts[4];
plots:-display(Plts);
NULL;

This is a workaround I use...

Ronan 1506
December 22 2021
1 2

I came up with this earlier in the year. Don't know how well it works for you graphics.

At the time I was taking notes from you tube that were formated in Scientific Workplace.

Basically I put a fixed sized table into the document constrained to A4 width. It is not perfect

but can be useful.

restart

 

 

 

 

 

 

NULL
 

Download A4_format.mw

mtaylor doesn't seem to work for f(x,y)...

Ronan 1506
December 19 2021
2 1

I was typing my answer while @Kitonum answered. I had to work around the i2(x,y)  as shown. Using Maple 2021.2

restart

NULL

i2 := proc (x, y) options operator, arrow; -((1/2)*I)*(exp(I*x)*sin(x)/x-exp(I*y)*sin(y)/y)/(x-y) end proc

proc (x, y) options operator, arrow; -((1/2)*I)*(exp(I*x)*sin(x)/x-exp(I*y)*sin(y)/y)/(x-y) end proc

 (1)

``

i2a := i2(y, z)

-((1/2)*I)*(exp(I*y)*sin(y)/y-exp(I*z)*sin(z)/z)/(y-z)

 (2)

i2b := i2(y, x)

-((1/2)*I)*(exp(I*y)*sin(y)/y-exp(I*x)*sin(x)/x)/(y-x)

 (3)

i3_r := -(1/2*I)*(i2a-i2b)/(z-x)

-((1/2)*I)*(-((1/2)*I)*(exp(I*y)*sin(y)/y-exp(I*z)*sin(z)/z)/(y-z)+((1/2)*I)*(exp(I*y)*sin(y)/y-exp(I*x)*sin(x)/x)/(y-x))/(z-x)

 (4)

``

> 

mtaylor(i3_r, [x,y,z], 10);

1/6+((1/1260)*I)*(x^5+x^4*y+x^4*z+x^3*y^2+x^3*y*z+x^3*z^2+x^2*y^3+x^2*y^2*z+x^2*y*z^2+x^2*z^3+x*y^4+x*y^3*z+x*y^2*z^2+x*y*z^3+x*z^4+y^5+y^4*z+y^3*z^2+y^2*z^3+y*z^4+z^5)-((1/90)*I)*(x^3+x^2*y+x^2*z+x*y^2+x*y*z+x*z^2+y^3+y^2*z+y*z^2+z^3)-(1/5670)*x^4*y*z-(1/5670)*x^3*y^2*z-(1/5670)*x^3*y*z^2-(1/5670)*x^2*y^3*z-(1/5670)*x^2*y^2*z^2-(1/5670)*x^2*y*z^3-(1/5670)*x*y^4*z-(1/5670)*x*y^3*z^2-(1/5670)*x*y^2*z^3-(1/5670)*x*y*z^4+(1/315)*x^2*y*z+(1/315)*x*y^2*z+(1/315)*x*y*z^2+((1/12)*I)*(x+y+z)-(1/5670)*x^4*z^2-(1/5670)*x^3*z^3-(1/5670)*x^2*z^4+(1/315)*x^2*z^2-(1/5670)*x*z^5-(1/5670)*x^5*z-(1/30)*x^2+(1/315)*x^4-(1/5670)*x^6-(1/30)*z^2+(1/315)*z^4-(1/5670)*z^6+(1/315)*x*z^3+(1/315)*x^3*z-(1/30)*z*x+(1/315)*y*x^3+(1/315)*y^2*x^2+(1/315)*y^3*x-(1/5670)*y*x^5-(1/5670)*y^2*x^4-(1/5670)*y^3*x^3-(1/5670)*y^4*x^2-(1/5670)*y^5*x-(1/30)*y*x-(1/5670)*z*y^5-(1/5670)*z^2*y^4-(1/5670)*z^3*y^3-(1/5670)*z^4*y^2-(1/5670)*z^5*y-(1/30)*z*y+(1/315)*z*y^3+(1/315)*z^2*y^2+(1/315)*z^3*y-(1/30)*y^2+(1/315)*y^4-(1/5670)*y^6

 (5)
> 

eval( (5), [x = 0, y = 0, z = 0]);

1/6

 (6)

Download series_expansion.mw

A potentially useful previous answer...

Ronan 1506
December 18 2021
0 0

Interesting question. This link may help you.

how to copy a plot into zoomed plot - MaplePrimes

Edit:- Another possible useful link

Plot inside one more plot - MaplePrimes

use the "is" command in Maple...

Ronan 1506
December 12 2021
2 0

7 nines as in your question
is(257^(1/2) <= 79999999/10000000);
                             false

8 nines
is(257^(1/2) <= 799999999/10000000);
                              true

 

Try this...

Ronan 1506
December 06 2021
1 1 
restart;
test := proc(x) local a; a := x^2; end proc;
test(5);
                               25

test(N*M);
                              2  2
                             N  M

I presume this is what you are trying to achieve.

This way....

Ronan 1506
December 06 2021
1 0

This shows how to get the answer.

Download Min_Length.mw

Look at this answer....

Ronan 1506
December 06 2021
1 0

There is an answer here. How to use a NULL result from fsolve as a condition to break a loop? - MaplePrimes

Is this what you want?...

Ronan 1506
December 05 2021
1 1

RMB on the plot and in the list of oprions at the bottom is export. Click and export as GIF. Then put the GIF file in powerpoint.

Powerpoint_GIF.pptx

3 4 5 6 7 8 9 Page 5 of 9