From cbdc205c3e7b5a3e8d54c39e3e393563229cac60 Mon Sep 17 00:00:00 2001
From: Kapil <kk733@cornell.edu>
Date: Sun, 12 Jan 2025 11:28:22 -0500
Subject: [PATCH] Update README.md

---
 README.md | 11 ++++++++++-
 1 file changed, 10 insertions(+), 1 deletion(-)

diff --git a/README.md b/README.md
index 76e72aa..e1c770f 100644
--- a/README.md
+++ b/README.md
@@ -80,8 +80,17 @@ cptmesh = cpt.mesh_sphere(name="sphere", radius=radius, center=(0, 0, 0), resolu
 cptmesh.keep_immersed_part(inplace=true)
 
 # declare it Julia mesh
-mesh = Mesh(cptmesh)
+mesh = Mesh(cptmesh)  
 ω = 1.03
 ζ = [0,0,1] # HEAVE: will be more verbose in future iteration. define it again even if defined in Capytaine.
 F = DiffractionForce(mesh,ω,ζ)
 A,B = calculate_radiation_forces(mesh,ζ,ω)
+```
+
+5. **Differentiability** :
+For differentiability with respect to mesh dimension, use `paper/MeshGradients_singlebody.jl`
+Differentiability needs an AD engine: use Zygote
+```julia
+using Zygote
+A_w_grad, = Zygote.gradient(w -> calculate_radiation_forces(mesh,dof,w)[1],ω)
+```